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FOREWORD 


General  Dynamics,  Convair  Aerospace  Division,  has  been 
involved  in  a  program  to  provide  technical  support  leading  to 
a  composite  wing  for  an  aircraft  transonic  improvement  program. 
The  data  derived  from  this  program  provide  an  aeroelastic 
design  of  a  complex  stiffness-critical  wing,  the  evaluation  of 
advanced  preliminary  scatLc  and  dynamic  wing  design/analysis 
procedures,  and  evaluation  of  the  reliability  characteristics 
of  critical  structural  components  utilizing  composite  materials. 

The  work  was  performed  under  Air  Force  Contract  F33615-70- 
C-1242,  "Composite  Wing  for  Transonic  Improvement  Program  - 
Aeroelastic  Analysis  and  Reliability  Program,"  under  the  spon¬ 
sorship  of  the  Air  Force  Flight  Dynamics  Laboratory,  Structures 
Division,  Air  Force  Systems  Command,  Wright- Pat ter son  Air  Force 
Base,  Ohio.  The  work  was  initiated  by  Mr.P.  A.  Parmley.  Mr. 

R.  T.  Achard  was  the  Air  Force  Project  Engineer,  and  Mr.  M.  E. 
Waddoups  was  the  program  manager  for  General  Dynamics. 

This  final  report  is  divided  into  three  volumes  as  follows: 

o  Volume  I  -  Composite  Wing  Aeroelastic  Response  Study 

o  Volume  II  -  Advanced  Analysis  Evaluation 

o  Volume  I?I  -  Structural  Reliability  Study 

The  work  reported  in  this  volume.  Volume  III,  concerns 
che  reliability  of  a  bonded  joint.  The  design,  analysis ,  and 
test  methods  have  been  developed  and  evaluated.  The  significant 
parameters  affecting  the  reliability  of  a  bonded  joint  are 
presented. 


This  technical  report  has  been  reviewed  and  is  approved. 


)/t/& 


Philip  A.  Parmley 
Chief,  Advanced  C^tfiposites  Branch 
Structures  Division 
Air  Force  Flight  Dynamics  Laboratory 
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SECTION  I 


INTRODUCTION 

1 . 1  OBJECTIVE 


The  objectives  of  this  program  were  (1)  to  develop  a  relia¬ 
bility-based  technology  for  the  design  of  large-scale  bonded 
joints  and  (2)  to  develop  a  data  base,  using  boron-epoxy-to-titan- 
ium  double-scarf  bonded  joint  specimens  to  explore  the  cnaracter 
of  the  failure  process  and  to  determine  the  joint  characteristics 
which  effect  structural  reliability. 


1.2  SCOPE  AND  METHODOLOGY 


Experimental  programs  in  composite  airframe  reliability  have 
followed  a  metals  technology-based  characterization  strategy. 

This  has  included  the  constant  amplitude  fatigue  (S-N)  character¬ 
ization  of  laboratory-sized  coupons,  spot  checks  of  the  effects 
of  specific  environments,  and  a  limited  check  of  the  effect  of 
classical  joint  design  parameters  such  as  f/t  and  adhesive  type 
upon  joint  S-N  lifetime.  This  approach  is  useful  in  screening 
operations,  but  it  is  not  satisfactory  for  the  design  and  optimi¬ 
zation  of  fatigue-critical  structures.  Therefore,  this  program 
was  conceived  to  investigate  the  following  critical  issuer: 

1.  The  load-history  sensitivity  of  large-scale  bonded 
joints  to  assess  the  mean  time  to  first  failure 

2.  The  implication  of  residual  strength  and  fatigu, 
lifetime  variability  on  the  determination  of  design 
safety  factors 

3.  The  character  of  the  damage  process  and  determination 
of  the  bounds  on  joint  reproducibility  that  would  be 
achievable  in  production  practice 

4.  The  determination  of  scale  effects  in  joint  design  and 
the  impact  of  scale  dependence  on  the  design  data  re¬ 
duction  scheme 

5.  The  eva’uation  of  the  random  load  history  characteriza¬ 
tion  technique  as  a  design  data  acquisition  tool. 
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A  central  issue  was  the  investigation  of  an  alternate  scheme 
for  the  service  life  design  characterization  of  structural  ele¬ 
ments.  The  problem  of  load-history  effects  on  the  lifetime  char¬ 
acteristics  of  a  structure  was  studied. 

An  experimental  program  was  conducted  using  1/5-  and  1/2- 
scale  specimens  ('both  5  inches  wide)  in  five  types  of  tests: 

(1)  strain  survey,  (2)  static,  (3)  constant  amplitude  (tension- 
tension  and  tension-compression),  (4)  residual  strength,  and  (5) 
lifetime  (fatigue  loaded  to  failure).  The  first  three  types  of 
tests  were  performed  using  only  1/5-scale  specimens.  A  strain 
survey  was  conducted  using  1/5 -scale  specimens  to  evaluate  the 
effects  of  specimen  dimensional  variations.  The  extent  of  these 
effects  and  their  impact  on  the  test  program  were  evaluated. 

Constant  amplitude  lifetime  tests  were  conducted  using  1/5- 
scale  specimens.  For  these  tests,  the  specimens  were  subjected  to 
a  prescribed  constant  amplitude  cyclic  loading  until  the  specimen 
failed.  The  purpose  of  these  tests  was  to  provide  a  basis  for 
evaluating  the  effectiveness  of  Miner's  linear  cumulative  damage 
rule  for  predicting  the  fatigue  life  of  bonded  joints. 

Residual  strength  tests  were  performed  using  1/5-  and  1/2- 
scale  specimens  that  had  been  subjected  to  a  prescribed  random 
load  history.  A  cycle-by-cycle  random  load-history  simulation 
for  one  lifetime  (1334  missions)  was  generated  using  exceedance 
data  projected  for  the  transonic  wing. 

This  load-history  data  was  stored  on  a  magnetic  tape.  A 
computerized  testing  system  read  the  magnetic  tape  and  transformed 
the  tape  input  into  a  closed  loop  controlled  random  load  on  the 
specimen.  Data  were  obtained  from  these  tests  for  characterizing 
the  lifetime  strength  of  compos ite-to-metal  bonded  joints,  of  the 
type  studied,  subjected  to  a  random  load  history. 

A  computer  program  was  developed  for  generating  a  random  load- 
history  simulation  on  a  cycle-by-cycle  basis;  this  simulation  pre¬ 
serves  the  cumulative  exceedance  statistics  for  a  given  fatigue 
spectrum.  This  program  was  based  on  the  techniques  implemented 
by  Dr.  J.  C.  Halpin  of  the  Air  Force  Materials  Laboratory  (Refer¬ 
ence  1). 

A  reliability-based  program  involves  the  determination  of 
variability  and  mean  lifetime  response  of  the  specimens.  Efforts 
to  understand  the  impact  of  fatigue  variability  upon  structural 
design  have  gained  appreciable  momentum  (References  2,  3,  4,  and 
5)  during  recent  years.  A  review  of  the  cited  efforts  provided  the 
initial  impetus  to  this  program.  A  primary  objective  in  the  program 
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was  to  provide  an  engineering  understanding  of  the  source  of  scat¬ 
ter  in  fatigue  lifetime  data.  Since  the  generation  of  fatigue 
failures  appears  to  be  a  random  process,  the  method  of  attack  for 
this  program  was  a  generation  of  statistical  data  and  analysis  of 
the  results. 

A  reliability  concepts  approach  was  used,  If  it  can  be  shown 
that  the  reliability  approach  is  feasible  for  characterizing  the 
strength  of  this  joint  type,  it  may  be  possible  to  develop  a  sound 
methodology  for  composite  bonded  joints  in  general.  The  experi¬ 
mental  and  analytical  procedure  used  to  achieve  the  stated  objec¬ 
tive  are  documented  in  this  report. 


1.3  IMPORTANCE  OF  JOINTS 


The  service  life  of  an  aircraft  is  intimately  dependent  on 
the  reliability  of  its  joints.  Experience  indicates  that  struc¬ 
tural  fatigue  failures  frequently  start  at  joints.  A  joint  failure 
in  primary  aircraft  structure  must  be  avoided  because  of  the 
possible  catastrophic  results.  Joints  traditionally  store  large 
amounts  of  energy  per  volume  of  material  compared  to  the  structure 
being  joined  and  exhibit  brittle  fracture  characteristics.  If  the 
energy  stored  in  a  joint  is  suddenly  released,  due  to  failure  of 
the  joint,  the  surrounding  structure  will  also  fail  unless  it  is 
able  to  dissipate  this  energy.  Therefore,  reliable  joints  are 
needed  to  maintain  the  structural  integrity  of  the  aircraft. 


1.4  FATIGUE-LIFE  PREDICTIONS 


Fatigue  is  a  random  phenomenon  that  is  not  fully  understood; 
nevertheless,  aircraft  designers  must  address  this  problem  to 
develop  a  satisfactory  structural  design.  The  designer  must  per¬ 
form  the  necessary  analysis,  tests,  and  evaluations  to  design  the 
structure  for  stiffness  and  for  static  and  fatigue  strength.  In 
addition,  the  designer  must  periodically  evaluate  the  structural 
integrity  of  the  aircraft  for  continued  usage.  This  task  is  of 
utmost  importance  for  the  safe  operation  of  the  aircraft;  there¬ 
fore,  the  fatigue  life  evaluation  must  be  as  realistic  as  possible. 

Safety  and  economics  are  both  important  factors.  From  a 
practical  standpoint,  the  fatigue-life  assessment  should  be  con¬ 
servative  for  safety  reasons  but  not  too  conservative  for  economic 
reasons.  If  the  fatigue-life  prediction  is  too  unconservative, 
the  maintenance  cost,  for  a  given  period  of  aircraft  usage,  may 
be  prohibitive. 
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It  is  imperative  that  the  fatigue  methodology  used  results 
in  realistic  service-life  predictions.  This  requires  that  the 
simulated  environment  characterize  actual  service  histories  as 
closely  as  possible.  Random  load  fatigue  testing  using  a  closed 
loop  test  system  and  a  simulated  environment  characterize  actual 
service  histories  as  closely  as  possible.  Random  load  fatigue 
testing  using  a  closed  loop  test  system  and  a  simulated  environ¬ 
ment  should  more  realistically  characterize  the  fatigue  process 
than  constant  amplitude  loading  methods.  Thus,  the  primary 
emphasis  in  this  test  program  was  to  evaluate  the  lifetime  and 
residual  strength  characteristics  of  the  bonded  joints  under  sim¬ 
ulated  service  loads  and  to  study  the  feasibility  of  using  random 
load  testing  in  the  fatigue-life  characterization  of  bonded  joints. 
Material  characterization  parameters  that  correlate  both  mean  and 
variance  for  static  strength  and  lifetime  were  determined. 

The  random  load  simulation  procedure  used  for  this  program 
is  described  in  Appendix  I. 


1.5  RANDOM  FATIGUE  PARAMETERS 


Structural  fatigue  under  realistic  random  loading  is  charac¬ 
terized  by  parameters  that  vary  during  the  service  lifetime.  The 
principal  parameters  that  vary  with  time  include: 

1.  Material  properties 

2.  Operating  loads 

3.  Operating  environments 

4.  Structural  strength. 

In  addition,  initial  factors  relating  to  structural  reliability 
exist,  such  as  manufacturing  deviations.  The  realistic  charac¬ 
terization  of  the  fatigue  process  in  composites  should  include  at 
least  the  random  variables  above. 
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1.6  LOAD-HISTORY  EFFECTS 

1.6.1  Metals 


The  load-sequence  effects  on  the  fatigue  life  of  metal  struc¬ 
tures  have  been  investigated  by  several  researchers  (References  6 
through  12).  The  load  history  simulation  for  a  given  service  load 
spectrum  is  important  for  realistically  characterizing  the  .fatigue 
process.  The  order  and  magnitude  of  the  fatigue  loading  have  a 
significant  effect  on  the  fatigue  life,  the  flaw  growth,  and  the 
residual  strength  of  the  structure.  The  sequence  effect  of  vari¬ 
ous  block  loadings  (simulating  the  same  load  spectrum)  on  the 
fatigue  life  of  2024  and  7075  aluminum  is  shown  in  Figure  1  (Page 
47  of  Reference  13).  The  low-high  sequence  produced  greater 
fatigue  damage  than  any  of  the  other  loading  sequences.  The  random 
block  sequence  was  the  next  most  damaging  load  sequence. 

Swanson  (Reference  14)  observed  that  cracks  are  more  likely 
to  occur  under  a  randomized  loading  than  a  deterministic  loading; 
he  also  observed  that  the  cracks  found  in  service  are  more  likely 
to  propagate  in  a  random  load  test  than  in  a  block  loading  test. 

The  literature  also  suggests  that  the  fatigue  damage  will  be  less 
representative  of  actual  service  failures  when  there  is  a  greater 
departure  from  the  true  randomness  in  loading. 

The  rate  of  crack  propagation  in  2024-T3  aluminum  is  dependent 
upon  the  loading  sequence  used  (Reference  15).  Crack  growth  is 
slower  for  a  high-low  sequence  than  for  a  low-high  sequence.  This 
difference  in  crack  growth  is  attributed  to  the  retardation  phenom¬ 
enon  occurring  in  the  crack  tip  areas.  That  is,  a  plastic  zone 
is  formed  by  high  tensile  stresses  in  the  crack  tip  region.  When 
the  tensile  stress  is  removed  or  reduced,  a  residual  compressive 
stress  field  is  left  in  the  plastic  zone.  This  residual  stress 
field  preloads  the  crack  and  tends  to  retard  the  rate  of  crack 
growth.  Compressive  stresses  may  be  produced  at  the  crack  tip  by 
gust  and  landing  conditions.  Such  stresses  tend  to  wipe  out  or 
reduce  the  crack  preload  which  in  turn  accelerates  the  rate  of 
crack  propagation.  The  progressive  creation  and  dissipation  of 
the  plastic  zone  continues  as  the  crack  propagates. 

The  retardation  phenomenon  has  been  investigated  for  D6ac 
steel  and  this  investigation  is  discussed  in  References  16  and  17. 
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1.6.2  Composites 


The  order  and  magnitude  of  the  fatigue  loading  can  significantly 
affect  the  resulting  fatigue  life  for  metals  (References  6  through 
12).  A  vigorous  development  program  vas  required  to  determine  the 
significance  of  the  load  sequence  for  metals. 

The  trend  in  fatigue  testing  is  to  use  a  random  loading  that 
simulates  the  service  history.  This  allows  the  load  sequence  effect 
to  be  accounted  for  in  the  laboratory. 

Fatigue  data  exists  for  composites  (References  18  through  25), 
but  this  data  is  inadequate  to  evaluate  the  significance  of  the 
load  sequence  on  the  fatigue  life  of  composites  in  the  same  way 
as  metal.  There  is  no  need  to  investigate  the  load  sequence  effect 
for  composites  if  random  load  fatigue  testing  is  used  because  the 
effect  will  be  reflected  in  the  test  results. 


1.7  FLY-IT-IN-THE-LAB  APPROACH 


The  desire  for  realism  in  fatigue  testing  can  be  accomplished 
by  the  "fly-it-in-the-lab  approach.”  This  approach  concerns  the 
random  simulation  of  the  service  load-environment  and  the  applica¬ 
tion  of  this  simulation  to  a  test  specimen,  component,  or  assembly 
using  a  computerized  closed  loop  test  system. 

A  random  load-history  simulation  procedure  (Reference  26)  has 
been  developed  to  preserve  the  composite  cumulative  lead  exceed¬ 
ances  and  the  waveform  defined  by  N0/Np.  The  environment  simu¬ 
lation  has  not  been  developed  into  an  integrated  random  load  test 
system;  however,  such  a  simulation  appears  technically  feasible. 

If  a  random  load  history  and  compatible  environment  can  be 
simulated  and  adapted  to  a  closed  loop  test  system,  this  methodology 
will  permit  a  practical  characterization  of  the  fatigue  process  in 
the  laboratory.  The  load-environment  simulation  can  be  generated 
by  mission  or  by  service  life.  The  degree  of  realism  depends  upon 
the  accurate  definition  of  the  expected  load-environment  history 
and  the  resulting  simulation  of  this  history. 
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1 . 8  REPORT  ORGANIZATION 


This  report  is  organized  into  eight  major  sections,  seven 
appendices,  and  a  list  of  references.  The  test  plan  and  objectives 
are  defined  in  Section  II. 

The  specimen  design,  fabrication,  and  inspection  details  are 
discussed  in  Section  III. 

Random  fatigue  testing  is  discussed  in  Section  IV.  This  in¬ 
cludes  the  documentation  of  the  random  load-history  simulation 
procedure  used  to  generate  the  test  tape  employed  in  the  experi¬ 
mental  program.  The  test  tape  statistics,  load  calibration,  and 
truncation  load  levels  are  also  discussed. 

The  execution  of  the  test  program  and  the  test  results  gene¬ 
rated  are  discussed  in  Section  V.  The  test  results  of  Section  V 
are  evaluated  in  Section  VI.  This  section  concerns  the  effective¬ 
ness  of  Miner's  rule,  the  statistical  evaluation  of  the  test  re¬ 
sults,  scaling  factors,  sample  size,  failure  modes,  effect  of 
residual  strength, and  fatigue  life  variability  on  required  design 
safety  factors,  etc. 

The  reliability  aspects  of  bonded  joints  are  discussed  in 
Section  VII,  and  the  principles  described  are  used  to  analyze  the 
test  results  of  Section  V.  Wearout  models,  residual  strength, 
Weibull  distribution,  sample  size,  etc.  are  discussed. 

The  conclusions  and  recommendations  of  this  investigation  are 
given  in  Section  VIII. 

Section  VIII  is  followed  by  seven  appendices  and  a  list  of 
references.  The  appendices  contain  the  details  of  the  spectrum 
simulation  technique  used,  the  specimen  warpage  problem,  a  glossary 
of  technical  terms,  photographs  of  failed  specimens,  statistical 
properties  of  the  Weibull  distribution,  and  a  discussion  of  a 
"least-of-two"  statistical  analysis. 
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SECTION  II 
TEST  PLAN 


2.1  TEST  OBJECTIVE  j 

The  objective  of  the  experimental  program  was  to  generate  ] 

statistical  data  characterizing  the  residual  strength/lifetime  ■ 

characteristics  of  boron  laminate-to-titanium  double-scarf  adhes-  ! 

ive  bonded  joints  using  simulated  flight-by-flight  loading. 

The  experimental  program  is  summarized  in  Table  I.  Two  sizes 
of  bonded  specimens  were  used.  These  were  1/5-  and  1/2-  scale 
boron  laminate-to~titanium  double-scarf  adhesive  bonded  joints. 

Specimen  identifications  and  test  assignments  are  summarized  in 
Tables  II  and  III  for  1/5-  and  1/2-scale  specimens, respectively . 

The  F  and  K  prefixes  on  the  specimen  numbers  in  Tables  II  and  III 
refer  to  manufacturing  period.  The  K  specimens  were  manufactured 
before  May  1971  and  the  F  specimens  were  manufactured  after 
October  1971. 


2.2  TEST  SYSTEMS 


A  closed  loop  test  system,  shown  in  Figure  2,  was  used  for 
random  load  tests.  A  Baldwin-Tate-Emery  Universal  test  machine 
was  used  for  static  and  residual  strength  tests  (Figure  3).  Photo¬ 
graphs  of  the  system  are  shown  in  Figures  3  and  4  and  a  description 
of  the  apparatus  used  is  given  in  the  following  paragraphs.  The 
essential  elements  of  this  system  are: 

1.  Magnetic  tape 

2.  Varian  620/ i  computer 

3.  Servo  amplifier 

4.  Servo  valve 

5.  Hydraulic  ram 

6.  Load  cell 

7.  Test  fixture  (3  for  1/5-scale  specimens  and  1  for  1/2- 
scale  specimens) 
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Table  I  EXPERIMENTAL  TEST  PROGRAM 
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Notes:  ^hese  speciraen  also  used  for  static  and  constant  amplitude  tests. 

Total  number  of  specimens  for  test  program  -  81  (1/5  scale)  10  (1/2  scale) 
All  tests  were  conducted  in  a  normal  laboratory  environment. 
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Table  II  SUMMARY  OF  1/5-SCALE  SFECIMENS 
AND  TEST  ASSIGNMENTS 


Specimen  I.D. 

Test  Assignment 

F504417 

FTF 

F504418 

FTF 

F504419 

FTF 

F504421 

FTF 

F504622 

FTF 

K012437* 

S 

K012438 

CA(T-C) 

K0124394 

S 

K012440* 

s 

K0124414 

s 

K012442 

CA(T-C) 

K012443 

FTF 

K012444 

CA(T-T) 

K012445 

CA(T-C) 

K012446* 

s  ; 

K012459* 

S 

K012461 

CA(T-T) 

K012462 

CA(T-C) 

K012463* 

S 

K012464 

FTF 

K012465 

CA(T-T) 

KOI 2474* 

S 

K012475 

CA(T-C) 

K012476 

CA(T-T) 

K012477 

CA(T-T) 

K012478 

CA(T-T) 

K012479 

FTF 

K012480 

FTF 

K012481 

FTF 

K012482 

CA(T-T) 

K012483 

CA(T-T) 

K012484 

FTF 

K012485 

FTF 

K900441 

R(17.) 

K900442 

CA(T-C) 

K900443 

FTF 

K900445 

CA(T-T) 

K900446 

FTF 

K900447 

R(10%) 

K900448 

R(17o) 
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Table  II  SUMMARY  OF  1/5-SCALE  SPECIMENS 

AND  TEST  ASSIGNMENTS  (Continued) 

Specimen  I.  D. 

Test  Assignment 

K900449 

R(50%) 

K900450 

R(50%) 

K904950 

R(10%) 

K904951 

R(50%) 

K904952 

R(50%) 

K904953 

R(50%) 

K904954 

R(5C%) 

K904955 

R(l%) 

K904956 

R(50%) 

K904957 

R(10%) 

K904953 

R(l%) 

K904959 

R(10%) 

K905015 

R(507o) 

K905016 

R(17o) 

K905017 

R(107o) 

K905018 

R(507o) 

K905020 

R(50%) 

K905021 

R(107o) 

K905022 

R(507o) 

K905023 

R(507c) 

K905362 

R(507o) 

K905363 

R(507o) 

K905364 

FTF 

K905365 

R(507o) 

K905366 

R(507o) 

K905367 

FTF 

K905368 

FTF 

K905369 

R(507.) 

K905371 

R(507o) 

K905467 

R(50%) 

K905468 

FTF 

K905469 

R(507o) 

K905470 

FTF 

K994565 

FTF 

K994566* 

CA(T-T) 

K994567 

CA(T-T) 

K9945684 

CA(T-C) 

K994569* 

S 

K994578 

S 

K994587 

CA(T-T) 

K994588n 

S 
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Table  II  (Concluded) 


Notes 

S 

CA(T- 

CA(T- 

R(X7o) 

FTF 

* 

A 


Notes : 


=  Static  test  to  failure 

T)  =  Constant  amplitude  (Tension-Tension)  test 
to  failure 

C)  =  Constant  amplitude  (Tension-Compression)  test 
to  failure 

-  Random  Loading  (X%  Specified  lifetime)  and 
static  test 

=  Fatigue  to  failure 

=  Specimen  rejected  due  to  excessive  dimensional 
variations 

=  Specimen  also  and  for  strain  survey 


Table  III  SUMMARY  OF  1/2-SCALE  SPECIMEN 
AND  TEST  ASSIGNMENTS 


Specimen 

Test  Assignment 

F504619 

FTF 

F504620 

FTF 

F504621 

FTF 

F504622 

FTF 

F504623 

FTF 

F504618 

R(10%) 

F504624 

R(107o) 

F504625 

R(10%) 

K012453 

R(107o) 

K012454 

R(107o) 

R(X%)  =  Random  Loading  (X%  lifetime)  and  static  test 

FTF  =  Fatigue  to  failure 
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Figure  3  Closed  Loop  Test  System  (Test  Tape,  Varian 
620/i  Computer,  Servo  Amplifiers) 


Figure  4  Closed  Loop  Test  System  (Servo  Valve,  Hydraulic  Ram, 
Load  Cell,  Test  Fixture  and  Test  Specimen) 
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8.  Test  specimen 

9.  Data  recording  equipment 

10.  Strain  gauges. 

A  random  load  history,  simulating  a  given  service  spectrum, 
is  generated  and  stored  on  magnetic  tape.  This  tape  is  analyzed 
by  a  Varian  620/i  digital/analog  converter.  The  digital  records 
on  the  tape  are  converted  into  an  analog  signal  that  is  sent  to 
the  servo  amplifier. 

A  servo  valve  controls  the  flow  of  oil  for  driving  the 
hydraulic  ram.  The  ram  puts  out  an  arbitrary  load  signal.  This 
signal  goes  through  the  load  cell  back  to  the  servo  amplifier. 

The  signal  from  the  load  cell  is  analyzed  and  adjusted  to  meet 
the  magnetic  tape  requirement  before  the  load  is  applied  to  the 
specimen.  The  servo  amplifier  compares  the  signal  from  the  load 
cell  and  the  random  load  tape.  A  compensating  load  signal  is 
sent  directly  to  the  servo  valve  from  the  servo  amplifier.  This 
signal  is  compatible  with  the  load  requirement  from  the  test  tape. 
The  servo  valve  adjusts  the  oil  flow  to  the  hydraulic  ram,  which 
in  turns  loads  the  load  cell  and  the  specimen.  After  each  load 
is  applied  to  the  specimen,  a  new  signal  is  sent  to  the  servo 
amplifier  and  r  spared  with  the  signal  from  the  magnetic  tape. 

This  process  continues  until  the  desired  number  of  loads  are 
applied  to  the  specimen  or  unui  the  specimen  fails.  The  closed 
loop  test  system  is  used  for  both  the  constant  amplitude  tests 
and  the  random  tests.  The  Varian  620/i  digital/analog  converter 
is  used  only  for  the  random  fatigue  tests. 

Four  specimens  can  be  tested  simultaneously  using  the  closed 
loop  test  system.  Three  1/5-scale  specimen  and  one  1/2-scale 
specimen  tests  can  be  conducted  at  the  same  time.  This  is  possible 
using  a  servo  valve  and  servo  amplifier  for  each  test  fixture  used. 
The  loads  stored  on  the  magnetic  tape  are  in  nondimens ional  wave 
form.  The  desired  load  on  the  specimen  is  obtained  by  calibrating 
the  load  signal  output  for  each  servo  amplifier.  The  load  history 
for  each  specimen  is  recorded  and  the  load  signal  is  plotted. 

When  only  1/5-scale  specimens  are  tested,  a  test  rate  of 
about  5  Hz  is  possible;  however,  when  a  1/2-scaie  specimen  is 
also  tested,  the  test  rate  is  reduced  to  3  or  4  Hz  because  of 
system  limitations.  Since  the  recording  equipment  is  independent, 
different  lifetime  random  tests  can  be  run  at  the  same  time. 
Specimens  are  tested  one  at  a  time  or  in  multiples  until  the  test 
requirements  are  satisfied. 
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SECTION 


III 


f; 

i 


SPECIMEN  DESIGN,  FABRICATION, 
AND  INSPECTION 

3.1  SPECIMEN  DESIGN  DETAILS 


The  specimen  descriptions,  design  criteria,  quality  control 
acceptance  tests,  and  material  mechanical  properties  are  dis¬ 
cussed  in  this  section. 


3.1.1  Description 


Boron-laminate-to-titanium  adhesive  bonded  double  scarf  joint 
specimens  (1/5  and  1/2  scale)  were  used  in  the  experimental  test 
program.  Specimen  constituents  include 

1.  Two  titanium  (6AL-4V)  lug  ends  with  a  double  tapered 
splice  plate 

2.  Bo^on-epoxy  laminate  (Narmco  5505),  pre-cut  to  match 
template 

3.  General  Dynamics  Specification  FMS-1013  IA  adhesive 
(Relibond  398,  Reliable  Manufacturing  Company) 

4.  General  Dynamics  Specification  FMS-1014  IIA  adhesive 
primer. 

Drawings  of  the  1/5-  and  1/2-scale  specimens  are  shown  in  Figure 
5  and  6, respectively. 

The  boron-epoxy  laminate  was  bonded  to  the  titanium  lug  ends 
with  FMS-1013  IA  adhesive  to  FW6920025  and  FW6920044  specification. 
The  fabrication  and  inspection  techniques  are  discussed  in  subsec¬ 
tion  3.2  and  3.3,  respectively. 

The  FW6920025  (1/5  scale)  and  the  FW6920044  (1/2  scale)  speci¬ 
mens  had  the  same  f/t  value  based  on  the  average  taper  thickness 
of  the  boron-epoxy  laminate.  The  bond  areas  and  the  cross-sec¬ 
tional  areas  through  the  laminate  for  the  respective  specimens 
(FW6920025  and  FW6920094)  were  proportioned  so  that  the  smaller 
specimen  was  approximately  1/5-scale  size  relative  to  the  1/2-scale 
specimen.  \-j 


Pre-Cut  Detail 


Figure  5  One-Fifth  Scale  Fatigue  Specimen  -  Boron -to -Titanium 
Double  Scarf  Adhesive  Bonded  Joint 
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Pre-Cut  Detail  (2) 
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Pre-Preg  (Total) 


1.50R  Typ 


Figure  6  One-Half  Scale  Fatigue  Specimen  -  Boron-to-Titanium 
Double  Scarf  Adhesive  Bonded  Joint 


3.1.2  Design  Criteria 


The  design  development  tests  leading  to  the  joint  configura¬ 
tions  used  in  this  program  are  documented  in  Reference  27.  The 
original  design  goal  was  to  develop  a  full-scale  bonded  joint  with 
the  capacity  to  transmit  66,000  Ibs/in  at  ultimate  load.  The 
reproducibility  characteristics  of  the  bonded  joints  were  not 
known  before  this  program;  therefore,  joint  allowables  were  estab¬ 
lished  by  test.  Special  factors  of  safety  were  used  in  the 
Reference  27  program;  hence,  ultimate  load  was  equal  to  1.5  times 
the  limit  load.  Laminate  design  allowables  based  on  the  "B" 
criteria  of  MIL-HDBK  5  were  used  to  determine  allowable  stress 
levels . 

The  full-scale  test  results  of  Reference  27  indicated  a  pos¬ 
sible  scale  effect.  These  results  were  used  to  develop  1/5-  and 
1/2-scale  fatigue  specimens  for  this  program  compatible  with  the 
full-scale  specimens  of  Reference  27. 

The  dimensions  of  the  scaled  specimens  are  listed  in  Table 
IV.  The  1/5-  and  1/2-scale  specimens  have  the  same  ratio  of 
laminate  thicknesses  as  the  ratio  of  scarf  length.  This  is  re¬ 
flected  in  the  i/t  ratio. 


3.1.3  Acceptance  Tests 


The  material  used  for  this  program  was  Narmco  5505  boron- 
epoxy  tape  per  General  Dynamics  specification  FMS-2001A  manufac¬ 
tured  by  Narmco  Materials  Division  of  the  Whittaker  Corporation. 
The  acceptance  data  obtained  at  General  Dynamics  for  sample 
batches  of  material  used  in  this  program  are  shown  in  Table  V 
and  compared  with  the  requirements  of  FMS-2001A.  Data  in  Table  V 
are  based  on  the  ultimate  strength  averages  of  three  specimens 
per  point.  Statistical  properties  for  the  data  of  Table  V  are 
summarized  in  Table  VI. 


3.1. A  Material  Mechanical  Properties 


Room  temperature  physical  properties  for  boron-reinforced- 
epoxy  lamina  are  summarized  in  Table  VII.  Material  correction 
factors  for  temperature  of  -65°F  to  A20°F  are  also  available. 
Design  allowables  are  summarized  in  Table  VIII  for  boron- 
reinforced-epoxy  composites  at  room  temperature.  A  design 
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Table  V  SUMMARY  OF  ACCEPTANCE  TEST  RESULTS  FOR 
BORON  COMPOSITE  FIBERS® 


BATCH 

ROLL 

DATE 

TEMP. 

TRIAL 

b 

d 

NO. 

NO. 

NO. 

TESTED 

NO. 

r(psi) 

0°(P8l> 

o-90o(psi) 

1 

Eim 

vmnmm 

EQEBKK 

R.T. 

2 

ESS 

I  I 

3 

16.3 

EBoKi 

81490 

401 

5 

8-17-70 

Avr. 

16.0 

254.6 

14.7  1 

4 

gngB 

ElffKK 

350°F 

5 

■m! 

QQ|K 

6 

KK 

11.8 

Avr  . 

6.7 

222.8 

12.2 

1 

251.4 

16.5 

R.T. 

2 

15.8 

234.2 

16.6 

3 

16.2 

218.6 

17.6 

81572 

401 

14 

9-1-70 

Avr. 

15.0 

234.7 

16.9 

4 

6.5 

195.1 

11.0 

350°F 

5 

6.8 

212.0 

12.7 

6 

6.7 

12.7 

Mznm 

203.5 

12.1 

1 

M*wim 

Br-ft  1 * 

EQOH 

R.T. 

2 

ilfl 

ES^H 

3 

IsMi 

215.5 

EkkKK 

81507 

402 

Sample 

8-17-70 

Avr. 

14.7 

226.6 

15.4 

4 

8.5 

216.2 

12.5 

350°F 

5 

8.2 

208.3 

12.4 

6 

7.4 

213.8 

12.3 

. 

__Avr._ 

8.0 

212.7 

12 A _ _ 

1 

■grip 

inm 

R.T. 

2 

1  1 

3  , 

IM 

15.9 

81618 

402 

28 

9-1-70 

msm 

17.0 

277. 2 

16.7 

4 

mum 

WWWMM 

ign| 

350°F 

5 

■ 

ISIS 

E 

6 

KB 

232.4 

BOuKK 

Avr. 

7.0 

231.6 

11.9  1 

1 

m Tsmwm 

EgnB 

R.T. 

2 

B 

SK 

3 

16.5 

E&kKI 

BH9H 

81571 

402 

32 

9-1-70 

Avr. 

15.9 

267.6 

16.1  J 

warn. 

EyQKB: 

II 

i 

ExE^H 

BEK 

wnwmm 

kSKKI 

mmm 

BBSBB 

221.9 

12,9  f 
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Table  V  SUMMARY  OF  ACCEPTANCE  TEST  RESULTS  FOR 
BORON  COMPOSITE  FIBERS® 

(Continued) 


T.T. 

NO. 

BATCH 

NO. 

ROLL 

NO. 

DATE 

TESTED 

TEMP. 

TRIAL 

NO. 

b 

T(psi) 

^(psi) 

<r90°(psi) 

■m 

1 

1F1« 

maam 

1 

2 

Elwm 

W- 

I 

1 

3 

im 

217.4 

Ehm 

749409 

402 

52 

9-18-70 

1 

Esm 

13.2 

248.0 

16.7 

4 

5.9 

179.0 

12.2 

350°F 

5 

6.1 

170.9 

11.2 

6 

5.8 

177.5 

10.9 

1S53HH 

5.9 

175.8 

EfKpaai 

Notes:  a  Acceptance  test  per  FMS-2001 

b  Horizontal  Shear  (Required  Values:  13  psi  at  R.T.  and 
5  psi  at  350°F) 

c  0°  Flexural  (Required  Values:  225  psi  at  R.T.  and  170  psi 
<a  350°F) 

d  90°  Flexural  (Required  Values:  13  psi  at  R.T.  and  8  psi  at 
350°F) 
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Table  VI  SUMMARY  OF  STATISTICAL  PROPERTIES  BASED  ON  ACCEPTANCE  TEST 
RESULTS  FOR  BORON  COMPOSITE  FIBERS 


Statistical 

Property 

r 

o 

o 

b 

<r$ 

10° 

R.T. 

350°F 

R.T. 

350°F 

R.T. 

350°F 

Mean  (psi) 

<^) 

15.29 

1 

6.81 

251.48 

211.42 

16.06 

12.15 

Variance 

( <r2) 

1.963 

0.442 

438.235 

357.919 

0.908 

H 

Standard 

Deviation 

(cr> 

1.401 

0.665 

20.93 

18.919 

J" 

0.953 

0.680 

Coefficient  of 
Variation 
( 

9.164% 

9.764% 

8.323% 

8.949% 

5.933% 

5.597% 

Notes:  Same  notes  apply  as  given  at  bottom  of  Table  V, 

t*  =  l  Xi 

i  =  1 


n 


0-=  v/£  (Xhu) 

i=  l 


O' 


V  =  (~)  100 


=  18  (Reference  data  Table  V). 


24 


n 


BggyqiT  tom  Viwry 


i 


» 

ft 


Table  VII  PHYSICAL  PROPERTIES  -  BORON-EPOXY 


Property 


i 


Density,  lb/in. 

Ply  Thickness,  in. 

Volume  of  Reinforcement,  % 

.  .  .  BTU  -  in. 

fhermal  Conductivity,  ~  ^2 - - 


Linear  Coefficient  of  Thfermal 
Expansion,  10“6  in./in./°F 


Dielectric  Constant,  8.2  GHz 

Loss  Tangent,  8.2  GHz 

Capacitance,  ft  ft  farad,  dry 

Capacitance,  ft  ft  farad,  humidity 
exposure 

Specific  Heat,  BTU/lb/°F 


Value 


0.073 

0.0051  to  0.0054 
50  (average) 

2.8  at  200°F^ 

(0°  Lamina) 

2.5  at  2 00°F paralle  1 
(0°  Lamina) 

13.0  at  200°F^2\  normal 
(0°  Lamina) 

27  at  R.t/3^ 

(0°  Lamina) 

0.007  at  R.T. ^ 

(0°  Lamina) 

221.5(4) 

(Unpainted  0°  Lamina) 
225.5(4) 

(Unpainted  0°  Lamina) 
0.327  at  200°F 


NOTES:  (1)  Data  available  for  temperature  range  -65F  to  420°F 
and  for  specific  laminates. 

(2)  Data  available  for  temperature  range  -50F  to  350°F 
and  for  specific  laminates. 

(3)  Data  available  for  temperature  range  R.T.  to  420°F 
and  frequency  range  8.2  to  11.0  GHz  and  for  specific 
laminates . 

(4)  Per  MIL-R-7705A(ASG) r  Data  also  available  for 
specific  laminates. 

(5)  From  Reference  28. 
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*  le  VIII  DESIGN  ALLOWABLES  -  BORON-EPOXY 
ROOM  TEMPERATURE,  "B"  BASIS  ’ 


Property 


_ _ Direction 

jarallel  to  Fiber  Itformal  fo 


Axial 

Modulus  of  Elasticity 

106  psi  y' 

30.0 

2.7 

Design  Limit  Tensile" 

Strain,  in. /in. 

0.0040 

0.0023 

Design  Ultimate  Tensile 
Strain,  in. /in. 

0.00600 

0.00345 

Design  Limit  Compressive 
Strain,  in. /in. 

0.0075 

0.0075 

Design  Ultimate  Compres¬ 
sive  Strain,  in. /in. 

0.01125 

0.01125 

Shear 

io6apSi°dulus  of  Ri§idity> 

0.65 

Poisson's  Ratio 

0.21 

Design  Limit  Shear  Strain 
in. /in.  ’ 

0.013 

■ 

Design  Ultimate  Shear 

Strain  in  /,•  _ 

0 .0195 

ultimate  intersection  diagram  is  shown  in  Figure  7  for  (O2/+45) 
laminate  configurations. 


3.2  FABRICATION  TECHNIQUE 


Reusable  titanium  lug  ends  were  used  to  build  the  1/5-  and 
1/2-scale  specimens.  Thirty  sets  of  the  1/5-scale  lug  ends  and 
three  sets  of  the  1/2-scale  lug  ends  were  used  in  the  program. 
Specimens  were  built,  tested,  and  rebuilt  until  the  required  test¬ 
ing  was  complete. 

The  fabrication-test  recycling  procedure  is  shown  in  Figure 
8.  Typical  pictures  of  the  specimen  elements,  tools,  and  appara¬ 
tus  are  shown  in  Figures  9  and  10. 

At  the  start  of  the  program,  tools  were  available  for  build¬ 
ing  five  1/5-scale  and  one  1/2-scale  specimens  at  the  same  time. 

As  the  testing  progressed,  it  became  apparent  that  a  faster  rate 
of  1/5-scale  specimen  recycling  was  needed  to  keep  pace  with  the 
testing.  Five  additional  tools  were  made  for  building  the  1/5- 
scale  specimen.  With  the  new  tools,  ten  1/5-scale  specimens  could 
be  built  at  the  same  time.  No  additional  tools  were  made  for 
building  the  1/2-scale  specimens  since  only  ten  specimens  were 
required  for  the  test  program. 


3.3  INSPECTION 


3.3.1  Fluorescent  and  Ultrasonic 


The  titanium  lug  ends  were  inspected  with  fluorescent  pene¬ 
trant  per  MIL-I-6866  type  1,  and  the  specimens  were  ultrasonic 
inspected.  Approximately  30  of  the  first  81  1/5-scale  specimens 
and  3  of  the  first  10  1/2-scale  specimens  were  inspected  as 
described.  It  was  found  that  the  rate  of  specimen  supply  could 
not  stay  up  with  the  needs  for  continuous  testing.  To  avoid 
testing  delays,  inspection  was  waived  on  the  remaining  specimens. 


3.3.2  Warpage 


The  first  25  1/5-scale  specimens  built  had  a  noticeable  de¬ 
gree  of  warpage.  The  degree  of  warpage  was  determined  by  placing 
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Figure  7  Boron-Epoxy  [C^/j:  45]  Design  Ultimate  Interaction  Diagram,  RT 


Figure  8  Fabrication-Test- Recycling  Procedure 
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the  specimen  on  a  solid  level  surface  and  taking  relative  dial  | 

gauge  surface  measurements  at  the  locations  shown  in  Figure  11.  | 

Relative  measurements  were  taken  with  respect  to  surface  location  ; 

2.  Surface  measurements  are  summarized  in  Table  IX  for  24  of  the 
25  specimens  measured.  i 

? 

Straightening  of  specimen  K012440  was  attempted  by  heating. 

Seven  of  the  24  specimens  evaluated  were  rejected  because  of  ex¬ 
cessive  warpage.  Specimen  K012440  was  also  rejected.  The  rejec¬ 
ted  specimens  were  static  tested  and  the  lug  ends  were  recycled. 

The  static  test  results  for  these  specimens  were  segregated  from 
the  acceptable  static  test  specimens  (Table  II). 

The  specimen  fabrication  procedure  was  analyzed  to  correct 
the  warpage  problem.  It  was  found  that  the  layers  of  bleeder 
cloth  between  the  top  surface  of  the  tool  and  the  bottom  surface 
of  the  specimen  prevented  the  lug  ends  from  resting  flat  against 
the  surface  of  the  tool.  This  problem  was  easily  corrected  by 
placing  cloth  plies  between  the  lug  ends  and  the  tool  to  compen¬ 
sate  for  the  bleeder  cloth.  All  specimens  built  after  this  change 
were  visibly  checked  for  flatness  by  placing  specimens  on  a  level 
surface  and  observing  clearances.  None  were  considered  unaccept¬ 
able.  All  1/2-scale  specimens  were  built  with  allowance  for  the 
bleeder  cloth  thickness;  therefore,  no  warpage  problems  occurred. 

The  impact  of  specimen  warpage  was  investigated  to  determine 
if  there  would  be  a  serious  effect  on  the  experimental  results. 

A  strain  survey  was  performed  to  evaluate  the  effects  of  specimen 
warpage.  The  strain  data  generated  and  the  analysis  performed 
are  given  in  Appendix  VII.  It  was  concluded  that  the  effects  of 
specimen  warpage  were  negligible.  This  is  documented  in  Appendix 
II. 
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Table  IX  SUMM/  "  DIAL  GAUGE  SURFACE  MEASUREMENTS 
FOR  INI  AL  GROUP  OF  1/5 -SCALE  SPECIMEN 


•s 

l 

[ 


I 


.i 

) 


k 


♦Specimen  rejected  due  to  excessive  dimensional  variations 
4Dial  gage  surface  locations  are  shown  in  Figure  11. 
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SECTION  IV 
RANDOM  FATIGUE  AND 
LOAD  SPECTRUM 


The  random  load  history  generation  procedures  used  for  this 
program  are  discussed  in  this  section.  These  procedures  were  used 
to  generate  the  test  tape  random  load  simulation  discussed  in  sub¬ 
section  4.3.  A  computer  program  (Reference  26)  was  developed  which 
simulates  the  random  load  history  and  stores  it  on  a  magnetic  tape. 

This  tape  was  used  in  a  computerized  closed  loop  test  system. 

Most  of  the  technical  terms  used  in  this  section  are  defined  in  j 

Appendix  III.  | 

j 

4.1  GENERAL  DESCRIPTION  OF  PROCEDURE  j 


The  principal  steps  and  essential  elements  of  the  random  load 
history  simulation  procedure  are 

1.  Define  Service  History 

a.  Mission  profiles 

b.  Mission  segment  types 

c.  Time  spent  in  each  segment 

d.  Number  of  missions 

e.  Exceedances  for  positive  and  negative  spectra 

2.  Analyze  exceedance  curves  using  method  by  Press  (Refer¬ 
ence  29) 

a.  Number  of  RMS  levels  for  each  mission  segment  type 

b.  a PMS  values 

c.  Clipping  ratios 

3.  Power  Spectrum  Analysis 


a.  Power  Spectral  Density  (PSD) 

b.  N0 

c.  N0 


4.  Random  Number  Mapping 

a.  Rice-Bendat-Kowalewski  probability  density  function 

b.  Cumulative  distribution  curve 

c.  Pseudo  random  number  generator 

d.  Acceptable  random  numbers 

5.  Load -by-Load  Creation 

a.  Use  mapping  procedure  repeatedly 

b.  Pj  =  Pmeanj  +  (-l)j+1  (fitted  random  no.)j  (^RMS)^ 

c.  Haversine  wave  form  connects  successive  loads 

d.  Bookkeeping 

6.  Statistical  Analysis  of  Random  Load  History  Simulation 

a.  Level  crossings 

b.  Mean  level  crossings  with  positive  slope  (N0) 

c.  Total  number  of  loads  for  simulation 

7.  Test  Tape  Generation. 


4.2  SERVICE  HISTORY  DESCRIPTION 


The  random  load  history  simulation  used  in  this  program  is 
based  on  the  F-lll  TIP  "B"  wing  spectrum , partially  documented  in 
Reference  30.  This  is  a  maneuver  load  dominated  spectrum  (as 
opposed  to  a  gust  critical  spectrum).  A  typical  mission  profile 
was  developed  from  this  data.  The  typical  mission  length  was 
assumed  to  be  3  hours.  For  a  4000-hour  airplane,  this  is  equiva¬ 
lent  to  1334  missions  per  service  lifetime.  It  was  assumed  that 
the  typical  mission  was  repeatedly  flown  1334  times  in  one  service 
life.  A  random  load  history  was  generated  to  simulate  the  loads 
occurring  in  one  service  life. 

A  typical  fighter  aircraft  flies  20  to  30  different  types  of 
missions.  This  program  was  more  concerned  with  studying  the  effect 
of  a  random  load  history  than  it  was  with  the  fact  that  a  single 
mission  was  used  to  represent  one  service  lifetime.  The  mission 
profile  and  the  corresponding  exceedance  curves  for  each  mission 
segment  type  are  discussed  in  paragraphs  4.2.1  and  4. 2. 2, re¬ 
spectively. 
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4.2.1  Mission  Profile 


A  mission  profile  describes  the  type  and  time  of  each  maneuver 
by  the  aircraft  in  a  single  mission.  Typical  maneuvers  are  climb, 
cruise,  terrain  following  radar  (TFR) ,  air-to-ground  cycles  (A/G), 
air-to-air  cycles  (A/A),  ground-to-air-to-ground  cycles  (GAG),  and 
touch-and-go  cycles  (T&G) .  A  maneuver  or  maneuvers  are  classified 
by  mission  segment  types.  Several  mission  segment  types  are  re¬ 
quired  to  characterize  a  single  mission.  Different  missions  are 
described  by  assembling  the  mission  segment  types  in  various 
orders . 

The  mission  profile  used  for  the  random  load-history  simula¬ 
tion  for  this  program  is  shown  in  Figure  12.  This  mission  profile 
represents  a  typical  F-lll  mission.  The  mission  segment  number, 
the  maneuver  description,  and  the  time  distribution  for  each  man¬ 
euver  is  summarized  in  Table  X. 


Table  X  SUMMARY  OF  MISSION  PROFILE  DATA 


Mission 

Segment 

No. 

Maneuver 

Description 

Time  Distribution 
Time 

%  Mission  (Hours) 

1 

Climb-Cruise 

30% 

0.90 

2 

TFR 

20% 

0.60 

3 

A/G 

15% 

0.45 

4 

TFR 

10% 

0.30 

5 

A/A 

5% 

0.15 

6 

Cruise,  descent,  T&G 

20% 

0.60 

Mean  wing 

root  B.M.  (lg)  =  2.1  x 

106  in  lb 
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Figure  12  Mission  Profile 


An  exceedance  curve  defines  the  number  of  times  per  unit  time  a 
given  delta  load  level  is  exceeded  for  a  given  mission  segment 
type.  The  exceedance  data  for  this  report  is  based  on  the  F-lll 
TIP  MB"  wing  spectrum.  Cumulative  exceedances  per  mission  were 
plotted  against  the  wing  root  (AB.M.)^  for  analyses  purposes. 

The  positive  spectrum  cumulative  exceedance  curves  used  to  gener¬ 
ate  the  test  tape  are  given  in  Figures  13  through  18.  For  this 
program,  the  random  load  history  used  in  the  test  was  tension; 
hence,  for  computational  convenience  it  was  assumed  that  the  nega¬ 
tive  spectra  was  the  same  as  the  positive  spectra. 

The  exceedance  curves  of  Figures  13  through  18  are  analyzed 
using  the  method  of  Press  (Reference  29).  Results  are  shown  on 
the  corresponding  figures.  Press's  method  is  described  in  Appen¬ 
dix  I.  A  summary  of  exceedance  curve  parameters  used  for  generat¬ 
ing  the  test  tape  are  given  in  Table  XI. 


4.3  LOAD  HISTORY  SIMULATION  FOR 
RANDOM  FATIGUE  TESTS 


This  section  describes  the  input  data  and  its  applies  ':ion  in 
the  development  of  the  random  load  history  tape  employed  in  the 
random  fatigue  tests.  The  random  load  simulation  is  statistically 
analyzed, and  the  results  are  compared  to  the  input  composite  cumu¬ 
lative  exceedances. 

The  computer  program  of  Reference  26  is  used  to  generate  the 
random  load  history  tape.  The  procedures  described  in  the  previous 
sections  are  used. 


4.3.1  Description  of  Input  Data 


A  single  mission,  defined  by  the  mission  profile  of  Figure 
12,  was  repeatedly  flown  1334  times  to  simulate  one  lifetime 
(4000  hours).  The  input  data  is  described  as  follows: 

1.  6  mission  segments 


2.  Nq/Nj-,  =  0.85  (assumed) 


3.  N0  =  5  Hz 


4.  Mean  wing  root  bending  moment  -  2.1  x  10^  in.  lb  (lg, 
limit) 


Cumulative  Exceedances /Miss ion /Segment 


F-lll  TIP  'B'  Wing 
Cumulative  Exceedances 
vs 

( A  B.M. )2 

100.0  -  Segment  1  (Climb  -  Cruise) 

Note:  (/l B.M.) ^  w/R  to  lg  Baseline 

B.M.  at  lg  =  2.1  x  10^  In. -Lb  at  Wing  Root 

K  Limit  B.M.  =  20.8  x  10^  In. -Lb  at  Wing  Root 

Limits  B.M.  =  (20.8  -  2.1)  x  10&  =  18.7  x  10&  In. -Lb 
(dB.M.)i  -  (B.M.)i  -  2.1 


Clipping  Level  « 
Limit  (4  B.M.)2 


(72  =  3.18 
Clip  ='5.88 


<7i  =  4.93 
Clip  =  3.80 


(  4B.M.  x  10-6)2  -  In. 2 -Lbs 2 

NOTE:  Upper  truncation  is  defined  by  the  clipping  level. 
No  central  truncation  was  used  other  than  that 
implied  within  the  input  data.  Test  machine  ac¬ 
curacy  provides  central  truncation  at  i/5"  of  the 
peak  load. 

Figure  13  Exceedance  Curve  Mission  Segment  1  (Climb-Cruise) 
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Figure  14  Exceedance  Curve  Mission  Segment  2  (TFR) 
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F-lll  TIP  'B *  Wing 
Cumulative  Exceedances 


vs 

( A  B.M. )2 


Segment  3  (A/G) 


Note:  (A  B.M.)  W/R  to  lg  Baseline 


Clipping  Level 
Limit  (A  B.M.) 


cr  2  =  3.62 

Clip  =  5.16 


< 7\  =  5.41 
Clip  =3.46 


(A  B.M.  x  10'°)“  ~  In. ''-Lbs' 


Figure  15  Exceedance  Curve  Mission  Segment  3  (Air-t<  -Ground) 
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Cumulative  Exceedances /Mission/Segment 


F-lll  TIP  'B '  Wing 
Cumulative  Exceedances 
vs 

(  AB.M,)2 
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Figure  16  Exceedance  Curve  Mission  Segment  4  (TFR) 
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igure  17  Exceedance  Curve  Mission  Segment  5  (Air-to-Air) 
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F-lll  TIP  'B'  Wing 
Cumulative  Exceedances 
vs 

( A B.M.)2 

Segment  6  (Cruise  -  T  &  G) 

Note:  (4B.M.)2  t0  Baseline 
T  &  G  =  Touch  and  Go 


(4  B.M.  x  1(T6)2  ~  In.  2-Lbs2 


Figure  18  Exceedance  Curve  Mission  Segment  6  (Cruise-T&G) 


Table  XI  SUMMARY  OF  INPUT  EXCEEDANCE  DATA  FOR  RANDOM  LOAD  HISTORY  SIMULATION 


5.  Truncation  levels, 

upper  =  17.0  x  10^  in4 lb 
lower  =  1.7  x  10^  in. lb 

6.  No  GAG  cycles 

7.  Maximum  clipping  value  =  +  6.0 

8.  RMS  correction  factor  =  1.1347 

The  exceedance  data  for  characterizing  the  exceedance  curve 
of  Figures  13  through  18  are  summarized  in  Table  XI.  Shown  in 
this  table  are: 

1.  Mission  segment  numbers 

2.  Number  of  positive  and  negative  RMS  levels  and  the  order 
in  which  they  are  used  in  the  random  load  simulation 

3.  Coordinates  for  the  RMS  levels  depicted  in  the  exceedance 
curves  for  each  mission  segment  type 

4.  kMS,  clipping  ratio,  and  N0  values  for  each  RMS  level 
of  each  mission  segment 

5.  Distribution  of  time  spent  in  each  RMS  level. 

Applications  of  the  input  data  and  various  aspects  of  the  genera¬ 
tion  procedure  are  discussed  in  the  following  section. 

4.3.2  Application  of  Input  Data 


The  RMS  values  shown  in  Table  XI  were  determined  using  the 
method  of  Press  (Appendix  I).  The  random  load  history  was  gene¬ 
rated,  load-by-load,  using  the  RMS  values  in  the  order  listed  in 
Table  XI. 

In  mission  segments  1,  3,  5,  and  6, all  loads  were  generated 
for  a  given  RMS  value  (positive  and  negative  spectra)  before  going 
to  a  different  RMS  value.  The  total  time  spent  in  these  segments, 
on  a  N0  basis,  was  small  compared  to  mission  segments  2  and  5 
(TFR  and  A/A  respectively).  For  mission  segments  2  and  5,  the 
RMS  values  were  not  completely  decomposed  at  the  same  time;  in¬ 
stead,  the  time  spent  in  each  RMS  value  was  distributed  so  that 
RMS  values  could  be  partially  decomposed  rather  than  all  at  once. 
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The  total  time  spent  in  a  given  RMS  level  is  still  preserved  for 
each  mission  segment  (Table  XI). 

An  N0/Np  value  of  0.85  was  assumed  because  of  the  lack  of 
wave  form  data.  This  value  is  based  on  F-lll  load  experience. 

No  GAG  cycles  were  used  in  the  random  load  simulation;  how¬ 
ever,  the  computer  program  of  Reference  26  has  the  capability  to 
include  such  cycles  in  the  load  simulation. 

A  correction  factor  of  1.1347  was  applied  to  the  input  RMS 
values  so  that  the  output  exceedances  and  RMS  values  would  match 
the  input  data.  The  correction  factor  was  determined  by  fitting 
a  fifth  degree  polynomial  through  the  results,  for  a  simplified 
case,  based  on  N0/Np  =  0.0,  0.2,  0.4,  0.6,  0.8,  and  ’.O. 

In  generating  the  test  signal,  the  random  numbers  were  alter¬ 
nately  multiplied  by  positive  and  negative  numbers.  For  the 
Gaussian  case  (N0/Np  =  0.0)  there  is  an  equal  probability  of 
selecting  a  positive  or  negative  delta  load.  The  sign  alteration 
is  deterministic  but  the  random  numbers  have  an  equal  probability 
of  being  positive  or  negative.  Without  some  compensation,  the 
procedure  generates  an  actual  N0/Np  =  0.5  for  the  Gaussian  case. 
Inerefore,  a  load  signal  adjustment  is  required  to  preserve  the 
input  exceedances.  No  adjustment  is  required  when  N0/Np  equals 
one  (Rayleigh  case). 

The  correction  factor,  in  effect,  requires  that  a  greater 
number  of  loads  be  created  to  simulate  the  load  history.  The  use 
of  such  a  factor  was  justified  in  this  case  since  the  output  ex¬ 
ceedances  agreed  with  the  input  exceedances  and  the  essential 
wave  form,  N0/Np,  was  preserved. 


4.3.3  Statistical  Analysis  of  Load  Simulation 


A  statistical  analysis  was  performed  to  justify  the  random 
load-history  simulation  for  the  actual  load  history.  A  composite 
cumulative  exceedance  curve  (Figure  19)  was  plotted  for  the  posi¬ 
tive  spectra  using  the  data  of  Table  XII.  The  composite  cumu¬ 
lative  exceedances  were  determined  by  summing  the  ordinates  to 
the  cumulative  exceedance  curves  of  Figures  13  through  18  for 
given  (  AB.M.)^  values.  The  resulting  cumulative  exceedances  are 
summarized  in  Table  XIII. 

The  BY4  computer  program  generates  random  loads  one  at  a 
time  and  stores  them.  Subroutines  are  employed  which  count 
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Table  XII  DATA  FOR  ACTUAL  COMPOSITE  EXCEEDANCE  CURVE 


1*1 
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B.M.  =  Resultant  wing  bending  at  SS  70.3 

*  =  Estimated  value  beyond  truncation  level 
**  =  TFR  loads  are  truncated  as  shown  due  to 
an  automatic  fly-up  command  that  is  g 
level  limited. 


Positive  Composite  Cumulative  Exceedances /Mission 
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- , 

< 

100.0 

[ 

Notes:  1.  Based  on  F-lll  TIP  * B *  Spectrum 

2.  The  Actual  Composite  Cumulative  Exceedance  Curve 
Shown  is  Based  on  the  Summation  of  the  Positive 
Exceedances  for  Mission  Segments  1  through  6 
(reference  Figure  13  through  18  and  Table  13) 
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Analysis  of  the  Random  Load  Simulation 
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Figure  19  Comparison  of  Load  Simulation  Composite  Cumulative 
Exceedances  Versus  Actual  Exceedances 


Load 


Figure  20  Random  Load  Statistics 
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Table  XIII  Summary  of  Composite  Cumulative  Exceedances 
Based  on  Statistical  Analysis  of  Simulated 
Random  Lead  Test  Tape 
(Continued) 

Notes:  a.  Obtained  using  computer  program  of  Reference  26 

b.  Statistics  based  on  analysis  of  119,991  load  records 
from  the  test  tape. 

•  Total  p’tmber  of  load  records  on  the  test  tape  = 

.  873,465  (1  lifetime) 

•  Total  number  of  positive  load  records  on  test  tape 
436,732  (1  lifetime) 

•  Equivalent  number  of  positive  load  records  for 
Rayleigh  case  (N0/Np)  =  383,258 


1.  The  number  of  loads  in  the  simulation 

2.  The  positive  slope  crossings  of  the  mean  and 
at  each  load  level  crossed 

3.  The  number  of  load  exceedances  at  given  load 
levels  above  and  below  the  mean  (Figure  20). 

The  load  statistics  are  analyzed  in  batches  of  40,000  loads. 
All  loads  on  the  tape  can  be  analyzed  by  pooling  the  data  for  each 
load  batch.  A  sample  output  of  the  load-time  history  simulation 
is  shown  in  Figure  21.  The  exceedances,  based  on  the  load  simula¬ 
tion  for  1  mission  and  300  missions,  are  also  compared  against  the 
actual  exceedances  in  this  figure. 

The  composite  cumulative  exceedances,  based  on  the  statis¬ 
tical  analysis  of  119,991  loads  on  the  test  tape,  are  summarized 
in  Table  XIII.  These  results  were  generated  using  the  BY4  com¬ 
puter  program  of  Reference  26.  To  compare  the  simulated  load 
history  against  the  actual  load  history,  the  composite  cumulative 
exceedances  of  Table  XII  were  plotted  against  the  actual  cumula¬ 
tive  exceedances  (Figure  19).  The  agreement  between  the  simulated 
load  history  and  the  actual  load  history,  as  shown  in  Figure  21, 
was  excellent.  A  value  of  N0/Np  equal  to  0.877  was  obtained  for 
the  simulated  load  history.  This  compares  with  an  input  value  of 
0.85. 

It  was  concluded  that  the  random  load  history  simulation  used 
was  a  valid  representation  of  the  actual  load  history. 

4.3.4  Test  Tape 


The  nndom  load  history  simulation  (digital)  is  stored  on  a 
magnetic  test  tape.  Successive  loads  on  the  tape  are  connected 
by  a  havers ine  wave  form.  The  shape  of  the  load  signal  is  in¬ 
variant  and  the  amplitude  is  set  by  calibrating  the  maximum  speci¬ 
men  load  to  the  maximum  signal  peak.  With  this  set-up, the  same 
test  tape  can  be  used  for  different  size  specimens  and  components. 
The  specimen  load  calibration  procedure  is  discussed  in  subsection 
4.4. 

The  test  tape  simulates  one  lifetime  of  random  load  history 
with  873,465  load  records.  The  numbers  of  loads  used  in  a  test 
are  counted  and  the  lifetime  is  computed  by  dividing  the  number 
of  loads  used  by  873,465  loads /I i fetime .  When  all  the  loads  have 
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Figure  21  Composite  Exceedance  Correlation  and  Sample  Random 
Load/Time  Signal  from  Test  Tape 


been  read  on  the  test  tape,  the  tape  is  rewound  and  the  process 
is  continued  until  the  desired  lifetime  is  reached  or  the  speci¬ 
men  fails. 


4.4  SPECIMEN  LOAD  CALIBRATION 


The  loading  signal  must  be  calibrated  with  respect  to  the 
maximum  bending  moment  on  the  tape  and  the  desired  maximum  load 
on  the  specimen.  Since  the  signal  on  the  tape  is  scaled  to  the 
maximum  peak,  the  signal  can  be  calibrated  for  axial  load,  shear, 
torsion,  etc.  (separately). 

The  design  limit  (D.L.)  load  for  the  1/5-  and  1/2-scale  spec 
mens  is  66. 7K  and  157. OK,  respectively.  The  loading  signal  is 
calibrated  to  these  loads.  These  maximum  loads  are  called  ”100 
percent  calibration  loads.” 

The  1/5-scale  specimen  was  designed  for  a  100K  ultimate  load 
The  maximum  limit  load  on  this  specimen  is  thus  66. 7K.  The  1/2- 
scale  specimen  100  percent  calibration  load  was  determined  using 
the  1/5-scale  limit  load  and  ratioing  the  1/5-  and  1/2-scale  spec 
mens  cross-sectional  areas  as  shown  below. 

Pj)L  (1/5  scale)  =  66. 7K  limit 

1/5-scale  cross  section:  5”  wide  by  40  plies  thick 

1/2-scale  cross  section:  5”  wide  by  94  plies  thick 

Both  specimens  have  the  same  width;  therefore,  the  calibration 
load  for  the  1/2-scale  specimen  can  be  computed  by  ratioing  the 
number  of  plies  for  the  1/5-  and  1/2-scale  specimens. 

94 

PDL  (1/2  scale)  =  —  x  66.7  =  157. OK 

With  the  above  calibration  loads>  each  specimen  is  subjected  to 
compatible  stress  histories.  This  puts  the  1/5-  and  1/2-scale 
specimens  on  the  same  baseline  and  allows  a  direct  comparison  of 
the  fatigue  results.  This  is  important  for  evaluating  the  scal¬ 
ing  effect. 

The  closed  loop  test  system  can  be  used  to  test  different 
specimen  sizes  s imultaneous ly  using  the  same  load  type.  The 
load  signal  is  calibrated  for  each  specimen  using  individual 
servo  amplifiers. 
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A  load  truncation  level  was  selected  on  the  following  basis: 


1.  Adequate  excursions  must  be  made  to  ensure  a  successful 
time  trace  for  failures  in  time. 

2.  The  truncation  level  does  not  exclude  the  series  of 
damage  producing  high  loads. 

The  truncation  value  would  be  exceeded  once  per  five  missions. 
This  is  equivalent  to  approximately  0.211  exceedances  per  mission 
(reference  Table  XII  and  Figure  19)  or  280  exceedances  per  life¬ 
time,  This  level  ensures  sufficient  excursions  to  the  peak  load 
to  develop  an  adequate  lifetime  distribution. 

For  this  program,  upper  and  lower  truncation  values  of  17.0 
x  106  in#  and  1.7  x  10^  in#,  respectively,  were  assumed.  These 
levels  satisfy  the  criteria  above.  The  limit  wing  root  bending 
moment  is  20.8  x  10^  in#.  Using  these  values  the  maximum  and 
minimum  load  on  the  1/5-  and  1/2-scale  specimens  was  computed  as 
shown  below. 

1/5  scale 

PTRU  =  x  66. 7K  =  54. 5K 

PTRL  =  2O  x  66*7K  =  5-45K 

1/2  scale 

PTRU  =  |o  X  157 -0K  =  128.0K 

pTRL  =  2O  x  157. OK  =  12. 8K 


SECTION 


V 


TEST  PROGRAM 


EXECUTION 


5 . i  INTRODUCTION 


A  description  of  the  tests  conducted  and  the  results  obtained 
are  given  in  th-'s  section.  The  types  of  tests  conducted  are  sum¬ 
marized  in  subsection  2.4. 


5.2  CONSTANT  AMPLITUDE  TESTS 


In  this  subsection,  the  constant  amplitude  tests  conducted 
are  described  and  the  results  of  these  tests  are  summarized. 


5.2.1  Tension-Tension 


The  purpose  of  the  tens  ion- tens  ion  tests  ras  to  generate 
constant  amplitude  load  versus  cycles  to  failure  data  needed  to 
predict  the  specimen  fatigue  life  for  a  given  loading  spectrum 
based  on  Miner's  r«*2e.  Twelve  tests  were  performed  using  1/5- 
scale  specimen  (R  =  G.10).  Three  sets  of  constant  amplitude  load¬ 
ing  were  used:  (1)  3  to  JOK.  (2)  4  to  40K,  and  (3)  5  to  50K. 

Only  one  specimen  (K012432)  was  tested  using  the  4  to  40K  loading 
because  of  the  failure  of  one  of  the  titanium  lugs  under  this 
loading.  The  lug  cracked  at  one  end  of  specimen  K994567  under  a 
3  to  30K  cyclic  loading.  The  cracked  lug  problem  is  discussed  in 
paragraph  5.2.3. 

Each  tension-tension  specimen  was  mounted  in  fixture  5  us 
shown  in  Figure  22.  Only  specimens  K012465  and  K99456G  had  strain 
gauges,  and  these  were  read  only  for  the  strain  surveys.  The  ap¬ 
plicable  cyclic  loading  was  continuously  applied  to  the  specimen 
at  3  Hz  until  failure  c'jurred,  The  specimens  tested,  the  cyclic 
loading  used,  the  cycles  to  failure,  and  the  test  completion  dates 
are  summarized  in  Table  XIV,  Phonographs  of  typical  failed  speci¬ 
mens  are  included  in  Appendix  IV. 


57 


t  Amplitude  Figure  23  Tension-Compression  Constant 

train  Survey  Amplitude  Setup  with  Center 

Support  Shown  Removed 
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Table  XIV  SUMMARY  OF  1/5- SCALE  SPECIMEN 
CONSTANT  AMPLITUDE  LOADING  TEST  RESULTS 


Specimen 

I.D. 

Type  Test 

Cyclic  Loading 
(KIPS) 

Cycles 

To 

Failure 

K994587f£ 

Tension- 

Tension 

3  to 

30 

202,667 

3-22-71 

K012444f 

113,938 

3-24-71 

K01246lf 

l,296,000a 

4-15-71 

K994567f 

3  to 

30 

598, 310^ 

5-3-71 

K012482 

4  to 

40 

101,750° 

5-4-71 

K012465f 

5  to 

50 

7,545 

3-17-71 

K9945664f 

5,458 

3-17-71 

K012476S 

5,210 

5-3-71 

K012477S 

31,800 

5-12-71 

K012478S 

11,071 

5-12-71 

K012483S 

2,205 

5-12-71 

K900445$g 

Tension- 

Tension 

5  t< 

3 

50 

25,643 

5-14-71 

K012445 

Tension- 

Compression6 

-6  to 

30 

494, 580^ 

4-22-71 

K012442fg 

-8  to 

40 

26,674 

wasBBm 

K994568$g 

-10 

to  50 

2,665 

3-15-71 

K012475i»g 

1,790 

5-6-71 

K0124384 

1,839 

5-6-71 

K012462fjg 

3,193 

5-10-71 

K900442  g 

Tension-Compressione 

-10 

to  50 

2,959 

5-11-71 

Notes:  a 
b 
c 
d 
e 


No  specimen  failure 
Lug  cracked  out  at  one  end 
Complete  lug  failure  at  one  end 
Also  used  for  strain  survey 

Specimen  supported  at  the  center  by  two  channels  with 
teflon  rub  strips  next  to  specimen 


f  Built  prior  to  manufacturing  procedure  change  for 
leveling  the  titanium  lug  ends  in  the  tool 
(reference  paragraph  3.3.2), 
g  Photographs  of  failed  specimens  given  in 
Appendix  IV. 
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5.2.2  Tension-Compression 


Tension-compression  constant  amplitude  tests  were  performed 
to  (1)  characterize  the  fatigue  life  of  the  1/5-scale  specimens 
under  load  reversals  and  (2)  compare  the  fatigue  damaging  effects 
of  tension-tension  versus  tension-compression  cyclic  loading. 

Press-fit  bushings  were  used  in  the  lug  ends  to  minimize  the 
"play"  between  the  lugs  and  pins.  This  precaution  was  taken  to 
protect  the  lugs  against  the  dead  space  induced  impact  caused  by 
load  reversals. 

Specimens  were  mounted  in  fixture  6  using  the  center  supporc 
option.  The  center  support  was  needed  to  stabilize  the  specimen 
under  compressive  load.  The  test  set-up  is  the  same  as  that  shown 
in  Figure  23  except  specimen  K994568  was  the  only  one  with  strain 
gauges  and  these  gauges  were  read  only  for  the  strain  survey. 

The  tests  were  conducted  using  three  sets  of  constant  ampli¬ 
tude  loading:  (1)  -6  to  30K,  (2)  -8  to  40K,  and  (3)  -10  to  50K. 
The  lug  cracked  at  one  end  of  specimen  K012445.  This  problem  is 
discussed  in  paragraph  5.2.3. 

The  test  results  are  shown  in  Table  XIV.  Typical  specimen 
failures  are  shown  in  Appendix  II. 


5.2.3  Lug  Problem 


Lug  cracks  developed  in  specimens  K012445,  K994567,  and 
KOI  2482.  Two  of  the  specimens  are  shown  in  Figures  24  and  25. 

A  fatigue  analysis  was  performed  with  and  without  the  effects 
of  the  press  fit  bushings  (R  =  -0.20)  using  Miner's  rule.  No  lug 
failure  could  be  justified  by  this  analysis. 

Similar  lug  failures  have  been  observed  at  low  stress  levels 
in  D6ac  specimens  (Reference  31).  In  this  case,  lug  failures  oc¬ 
curred  at  low  stress  levels  even  though  the  lug  net  section  was 
more  than  ten  times  larger  than  the  minimum  cross  section.  This 
phenomenon  was  unexplainable,  but  it  could  be  eliminated  by  in¬ 
creasing  the  lug  net  section  stress. 

To  protect  the  titanium  lugs,  the  cyclic  load  levels  were  in¬ 
creased  from  3  to  30K  to  5  to  50K  and  from  -6  to  30K  to  -10  to 
50K.  Changing  the  cyclic  loading  did  not  affect  the  purpose  of 


60 


the  constant  amplitude  tests.  No  lug  failures  were  experienced 
after  these  load  changes.  Cracked  or  failed  lug  ends  were  not 
repaired  and  reused  in  future  specimens. 


5.3  STATIC  TESTS 


Eleven  1/5-scale  specimens  were  static  tested  using  the 
Baldwin-Tate-Emery  Universal  test  rachine.  Seven  of  the  eleven 
specimens  tested  were  previously  rejected  due  to  excessive  war- 
page.  Specimen  K012440  was  partially  straightened  by  heating 
under  load  before  the  static  test.  The  test  set-up  is  shown  in 
Figure  26  except  only  three  specimens  were  instrumented. 

The  static  test  results  are  shown  in  Table  XV  for  all  speci¬ 
mens  tested.  The  average  failing  load  is  given  with  and  without 
the  results  for  the  rejected  specimens.  Pictures  of  typical 
failed  specim'  s  are  shown  in  Figure  27. 


5.4  RANDOM  LOAD  TESTS 


Random  load  tests  were  performed  to  characterize  the  fatigue 
life  of  1/5-  an  i  1/2-scale  specimens.  The  random  load-history 
simulation  (stored  on  magnetic  tape)  for  the  spectrum  described 
in  Section  IV  was  used  with  the  closed  loop  test  system  shown  in 
Figures  3  and  4, 

Residual  strength  tests  and  fatigue- to- failure  tests  were  per¬ 
formed.  The  residual  strength  tests  were  conducted  by  first  sub¬ 
jecting  the  specimen  to  a  prescribed  random  load  history  simulating 
a  fractional  service  lifetime  and  then  determining  the  static 
strength  of  the  specimen. 

The  desired  random  load  history  was  put  on  the  specimens  using 
three  test  fixtures  for  the  1/5-scale  specimens  and  one  test  fix¬ 
ture  for  the  1/2-scale  specimens.  The  test  system  used  permitted 
the  simultaneous  testing  of  three  1/5-scale  specimens  and  one  1/2- 
scale  specimen.  When  all  fixtures  were  used  simultaneously,  a  test 
rate  between  3  and  4  Hz  was  possible.  A  5  Hz  test  rate  could  be 
achieved  with  only  one  or  two  specimens.  The  number  of  fixtures 
used  at  the  same  time  depended  upon  the  specimen  supply.  Due  to 
the  specimen  recycling  time,  only  two  or  three  specimens  were 
usually  tested  at  the  same  time.  However,  all  fixtures  were  used 
dui ing  some  periods. 
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Table  XV  SUMMARY  OF  1/5  SCALE  SPECIMEN 
STATIC  TEST  RESULTS 


SPECIMEN 

I.D. 

Date 

Tested 

Failing 

Load  (KIPS) 

Average 

Failing  Load(KIP$ 

K994588  d ,  e ,  f 

3-17-71 

86.6 

1 

1 

■ 

K012439  d  ,e ,  f 

3-17-71 

109.8 

I 

1 

1 

K012441  d ,e , f 

3-17-71 

106.4 

1 

1 

1 

1  100.9 

K012440a 

3-30-71 

119.0 

■ 

■ 

■ 

K994578b 

4-1^ 

t-71 

89.2 

K012459b 

110.4 

KOI 2446 b 

111.2 

K012463b 

98.0 

K012437b 

108.5 

K012474b 

99.0 

K994569b 

4-1^ 

i-71 

92.0 

102.7 

Notes:  a  Specimen  was  straightened  to  reduce  dimensional  vari¬ 
ations  prior  to  testing 

b  Rejected  specimen  with  excessive  dimensional  variations 
c.  Specimen  are  listed  i.  the  order  tested 


d  Specimen  also  used  for  static  strain  survey 

e  Built  prior  to  manufacturing  procedure  change  for 
leveling  the  titanium  lug  ends  in  the  tool 
(ref.  paragraph  3.3.2). 

f  Photographs  of  failed  specimens  given  in  Appendix  IV. 
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Details  of  the  1/5-scale  and  1/2-scale  specimen  tests  and 
the  results  are  discussed  in  paragraphs  5.4.1  and  5.4.2  ,  respec¬ 
tively. 


5.4.1  One-Fifth  Scale  Specimens 


5.4. 1.1  Residual  Strength  Tests 

Residual  strength  tests  were  performed  using  specimens  that 
were  first  subjected  to  a  1,  a  10,  and  a  50  percent  lifetime  of 
random  loading.  The  percentage  lifetime  load  history  is  based  on 
the  use  of  a  fractional  amount  of  the  total  load  records  stored 
on  the  magnetic  tape. 

Thirty-one  residual  strength  tests  were  performed.  The  speci¬ 
men  breakdown  is  as  follows: 

1.  5  at  1  percent  lifetime 

2.  6  at  10  percent  lifetime 

3.  20  at  50  percent  lifetime. 

The  test  set-up  for  the  random  loading  is  shown  in  Figure  22 
(no  strain  gauges).  A  Baldwin-Tate-Emery  Universal  test  machine 
was  used  to  evaluate  static  strength.  Table  XVI  summarizes  the 
specimens  used,  the  types  of  testing,  the  residual  strengths,  and 
the  test  completion  dates.  Photographs  of  typical  specimen 
failures  are  included  in  Appendix  IV. 

5 . 4 . 1 . 2  Fatigue-to-Failure  Tests 

Twenty  fatigue- to-failure  tests  were  conducted  using  the 
random  load-history  simulation  described  in  Section  IV.  These 
tests  were  performed  to  determine  the  lifetimes  of  1/5-scale 
boron-to-titanium  double  scarf  adhesive  bonded  joints  under  random 
loading.  This  information  was  needed  for  a  statistical  data  base. 

For  these  tests,  a  specimen  was  placed  in  a  fixture,  as  shown 
ir.  Figure  22,  and  subjected  to  a  continuous  random  loading  until 
failure.  The  specimens  used,  the  lifetimes  to  failure,  and  the 
test  completion  dates  are  shown  in  Table  XVII.  Photographs  of 
typical  specimen  failure  are  presented  in  Appendix  IV. 
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TABLE  XVI  SUMMARY  OF  1/5  -SCALE  SPECIMEN 
RESIDUAL  STRENGTH  TESTS  AFTER  RANDOM  LOADING 


Specimen11 

I.D. 


Type  Test 


Residual 

Strength 

(KIPS) 


Test 

Complete 

Date 


K900448  g 
K90044lg 
K904958g 
K905016g 
K904955g 


1%  lifetime3 


17o  lifetime3 


104.4 
114.0 

108.4 
109.1 

90.2 


5-20-71 

5-20-71 

5-26-71 

5-26-71 

5-26-71 


K900447 g 
K904950g 
K905017  g 
K904957g 
K904959g 
K905021 


107,  lifetime^ 

* 


*  b 

10%  lifetimeb 


79.4 

105.0 

110.6 

104.6 
106.0 

108.6 


5-21-71 

5-21-71 

5-26-71 

5-26-71 

5-26-71 

11-9-71 


K900449g 

K900450g 

K904951g 

K904952g 

K904953g 

K904954g 

K904956g 

K905015g 

K905022g 

K905023g 

K905018g 

K905020 

K905362 

K905363 

K905365 

K905366 

K905369 

K905371 

K905467 

K905469 


507o  lifetime2 


507,  lifetime' 


115.6 

73.5 

79.6 

86.0 

106.0 

96.0 

102.0 

95.0 

110.0 

111.9 

116.7 

110.4 
111.0 
105.0 

99.3 

113.0 

103.5 
92.5 


5-24-71 

5-24-71 

5-26-71 

5-26-71 

5- 26-71 

6- 1-71 


,  * 

6-1-71 

11-9-71 

k 


t 

11-9-71 

11-17-71 

11-17-71 

11-17-71 


Notes:  a  Equal  to  8735  records  from  random  load  tape 
b  Equal  to  87347  records  from  random  load  tape 
c  Equal  to  436,733  records  from  random  load  tape 
d  Specimen  Tailed  at  308,230  records  (35.37,  lifetime) 
e  Specimen  failed  at  347,577  records  (39.87c  lifetime) 
f  Specimen  are  listed  in  the  order  tested  for  each  type 
test . 

g  Photographs  of  failed  specimens  given  in 
Appendix  IV. 
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Table  XVII  SUMMARY  OF  1/5 -SCALE  SPECIMEN  RANDOM 
LOADING  FAT IGUE -TO-FAILURE  TEST  RESULTS 


Specimen 

I.D. 

Records 

to 

Failure 

Lifetime3 

Test 

Complete 

Data 

K012464c,d 

1,067,867 

1.22 

4-23-71 

K012443  c,d 

883,438 

1.01 

4-27-71 

K994565 

1,947,900 

2.23 

4-28-71 

K0 12480 d 

296,434 

0.34 

5-7-71 

K012479d 

634,950 

0.73 

5-10-71 

KQ12484d 

401,362 

0.46 

5-12-71 

K900443 d 

718.310 

0.82 

5-12-71 

K012481d 

633,168 

0.73 

5-14-71 

K012485d 

565,738 

0.65 

5-20-71 

K900446 

768,087 

0.88 

5-21-71 

K905368 

3,549,533 

4.06 

11-19-71 

K905367 

3,280,383 

3.76 

11-22-71 

K905470 

708,000. 

0.81 

11-24-71 

K905364 

4,446 , 134b 

5.07 

12-1-71 

F504622 

1,227,120 

1.40 

12-7-71 

F504418 

2,112,026 

2.41 

12-16-71 

K905468 

3,203,241 

3.66 

12-19-71 

F504421 

2,885,579 

3.30 

12-21-71 

F504419 

1,256,040 

1.43 

1-4-72 

F504417 

1,108,994 

1.27 

1-6-72 

Notes:  a  One  lifetime  =  873,465  records  (one  time  through  the 
random  load  tape) 

b  Random  load  test  was  stopped  and  specimen  was  static 
tested  to  109. 5K 


c  Built  prior  to  manufacturing  procedure  change  for 
leveling  the  titanium  lug  ends  in  tool  (ref.  para¬ 
graph  3.3.2). 

d  Photographs  of  failed  specimens  given  in  Appendix  IV. 


67 


5.4.2  One-Half  Scale  Specimens 

Five  10-percent  lifetime  residual  strength  tests  and  five 
fatigue-to-failure  tests  were  performed.  The  tescs  employed  the 
same  random  load  history  (with  a  different  load  calibration  than 
the  1/5-scale  specimens  to  maintain  the  same  stress  level  in  both 
specimens) and  procedure  used  for  the  1/5-scale  specimens.  The 
test  set-up  for  putting  the  random  loading  on  the  specimens  and 
performing  the  fatigue-to-failure  tests  is  shown  in  Figure  28. 
These  tests  were  conducted  to  make  a  preliminary  evaluation  of 
the  scaling  effect  between  the  1/5-  and  1/2-scale  specimens. 

The  residual  strength  and  the  fatigue-to-failure  test  results 
are  given  in  Tables  XVIII  and  XIX,  respectively.  Typical  speci¬ 
men  failures  are  shown  in  Appendix  V. 


1 


{ 


Table  XVIII.  SUMMARY  OF  1/2 -SCALE  SPECIMEN 
RESIDUAL  STRENGTH  TESTS  AFTER  RANDOM  LOADING  (10%  LIFETIME)3 


Specimen 

I.D. 

Residual 

Strength 

(KIPS) 

Test 

Complete 

Date 

K012453 

236.5 

12-3-71 

K012454 

240.0 

12-3-71 

F504624b 

168.5 

2-8-72 

F504625b 

206.0 

2-8-72 

F504618b 

219.0 

2-8-72 

Notes:  a  Equal  to  87347  records  from  random 

load  tape.  Test  fixture  number  12 
used  for  random  loading. 

b  Photographs  of  failed  specimens  given 
in  Appendix  V. 
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Table  XIX  SUMMARY  OF  1/2 -SCALE  SPECIMEN  RANDOM 
LOADING  FAT IGUE-TO-FAI LURE  TEST  RESULTS 


SPECIMEN 

I.D. 

RECORDS 

TO  FAILURE 

LIFETIMES 

TEST  COMPLETE 
DATE 

F504620 

450,258 

0.57 

12-10-71 

F504622  b 

685,366 

0.78 

12-23-71 

F504619 b 

387,884 

0.44 

12-30-71 

F504621 b 

198,220 

0.23 

1-4-72 

F504623  b 

730,151 

0.84 

1-13-72 

Notes:  a  One  lifetime  *  873,465  records  (one  time 

through  the  random  load  tape) .  Test  fixture 
number  12  used  for  random  loading. 

b  Photographs  of  failed  specimens  given  in 
Appendix  V. 
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SECTION  VI 


EVALUATION  OF  TEST  RESULTS 


The  test  results  compiled  from  the  tests  described  in  Section 
V  are  analyzed  in  this  section. 

An  assessment  of  Miner's  rule  for  computation  of  life  under 
random  loadings  was  made.  The  major  portion  of  the  test  program 
was  devoted  to  evaluation  of  the  lifetime/residual  strength  charac¬ 
teristics  of  the  1/5-  and  1/2-scale  bonded  joints.  These  evalua¬ 
tions  included  residual  strength  and  lifetime  distribution  func¬ 
tions  and  an  evaluation  of  scale  effects. 


6.1  FATIGUE  ANALYSIS  FOR  ONE-FIFTH  SCALE  SPECIMEN 


A  fatigue  analysis,  based  on  Miner's  linear  damage  rule,  is 
presented  in  this  section  for  the  1/5-scale  specimen.  This  analy¬ 
sis  employs  a  load  versus  cycles  to  failure  curve  based  on  static 
and  constant  amplitude  fatigue  tests  and  the  same  loading  spectrum 
used  for  the  random  load  test  tape.  The  predicted  fatigue  life 
is  compared  against  the  mean  lifetime  obtained  for  the  1/5-scale 
fatigue-to-failure  tests  (Reference  Table  XVII).  This  comparison 
is  used  to  evaluate  the  effectiveness  of  Miner's  rule  for  predict¬ 
ing  the  fatigue  life  of  composite  bonded  joints. 


6.1.1  Procedure  Used 


The  fatigue  analysis  is  based  on  the  following  procedure: 

1.  Divide  up  the  cumulative  exceedance  curve  versus 
A  B.M.  into  six  segments  as  shown  in  Figure  29. 
Determine  the  total  bending  moment  for  each  segment 
using, 

(B.M.^  =  (  AB.M.)i  +  2.1 

where  (B.M.).  =  total  bending  moment  (xlO  ^  in .#) 

at  the  wing  root  for  segment  i. 

(AB.M.).  =  delta  bending  moment  measured  from 

1  the  vertical  axis  of  Figure  10  to  the 
center  of  segment  i. 
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Composite  Cumulative  Exceedance/Mission 


100 


£ 

w 

* 


o 

c 

o 

a 

•H 

o 

o 

a 

w 


1  x  10^ 


y  Goes  Co  1 
-y-~  for  11  St 
^  at  0.250 

02. 7K  (Average 

1/5  Sc 

ale  Boron  -  Tita 

0 

nium  Specimens  _ 

atlC  IwSLb  y 

Cycles 

R  -  .1 

/-Single  T 
/  Cyclic  L 

!St  (A  K  to  A0  K 
oading) 

Average 
(5K  to 

...  _ 

for  6  Tests  / 
OK  Cyclic  Losdin 

L  _  _ 

i) 

\  1 

HH 

HH 

HH 

\ 

Average  for 
(3 K  to  30K  Cycl 

_ 1 _ L-JL-0-LLI 

4  Tests 
ic  Loading) 

_ I— i _ !— LLili 

1  x  103  1  x  104  lx  105  1  x  106 

Number  of  Cycles  Co  Failure 


1  x  10/ 


Figure  30  Load  Vs.  Cycles  to  Failure 


2.1 


~6 

bending  moment  (xlO  in.#)  at  the 
wing  root  of  Ig. 

The  total  maximum  limit  bending  moment  at  the  wing  root 
is  20.8  x  10^  in.  lbs.  An  upper  truncation  bending 
moment  of  17.0  x  10^  in.  lbs  was  used  for  the  test  tape; 
therefore,  the  last  segment  in  Figure  29  is  at  14.9  x 
lO^  in.  lbs  (14.9  +  2.1  =  17.0).  The  lower  truncation 
bending  monent  for  the  test  tape  was  1.7  x  10^  in.  lbs. 

2.  Using  Figure  29,  determine  the  total  number  of  times  per 
mission  that  the  bending  moment  for  each  segment  occurs. 
The  total  number  of  occurrences  is  equal  to  the  differ¬ 
ence  between  the  ordinates  of  the  curve  shown  in  Figure 
29  at  the  beginning  and  end  of  each  segment.  The  total 
number  of  cycles  (n^)  for  a  given  bending  moment  is 
equal  to  the  total  number  of  occurrences  for  each  seg¬ 
ment.  There  are  1334  missions  in  a  lifetime.  The  total 
number  of  occurrences  in  a  lifetime,  for  a  given  bending 
moment,  is  equal  to  1334  times  the  number  of  occurrences 
per  mission. 


3,  Determine  the  load  on  the  specimen  for  each  segment  shown 
in  Figure  29  using 


P 


load 


(B.M.)i  x  66. 7K 

2078 


where  66. 7K  is  the  limit  load  for  the  1/5-scale  specimen. 
A  P^oa(j  value  is  computed  for  each  segment  shown  in 

Figure  29. 


4.  Develop  a  load  versus  cycles  to  failure  curve  using  the 

test  results  summarized  In  Table  XIV.  This  curve  is  shown 
in  Figure  30, 


5.  Determine  the  number  of  cycles  (N^)  to  failure  using 
Figure  30,  corresponding  to  the  ?i0ad  value  for  each 
segment  shown  in  Figure  29. 

n . 

6.  The  fatigue  damage  for  each  segment  of  Figure  30,  (— ) , 
is  determined  separately  and  then  summed  to  get  i 
the  total  fatigue  damaged  .  The  predicted  fatigue 

"N£ 

life  of  the  1/5-scale  fatigue-to-failure  specimen  is 
computed  as  follows: 
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Fatigue  Life  = 


6.1.2  Fatigue  Life  Computations  Using  Miner's  Rule 


The  fatigue  analysis  performed  is  shown  in  Table  XX.  The 
results  are  evaluated  in  paragraph  6.2.3  to  determine  the  effec¬ 
tiveness  of  Miner's  rule  for  predicting  the  fatigue  life  of  bonded 
composite  joints  under  random  load. 


6.1.3  Effectiveness  of  Miner's  Rule 


A  fatigue  life  of  5.12  lives  was  predicted  for  the  1/5-scale 
fatigue-to-failure  specimen  based  on  Miner's  rule.  This  compares 
with  a  mean  lifetime  of  1.812  (Reference  Table  XVII)  based  on 
twenty  1/5-scale  fatigue-to-failure  specimens  subject  to  a  random 
fatigue  loading  (Test  Tape).  Miner's  rule  in  this  case  overpre¬ 
dicts  the  fatigue  by  a  factor  of  2.82. 

The  predicted  fatigue  life,  based  on  the  analysis  presented, 
should  be  a  conservative  prediction  (shorter  life)  in  view  of  the 
following: 

1.  The  analysis  conservatively  assumes  that  one  random 
load  occurrence  is  equivalent  to  one  load  cycle.  The 
definition  of  a  load  cycle  is  not  clear  for  a  random 
loading.  If  one  assumes  three  random  loads  constitutes 
one  load  cycle,  then  an  equivalent  number  of  load  cycles 
can  be  computed  using 


where  NE  =  equivalent  number  of  load  cycles  (Reference 
Figure  31). 

n^  =  total  number  of  loads  in  random  history 
-  1  =  number  of  1/2  load  cycles. 
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5 


obtained  from  Figure  30  for  given  Pload  value, 
erence  discussion  in  paragraph  6.2.3. 


Load 


Figure  31  Equivalent  Random  Load  Cycle 


Load 


Figure  32  Stress  Ratio  (R)  Variation  for  Random  Load  Test  Tape 


lb 
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If  the  equivalent  load  cycle  concept  is  used  in  the 
fatigue  analysis,  less  fatigue  damage  would  result; 
consequently,  a  longer  fatigue  life  would  be  predicted. 


2. 


The  stress  ratio,  R  (minimum  stress /maximum  stress  in 
a  fatigue  cycle),  for  the  equivalent  random  load  cycle 
is  larger  over  the  load  history  (more  positive)  than  R 
for  the  constant  amplitude  cycles.  The  test  tape  random 
load  history  had  a  mean  of  2.1  x  10^  in.  lbs  with  an 
upper  and  lower  truncation  of  17.0  x  10”  in.  lbs  respec¬ 
tively.  The  range  of  possible  values  for  Pmax  and  Pm^n 

is  2.1  x  106  in.  lbs  to  17.0  x  106  in.  lbs  and  1.7  x  106, 
respectively.  This  means  that  R  is  equal  to  0.10  only 
when  Pmax  and  Pm^n  are  equal  to  the  truncation  value  in 
the  load  cycle  The  other  possible  ratios  of  Pmax  to 
Pmin  will  be  greater  than  0.10  over  the  service  history 
(Figure  32). 


Fatigue  damage  decreases  as  R  increases  (becomes  more  posi¬ 
tive).  Since  the  fatigue-to-failure  specimens  for  the  random 
fatigue  loading  experienced  R  values  greater  than  0.10  during 
testing,  they  should  exhibit  a  longer  fatigue  life  than  the  case 
where  R  is  a  constant  0.10  during  testing.  The  2.82  factor  be¬ 
tween  the  predicted  fatigue  life,  based  on  Miner's  rule,  and  the 
actual  mean  life  obtained  by  random  testing  is  considered  to  be 
low  in  view  of  the  discussion  above. 


The  results  presented  in  this  section  indicate  that  Miner's 
rule  is  inadequate  for  predicting  the  fatigue  life  of  bonded 
composite  joints  under  random  loading. 


6.2  EVALUATION  OF  THE  RANDOM  FATIGUE  DATA 


Two  steps  were  executed  in  the  reduction  of  the  random  fatigue 
data.  First, the  load  and  lifetime  traces  were  plotted  on  Weibull 
coordinates  to  observe  correlation  of  the  observations  with  a 
Weibull  distribution.  Second,  a  first-order  analysis  of  the  ob¬ 
served  phenomena  was  made  on  the  basis  of  a  damage  growth  hypoth¬ 
esis.  Positive  correlation  has  been  achieved  with  the  damage 
growth  model. 

The  ranked  raw  data  is  given  in  the  following  tables: 

1.  Pocied  static  and  1  percent  lifetime  data 
1/5  scale  -  Table  XXI 
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Table  XXI  RANKING  OF  POOLED  STATIC  AND  1%  LIFETIME  DATA 

FOR  1/5 -SCALE  SPECIMENS 


Specimen 

Test 

(KIPS) 

Rank 

(ra) 

Probability  of  Failure 

Jd3L 

n+1 

K012440 

S 

119.0 

,, 

0.941 

K900441 

n 

114.0 

15 

0.883 

K012446 

s 

111.2 

14 

0.824 

K012459 

s 

110.4 

13 

0.765 

K012439 

s 

109.8 

12 

0.706 

K905016 

1% 

109.1 

11 

0.647 

K012437 

S 

108,5 

10 

0.588 

K904958 

17c 

108.4 

9 

0.529 

K012441 

S 

106.4 

8 

0.471 

K900448 

17. 

104.4 

7 

0.412 

KOI 24 74 

S 

99.0 

6 

0.353 

K012463 

S 

98.0 

5 

0.294 

K994569 

S 

92.0 

4 

0.236 

K904955 

1% 

90.2 

3 

0.177 

K994578 

S 

89.2 

2 

0.118 

IK994588 

S 

86 . 6 

1 

0.059 

Notes:  S  =  Static  test,  1%  =  1%  lifetime 
Mean  Failiuj-,  Load  =  103.51  K 
Reference  Tables  XV  and  XVI. 
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2.  Ten  percent  lifetime  residual  strength  data, 

1/5  scale  -  Table  XXII 

3.  Fifty  percent  lifetime  residual  strength  data, 

1/5  scale  -  Table  XXIII 

4.  Lifetime  data,  1/5  scale  -  Table  XXIV 

5.  Ten  percent  lifetime  data,  1/2  scale  -  Table  XXV 

6.  Lifetime  data,  1/2  scale  -  Table  XXVI 

Residual  strength  and/or  lifetime  was  plotted  against  the  proba¬ 
bility  of  failure  on  Weibull  paper  for  the  applicable  cases. 
Straight  lines  were  drawn  through  each  batch  of  data  points 
(reference  Figures  33,  34,  and  35).  It  should  be  noted  that  the 
pooled  static  and  one  percent  lifetime  data  randomly  oscillate 
about  the  best  fit  lines.  The  high  load  tail  of  the  10  and  50 
percent  lifetime  data  is  approximately  linear;  however,  the  low 
outliers  do  not  fit  the  same  line.  This  factor  is  significant 
in  the  second  pass  data  reduction.  A  summary  of  the  least  squares 
fit  Weibull  data  is  tabulated  in  Table  XXVII. 


6.3  RESIDUAL  STRENGTH  CHARACTERIZATIONS 


Numerous  cumulative  damage  theories  have  been  developed  for 
metallic  materials.  The  design  problem  encountered  by  using  these 
various  theories,  which  include  Miner’s  rule,  is  that  a  design 
ba  ed  on  only  theory  can  consistently  expect  more  than  a  decade 
of  disagreement  between  the  theoretical  life  and  the  empirically 
observed  life.  In  addition,  application  of  a  theory  such  as  Miner’s 
rule  requires  a  very  large  data  base  for  the  material  under  con¬ 
sideration.  Such  a  data  base  is  available  for  certain  types  of 
metallic  materials,  but  only  limited  data  has  been  developed  for 
advanced  composite  materials.  Classically,  development  of  a  large 
data  base  requires  establishing  static  and  fatigue  S/N  data  for 
several  stress  concentrations  and  stress  ratios.  A  data  base 
development  approach  such  as  this  would  not  be  economically  feasible 
for  advanced  composite  materials  because  the  "tailoring”  capability 
of  the  basic  lamina  allows  for  too  many  possible  laminate  config¬ 
urations  to  characterize  in  the  classical  manner.  Thus,  either 
lack  of  funds  to  generate  a  large  data  base  or  lack  of  confidence 
in  current  cumulative  damage  theories  to  predict  life  could  be 
reasons  to  justify  investigation  of  a  new  method  for  characteriz¬ 
ing  cumulative  damage. 
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Table  XXII  RANKING  OF  10%  LIFETIME  DATA  FOR 
1/5-SCALE  SPECIMENS 


Specimen 

Residual 

Strength 

(KIPS) 

Rank 

(ra) 

Probability  of  Failure 

_£l 

n+1 

K905017 

110.6 

6 

0.857 

K905021 

108.6 

5 

0.715 

K904959 

106.0 

4 

0.571 

K904950 

105.0 

3 

0.428 

K904957 

104.6 

2 

0.286 

K900447 

79.4 

1 

0.143 

Notes:  n  =  Number  of  specimen  ranked 

Mean  Residual  Strength  =  102. 37K 
Reference  Table  XVI. 
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Table  XXIII  RANKING  OF  50%  LIFETIME  DATA  FOR 
1/5-SCALE  SPECIMENS 


Specimen 

Residual 

Strength 

(KIPS) 

Rank 

(ra) 

Probability  of  Failure 

JSQL 

n+1 

K905362 

116.7 

20 

0.954 

X900450 

115.6 

19 

0.905 

K905371 

113.0 

18 

0.857 

K905020 

111.9 

17 

0.810 

K905365 

111.0 

16 

0.762 

K905363 

110.4 

15 

0.715 

K905018 

11C.0 

14 

0.667 

K904956 

106.0 

13 

0.619 

K905366 

105.0 

12 

0.571 

K905467 

103.5 

11 

0.524 

K905022 

102.0 

10 

0.476 

K905369 

99.3 

9 

0.428 

K905015 

96.0 

8 

0.381 

X905023 

95.0 

7 

0.333 

K905469 

92.5 

6 

0.286 

K904954 

86.0 

5 

0.238 

K904952 

79.6 

4 

0.191 

K904951 

73.5 

3 

0.143 

K900449 

a54.5 

2 

— — b 

K904953 

a54 .5 

1 

— — b 

Notes:  aMaximura  load  in  spectrum 

^Specimen  failed  before  reaching  a  50%  lifetime 
n  =  Number  of  specimen  ranked 

Mean  Residual  Strength  *  101. 5K  (Excluding  K900449  &  K904953) 
Reference  Table  XVI. 
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Table  XXIV  RANKING  OF  FATIGUE -TO -FA I LURE  DATA 

FOR  1/5 -SCALE  SPECIMENS 


Specimen 

Lifetime 

Rank 

(ra) 

Probability  of  Failure 
ra 
n+1 

K905364 

5.07 

20 

0.954 

K905368 

4.06 

19 

0.905 

K905367 

3.76 

18 

0.857 

K905468 

3.66 

17 

0.810 

F504421 

3.30 

16 

0.762 

F504418 

2.41 

15 

0.715 

K994565 

2.23 

14 

0.667 

F504419 

1.43 

13 

0.619 

F504622 

1.40 

12 

0.571 

F504417 

1.27 

11 

0.524 

K012464 

1.22 

10 

0.476 

K012443 

1.01 

9 

0.428 

K900446 

0.88 

8 

0.381 

K900443 

0.82 

7 

0.333 

K905470 

0.81 

6 

0.286 

K012479 

0.73 

5 

0.238 

KO 12481 

0.73 

4 

0.191 

K012485 

0.65 

3 

0.143 

K012484 

0.46 

2 

0.095 

KOI  2480 

0.34 

1 

0.048 

Notes:  n  =  Number  of  specimen  ranked 
Mean  lifetime  =  1.812 
Reference  Table  XVII. 
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Table  XXV  RANKING  OF  10%  LIFETIME  DATA  FOR 
1/ 2-SCALE  SPECIMENS 


Specimen 

Residual 

Strength 

(KIPS) 

Rank 

(ra) 

Probability  of  Failure 

J&. 

n+1 

K012454 

240.0 

5 

0.834 

K012453 

236.5 

4 

0.668 

F504618 

219.0 

3 

0.500 

F504625 

206.0 

2 

0.333 

F504624 

168.5 

1 

0.167 

Notes:  n  *=  number  of  specimen  ranked 


Mean  residual  strength  =  214. OK 
Reference  Table  XVIII 


Table  XXVI  RANKING  OF  FAT I GUE -TO- FAILURE  DATA 
FOR  1/2-SCALE  SPECIMENS 


Specimen 

Lifetime 

Rank 

Probability  of  Failure 

—Sdl 

n+1 

F504623 

0.84 

5 

0.834 

F504622 

0.78 

4 

0.666 

F504620 

0.57 

3 

0.500 

F504619 

0.44 

2 

0.333 

F504621 

0.23 

1 

0.167 

Notes:  n  =  Number  of  specimen  ranked 

Mean  lifetime  =  0.572  Reference  Table  XIX 


Probability  of  Failure  Probability  of  Failure  Probability  of  Failure 


Residual  Strength  (KIPS) 

Figure  33  Probability  of  Failure  Versus  Residual  Strength  Plot  on  Weibull  Paper 
for  1/5  Scale  Static  and  10%  Lifetime 


Figure  34  Probability  of  Failure  Versus  Residual  Strength  and/or  Lifetime  Plot  on 
Weibull  Paper  for  1/5  Scale  50%  Lifetime  and  Fatigue  to  Failure  v 


Residual  Strength  (KIPS)  Lifetimes 

Figure  35  Probability  of  Failure  Versus  Residual  Strength  and/or  Lifetime  Plot  on 
Weibull  Paper  for  1/2  Scale  10%  Lifetime  and  Fatigue  to  Failure 
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Table  XXVII  SUMMARY  OF  WEIBULL  PARAMETERS  FITTED  TO 
1/5 -AND  1/2-SCALE  TEST  RESULTS 


Number  of 

Weibul] 

- - - 

.  Parameters*5 

Reference 

Test 

Scale 

Specimens 

a 

0C 

?s 

Table 

Static 

1/5 

16a 

11.58 

108.0 

11% 

XXI 

10%  Lifetime 

1/5 

6 

7.38 

108.9 

16% 

XXII 

507,  Lifetime 

1/5 

20 

4.94 

105.8 

237 

XXIII 

Fatigue-to- 

Failure 

1/5 

20 

1.44 

1.99 

68% 

XXIV 

10%  lifetime 

1/2 

5 

7.05 

227.6 

17% 

XXV 

Fatigue-to- 

Failure 

1/2 

5 

1.95 

0.67 

53% 

XXVI 

. 

Notes:  a  Data  pooled  for  eleven  static  and  five  17-  lifetimes, 
b  Based  on  procedure  described  in  Section  VII. 

c  See  Appendix  VII  to  account  for  the  least-of-two 
scale  effects 


The  Weibull  distribution  seems  to  fit  the  high  probability 
of  failures  rather  well;  however,  the  lower  probabilities  of 
failure  tend  to  deviate  from  the  straight  line  on  the  Weibull 
plots  (Figures  33,  34,  and  35).  The  straight  line  through  the 
data  points  on  tne  Weibull  paper  must  be  extended  beyond  the  data 
to  the  probability  of  failure  level  desired.  Extrapolated  results 
on  this  basis  have  a  degree  of  uncertainty.  If  the  data  points, 
at  the  lower  probability  of  failure  levels,  tend  to  curve  away 
from  the  vertical  axis  on  the  Weibull  plot,  then  the  straight  line 
extrapolation  at  the  lower  end  should  be  conservative.  The  curva¬ 
ture  of  the  first  failure  and  of  the  residual  strength  distribu¬ 
tions  correlates  with  the  proposed  wearout  model. 

The  slope  of  the  straight  line  through  the  data  plotted  on 
the  Weibull  paper  decreased  with  lifetime.  Since  a  decreases  as 
the  slope  decreases,  the  coefficient  of  variation  (or  scatter) 
increases.  The  slopes  of  the  straight  lines  through  the  10  per¬ 
cent  lifetime  data  on  the  Weibull  plot  for  the  1/5-  and  1/2-scale 
specimens  were  approximately  equal  (Figures  33  and  35).  The  same 
tendency  was  observed  for  the  1/5-  and  1/2-scale  fatigue-to-failure 
cases.  For  the  10  percent  lifetime  case,  a  was  7.38  and  7.05  for 
the  1/5-  and  1/2-scale  specimens,  respectively  (Table  XXVII).  In 
the  fatigue-to-failure  case,  a  was  1.44  and  1.95  for  the  1/5-  and 
1/2-scale  specimens,  respectively.  Although  a  limited  number  of 
specimens  were  used,  the  results  imply  that  a  is  independent  of 
specimen  size.  The  shape  parameter  can  be  defined  using  a  given 
number  of  small  scale  tests  and  the  same  value  can  be  applied  to 
different  specimen  sizes.  The  mean  residual  strength  or  mean 
lifetime  can  be  defined  for  the  larger  specimens  with  a  limited 
number  of  tests.  The  scaling  factor  between  the  different  speci¬ 
men  sizes  can  be  determined  by  applying  the  a  for  the  small 
specimen  to  the  larger  specimen  and  using  the  respective  mean 
values . 

The  fatigue-to-failure  data  for  the  1/5-scale  specimens  were 
pooled  and  assumed  to  represent  a  homogeneous  population.  Analy¬ 
sis  of  the  a  ita  showed  that  the  specimens  fabricated  and  tested 
before  5-21-71  and  those  fabricated  and  tested  after  11-19-71 
(after  a  six-month  program  delay)  indicated  that  the  population 
was  not  homogeneous.  The  sign  test  (a  nonparametric  test,  Refer¬ 
ence  32)  indicates  only  a  0.0107  chance  of  the  fatigue-to-failure 
data  being  from  a  single  population.  Analysis  of  each  data  set, 
assuming  a  Weibull  distribution,  yielded  the  following  results: 

1.  Specimens  fabricated  and  tested  before  5-21-71,  a f  = 

2.07  and  t  =  1.03 
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2.  Specimens  fabricated  and  tested  after  11-19-71,  ot  .  = 

1.75  and  t  =  3.14. 

The  trend  was  for  average  fatigue  lifetime  to  increase  in 
time, which  is  suggestive  of  a  learning  curve. 

Even  though  the  empirical  data  seemed  to  segregate  into  two 
groups,  the  data  were  pooled  for  final  analysis.  Two  reasons 
are  used  to  justify  the  decision: 

1.  In  the  analysis  of  laminate  data,  Reference  33,  it  was 
determined  that  an  unconservative  analysis  of  the 
specimen  variation  would  K e  obtained  if  batch- to-batch 
variation  was  not  included. 

2.  The  test  program  was  designed  to  provide  reliability 
data  relevant  to  the  production  design  of  large-scale 
bonded  joints  and  the  batch- to-batch  variation  is  a 
significant  portion  of  the  expected  production  varia¬ 
bility. 

Although  the  effect  of  the  fatigue  shape  parameter  is  sig¬ 
nificant,  the  pooled  value  of  =  1.44  is  conservative. 

All  of  the  remaining  analyses  were  executed  with  the  assump¬ 
tion  of  a  homogeneous  fatigue-to-failure  population. 


6.4  RESIDUAL  STRENGTH  -  LIFETIME  DISTRIBUTION 
BASED  UPON  PROPAGATION  OF  INITIAL  FLAWS 


It  has  been  postulated  (Reference  34)  that  flaw  growth  can 
be  characterized  in  terms  of  residual  strength  capacity;  hence, 
life  and  residual  strength  distributions  are  related  by  a  flaw- 
growth  equation  to  the  initial  static  distribution.  In  this 
format,  the  specification  of  the  reliability  of  structural  com¬ 
ponent  requires  the  definition  of  a  residual  strength-lifetime 
and  lifetime-maximum  load  relationship.  The  empirical  determina¬ 
tion  of  these  relationships  was  obtained  through  simulated  service 
testing  using  the  closed  loop,  digitally  controlled  random 
fatigue  facility. 

A  detailed  fracture  mechanics  analysis  of  the  test  specimen 
has  not  been  exercised.  All  experimental  observations  have  been 
made  with  respect  to  the  residual  strength  capacity  of  the  joints, 
Fft(t).  A  summary  of  the  1/5-scale  random  fatigue  data  is  shown 
in  Table  XXVTI.  A  significant  broadening  of  the  residual  strength 
can  be  seen. 
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The  bonded  joint  is  heterogeneous  enough  to  defy  precise 
definition  of  the  microscopic  stress  field  and  geometry  of  the 
flaw  field.  It  is  assumed  that  the  dispersion  of  filaments  in 
the  polymer  matrix  combined  with  voids  and  other  micro-defects 
may  be  represented  by  a  flaw  field  C  (where  C  is  a  symbolic 
representation  of  a  dispersion  of  flaws).  Rivlin  and  Thomas 
(Reference  35)  proposed  that  during  the  quasi-static  growth  of 
a  flaw  at  which  the  external  work  change  is  zero  that  the  change 
in  the  free  energy  of  deformation  is  given  by 


where 

W  =  strain  energy 
T  =  characteristic  tearing  energy 
t^  =  reference  thickness.  / 

The  criterion  is  similar  in  form  to  Griffith’s  (Reference  36), 
but  T  is  not  a  surface-free  energy.  The  structure  will  remain 
stable  until  the  following  inequality  is  violated 

-rh  [I!]  ■  T  <  Tc  <2> 

where  T  is  the  critical  energy  associated  with  fracture  and  is 
a  material  property.  If  the  structure  is  elastic,  a  first-order 
approximation  to  the  energy  available  to  drive  the  flaw  becomes 

T  a  CW  a  C  (3) 

K 

where 

K  =  specimen  stiffness 
F  =  applied  load. 

The  instability  condition  becomes  equivalent  to  the  Griffith  con 
dition : 

Fn  VC"  «  VW  .  (4) 

R  c 

An  engineering  approximation  to  the  fatigue  crack  growth  problem 
may  be  gained  by  extending  the  quasi-static  analysis  to  the  time 
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dependent  fatigue  problem.  Assume  that  the  rate  of  flaw  size 
growth  dC(t)/dt  is  proportional  to  a  power  function  of  the  work 
available  to  create  new  surface  T(t)  where  both  quantities  are 
(t)  time  dependent.  This  assumption  yields 


“  T(t)r  (5) 

where  r  is  an  arbitrary  constant.  If  C(t)  is  slowly  changing 
compared  to  the  random  fluctuations  of_W(t),  the  steady  work  in¬ 
crease  W(t)  may  be  considered.  This  yields  the  expression 


T(t)  a  C(t)W(t)  (6) 

and  substituting  Equation  (6)  into  Equation  5  yields 

C(t)"r  dC(t)  =  A,  W(t)r  dt  (7) 

where  A]_  is  a  constant  of  proportionality.  Equation  7  may  be 
formally  integrated  to  yield 


C(t)_r  dC(t)  =  A 


W(t)r  dt 


(8) 


where  the  subscript  (o)  defines  the  initial  conditions.  The  in¬ 
tegral  on  the  right  side  of  Equation  8  cannot  be  evaluated  analy¬ 
tically.  The  test  spectrum  F(t)  was  developed  as  a  random  process 
that  preserved  the  expected  cumulative  statistics  for  the  aircraft 
studied.  Each  mission  segment  was  simulated  with  bursts  of 
random  loads  so  that  each  half  cycle  was  statistically  independent 
but  asymptotically  integrated  in  time  yielding  the  desired  cumu¬ 
lative  exceedances.  The  resulting  history  F(t)  is  a  stationary 
random  process;  hence, 


W(t)r  at 


(9) 


(10) 
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where  the  constant  A2  is  to  be  experimentally  determined.  The 
left  side  of  Equation  7  may  be  integrated  to  yield 

CoI_r+11  '  ctI"r+11  =  A1A2  lr-U  lt-t0]  (11) 

and  -  ince  A-^  and  A2  are  arbitrary  constants,  they  may  be  combined 
to  yield 

c0['r+1)  -  Ct('r+l!  =  A3[r-1]  (12) 


Equation  12  may  be  converted  to  a  residual  strength  form  by 
using  Equation  4.  Since 


Fd  VC~  a  VTcT  =  B, 

K  C  x. 


(13) 


Equation  12  may  be  rearranged  to  yield 

FR(°)2[r"1]  -  FR(t)2[r'1]  =  A3B12[r"1]  [r-1]  [t-tQl  (14) 
and  the  constants  A3  and  B^2^r  ^  may  be  combined  to  yield 


FR(o)2[r'1]  -  FR(t)2[r"1]  =  A4[r-1]  [t-t0]. 


(15) 


The  initial  (t  =  0.0)  residual  strength  distribution  is 
assumed  to  be  Weibull  (Reference  37) 


F 


P(FR(o)  >  Fr)  =  exp 


R 


:fr(o)1 


(16) 


where  a0  is  the  shape  parameter,  Fr(o)  is  the  scale  parameter  at 
time  t0,  and  Fr  is  the  random  variable.  To  produce  a  time-de¬ 
pendent  residual  strength  function,  it  is  noted  that  the  set  of 
values  Fr ( t )  >  FR  is  the  same  as  the  set  of  values  (Equation  15) 
where  1 

Fr(°)  >  {FR2tr"1]  +  A4[r-1]  [ t-tQJ }  [  .  (17) 
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Hence,  the  probability  that  FR(t)  >  Fr  is  just  the  probability 
that  Equation  17  holds.  Equation  16  gives  this  probability  as 

r  1 

P(FR(t)  >  FR)  -  P  Fr(o)  >  {FR2[r“1]  -1  A4[r-1]  [ t-t0i}  2[r'lJ 


=  exp 


_Fr_|  2!r_11  [~A^[r-l]  It-t,]]  _  f 

Fr(o)  FR(o)2[r"13 


(18) 


where 


af  =  HTTT  * 


(19) 


It  can  be  established  that  Equation  18  will  be  very  close  to 
a  Weibull  distribution  for  fatigue  lifetime  (t)  if  fatigue  failures 
are  obtained  at  a  load  level  Fjru  (the  truncation  load  in  the  fa¬ 
tigue  spectrum)  that  is  of  the  order  FR(o)/2  or  less.  By  assign¬ 
ing  a  specific  value  F^ru  for  Fr  in  Equation  18  and  rearranging 
the  expression 


P(FR(t)  >  Ftru)  =  exp 


t 

2 FT 


FR(o)2[r 

A4[r-1] 


V°) 


2[r-l] 


A4[ r-1]  t0 


the  desired  result  may  be  obtained.  Let  Equation  20  be  written 
in  the  form 

P(FR(t)  >  Ftru)  =  exp  -  j—  +  ERROR)  (21) 

( t5  ) 


where 


cb  " 


A  2f  r-1 1 
FR(o)Zlr 

A4  [  r  - 1  ] 


and  is  the  approximate  characteristic  failure  lifetime  for  the 
baseline  load  spectrum  F(t). 
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Under  the  conditions  that 


(1)  t  =  0  /. 


A4[r-1]  [o]  _ 


s  0 


FTRU  1 

*»<•>  2 


and  r  >  4 


the  error  term  in  Equation  21  would  be  less  than  0.0156,  which 
would  yield  a  mode  location  error  of  less  than  two  percent.  Thus, 
for  bonded  joints  which  exhibit  substantial  wearout,  the  lifetime 
distribution  using  Equation  20  would  be  essentially  Weibull  with 
parameters  and  1;^. 

Correlation  of  the  wearout  model  with  the  experimental  data 
follows  in  the  next  section. 


6.5  RANDOM  FATIGUE  DATA  EVALUATION 


A  reduction  and  evaluation  of  the  1/5-scale  specimen  wearout 
data  and  the  1/2-scale  specimen  scaling  experiment  are  presented 
in  the  following  paragraphs. 


6.5.1  One-Fifth  Scale  Wearout  Data 


A  qualitative  evaluation  of  the  character  of  the  wearout  data 
can  be  made  by  observing  Table  XXVII.  The  important  observations 
are  that  the  initial  static  distribution  broadened  significantly 
at  the  50  percent  lifetime  trace  while  only  a  small  change  in 
characteristic  strength  was  noted.  This  behavior  is  in  agreement 
with  the  assumptions  made  in  development  of  the  wearout  model. 

A  wearout  model  was  evaluated  by  global  fitting  the  pooled 
static  and  1  percent  lifetime  data,  the  10  percent  lifetime  data, 
and  50  percent  lifetime  data  and  the  lifetime  distribution  data. 
The  fit  was  accomplished  by  forming  an  error  function  of  the  form 


ERR  =  £  ( i FE ( t )  -  i_FR( t) ) ' 

i=l 
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where 


n  =  total  number  of  data  points 
FE(t)  =  experimental  value 

and  where  the  summation  was  over  the  number  of  points  in  the  ex¬ 
perimentally  established  data  set.  The  error  function  was  mini¬ 
mized  using  the  technique  in  Reference  38. 

A  summary  of  the  best  fit  results  are  shown  on  Table  XXVIII. 


The  computed  versus  observed  residual  strength  and  lifetime  dis¬ 
tributions,  as  determined  using  Equation  24, are  summarized  in 
Figure  36.  The  correlation  of  the  observed  distributions  yields 
material-related  parameters  (F^(o),  ac,  r,  A/t)  that  allow  screen¬ 
ing  in  terms  of  physical  properties.  The  observed  growth  rate  r 
is  typical  of  polymer  materials  and  indicates  that  the  adhesive 
or  adherand  polymer  constituents  are  dominating  the  joint  relia¬ 
bility  characteristics.  If  a  tighter  lifetime  distribution  is 
desired,  the  available  control  parameter  is  u  . 
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j  In  order  to  further  explore  the  meaning  of  the  data  presented,  | 

-  it  is  of  interest  to  determine  operational  load  levels  associated  | 

|  with  various  reliabilities  for  the  same  scale  bonded  joint.  The  ] 

l,  baseline  spectrum  was  defined  in  terms  of  a  random  history  F(t)  j 

l  truncated  at  a  load  level  Ftru*  The  information  to  be  derived  is 

\  the  lifetime  capability  of  the  joint  under  reduced  load.  No  ex¬ 

perimental  data  is  available;  therefore,  the  form  of  the  evaluation  • 
used  was  based  on  metal  fatigue  behavior.  In  the  range  of  cycles 
;  to  failure  from  10^  to  10^,  the  mean  life  behavior  (or  correspond¬ 

ingly  c)  is  linear  in  ln(life)  vs  ln(amplitude)  coordinates  for 
various  metals  (Reference  3).  For  aluminum,  the  ratio  of  the 
ln(life)  vs  ln(amplitude)  function  is  -2r. 

Consider  the  transformation  (Reference  39)  of  the  reference 
spectrum  (the  spectrum  used  under  this  program)  F(t)  by  the  func¬ 
tion, 

G(t)  =  B2F(t)  (25) 

and  since  the  spectrum  is  truncated  at  F-j-ru*  The  truncation  load 
for  the  transformed  spectrum  is 

GTRU  =  B2  FTRU’  ^ 


The  truncationAload  and  characteristic  life  are  related  as  fol¬ 
lows.  The  In  t  vs  In  G^ru  is  assumed  linear  with  a  slope  of  (-2r). 
The  In(life)  vs  ln(amplitude)  function  will  be 

In  t  =  2r  In  B2  +  constant  (27) 

where  the  assumed  relation  is  fixed  at  the  known  data  derived 
point  B2  =  1  and  t  = 


From  equations  22,  26,  and  27 


t 


£R(0)2[r~1] 

“  A4[r-1] 


ftru 

.gtru. 


(28) 


Utilizing  the  approximate  lifetime  relatio^  given  in  Equation 
21  (assuming  that  ERROR  =0.0  and  generalizing  tfc>  to  t)  together 
with  Equation  28  the  safe  load  capacity  of  the  specimen  can  be 
computed.  Two  cases  were  computed:  the  load  level  which  yields 
a  characteristic  life  of  4.0  lifetimes  and  the  load  level  which 
yields  a  0.99  probability  of  achieving  4.0  lifetimes.  The  re¬ 
sults  are  summarized  in  Table  XXIX,  where  the  probability  of 
surviving  one  lifetime  is  Ps  for  t=l  and  the  required  static 
safety  factor  (F.S.)  are  shown. 
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The  difference  in  the  estimate  of  t  =  2.03  for  the  baseline 
data  in  Table  XXIX  and  the  experimental  value  of  ^  =  1.99  in 
Table  XXVII  is  due  to  the  influence  of  the  residual  strength  data 
in  the  global  fit.  The  value  given  in  Table  XXVII  is  the  result 
of  a  one -dimensional  fit  only. 


Table  XXIX  ESTIMATED  JOINT  FATIGUE  PERFORMANCE 


Case 

Gtru 

KIPS 

£ 

Lifetimes 

Ps  for  t  =  1.0 
Dimensionless 

F.S.(1) 

Baseline  Data  (.368 
surviving  t  =  2.03) 

54.5 

2.03 

0.637 

1.54 

0.368  surviving 
t  =  4.0 

51.5** 

4.00 

0.816 

1.72 

0.99  surviving 
t  =  4.0 

37.0** 

220.4*** 

0.999 

2.27 

102.7  _  fl(static) 

(1.22)GTRU  Limit  Load 


**Equation  28  A 

***Equation  21,  generalized  for  t. 


6.5.2  Scale  Effects 


Size  effects  have  been  observed  in  metals  and  glass  (Refer- 
ences40  through  43).  Small-scale  specimens  generally  exhibit 
greater  strength  per  volume  of  material  and  longer  fatigue  lives 
than  their  full-size  counterparts.  The  small-scale  specimens 
tend  to  have  longer  fatigue  lives  than  full-size  specimens  al¬ 
though  both  are  subjected  to  the  same  stress  history. 

Composite  bonded  joint  specimens  also  exhibit  size  effects. 
The  random  load  fatigue  test  results  for  the  1/5-  and  1/2-scale 
specimens  were  used  to  evaluate  the  scaling  effect. 

Economically,  it  is  desirable  to  characterize  the  allowable 
static  and  fatigue  strength  of  full-size  joints  using  small-scale 
specimens.  To  achieve  a  realistic  characterization,  the  size 
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relationship  (or  scaling  effect)  between  tht  small-  and  full-size 
specimen  should  be  established  using  a  realistic  random  load 
history  and  compatible  environment.  If  small-scale  specimens  can 
be  used  to  develop  realistic  design  allowables,  then  a  practical 
design  methodology  can  be  developed. 

The  scaling  effect  was  evaluated  using  the  1/5-  and  1/2-scale 
specimen  test  results  for  the  10  percent  lifetime  and  the  fatigue- 
to-failure  cases.  Scaling  factors  (Sp)  were  computed  for  these 
cases  four  different  ways,  i.e.. 


1. 

sF  =  («1/5/^1/2) 

x  Cp 

(29) 

2. 

Sf  "  <-^l/5/^l/2) 

x  Cp 

(30) 

3. 

SF  =  (^1/5/^i/2) 

x  Cp 

(31) 

4. 

Sf  -  ^uAl2> 

x  cp 

(32) 

where 

It 1/ c,  /o  =  Mean  values  for  the  1/5-  and  1/2-scale  speci¬ 

men,  respectively 

#l/c,  01/2  =  Characteristic  values  for  1/5-  and  1/2-scale 

specimen  respectively  based  on  least  squares  fit 

/?l/2  =  Computed  characteristic  value 

(U\l2/fi l/5>  xy3J/5 
baj'r  :3  on  8^5  =  tt  \j2- 

^1/2  =  GomPute<^  characteristic  valueAbased  on  the 

'  maximum  likelihood  estimator,  /?  „ 

^MLE 

1 


base  on  a  1/5  =  a  1/2  with  -  residual  strength 
or  lifetime. 
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The  1/5-  and  1/2-scale  specimens  were  both  5  inches  wide  with 
a  thickness  of  40  and  94  plies,  respectively  (Figures  5  and  6). 
Since  the  specimen  widths  were  the  same,  the  mean  or  characteris¬ 
tic  residual  strengths  were  computed  on  a  per-ply  basis  to  put 
the  results  on  the  same  baseline.  The  ratio  of  the  ply  thickness 
for  the  1/2-  and  1/5-scale  specimens,  i.e.  94/40,  is  equal  to 
2.35.  This  multiplying  factor  converts  to  the  residual  strength 
data  for  the  1/5-  and  1/2-scale  specimens  to  the  same  baseline 
and  is  defined  as  Cp.  Mean  and  characteristic  lifetimes  for  the 
1/5-  and  1/2-scale  specimens  are  directly  comparable;  therefore, 

Cp  is  unity  for  fatigue-to-failure  comparisons. 

The  parameters  Ui  and  y3^  (i  =  1/5  and  1/2),  for  computing 
Sp,  were  taken  from  Tables  XXI,  XXIV,  XXVI,  and  XXVII.  Details 
are  given  in  the  following  sections  for  defining  [2  anc*  ^1/2* 

The  scaling  factors  were  computed  using  Equations  29  through  32. 
The  results  are  summarized  in  Table  XXX. 

6.5. 2.1  Definition  of  @i/2 

&1/2  was  computed  using  #1/5,  #l/2>  and  &1/2'  This  approach 
was  considered  to  determine  if  P\_/2  could  be  reasonably  estimated 
using  only  a  few  1/2-scale  specimens  to  define  lx\j2  and  a  larger 
number  of  1/5-scale  specimens  to  define  ^1/5  and  ^1/5.  If  feasi¬ 
ble,  this  could  reduce  the  number  of  1/2-scale  specimens  required 
to  evaluate  Sp  using  characteristic  values. 

The  mean  value  (H)  may  be  expressed  in  terms  of  /?  and  the 
gamma  function  (/")  as  follows: 


"i  ■  V  faT  +  l) 

If  a  1/5  and  a  1/2  are  equal, 


(i  =  scale  size). 


then 


"l/5 

^1/2 


£l/5 

*1/2  ’ 


Solving  for/9^/2  and  using  the  notation  /3  -^j2  t0  distinguish  the 
yS  from  the  least  squares  type, 


*1/2  =  x  ^1/5. 

a/5 

/?]/2  was  used  in  Equation  31  to  compute  Sp.  The  results  are 
summarized  in  Table  XXX. 
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Table  XXX  SUMMARY  OF  SCALING  FACTORS  BASED  ON 
1/5- and  1/2 -SCALE  SPECIMENS 
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6. 5.2.2  Definition  of  P-^/2 


If  the  a  value  is  known  for  a  small-scale  specimen,  then  the 
value  can  be  computed  for  larger  specimens  using  the  maximum 
likelihood  estimator,  i.e., 


1 


where  :  2  ~  al/5  (assumec0 

n  =  number  of  specimen 

x^  -  residual  strength  and/or  lifetime  for  given  specimen. 

A 

P 1/2  values  were  computed  for  the  10  percent  lifetime  and  the 
fatigue-to-failure  cases  using  0^/5  equal  7.38  and  1.44,  respec¬ 
tively  (Table  XXVII).  Scaling  factors  were  computed  using  P1/2 
in  Equation  32.  The  results  are  summarized  in  Table  XXX. 

6. 5. 2. 3  Discussion 

The  results  shown  in  Table  XXX  indicate  that  there  is  a 
scaling  effect  between  different  size  bonded  composite  joint  speci 
mens.  This  effect  is  significant  and  must  be  given  careful  con¬ 
sideration  in  developing  design  allowables. 

The  results  show  that  the  scaling  factors  increase  with  life¬ 
time.  In  fact,  the  scaling  factor  for  the  fatigue-to-failure 
case  was  about  three  times  larger  than  that  for  the  10  percent 
lifetime  case.  This  implies  that  the  governing  scaling  factor 
depends  on  the  fatigue  life  requirement  rather  than  the  static 
strength  requirement. 

The  scaling  factors  were  based  on  1/5-  and  1/2-scale  speci¬ 
mens.  The  lifetime  scaling  factor  for  a  full-size  specimen  was 
not  estimated  since  a  linear  extrapolation  would  have  been  re¬ 
quired.  Without  a  sufficient  number  of  full-size  specimens  tests, 
at  least  three  different  specimen  sizes  should  be  used  to  credibly 
extrapolate  the  scaling  factor  to  a  full-size  specimen. 

The  computed  scaling  factors  were  the  same  order  of  magni¬ 
tude  for  the  respective  cases  irrespective  of  the  method  used. 

The  two  largest  scaling  factors,  3.18  and  3.36  (Table  XXX),  were 
based  on  the  ratio  of  the  mean  lifetimes  and  the  ratio  p  1 75/^1  / 2 * 
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respectively.  Both  methods  should  be  used  to  compute  Sp.  The 
results  should  be  compared  and  the  larger  scaling  factor  should 
be  used  for  design. 

The  results  of  this  investigation  indicate  a  reasonable  po¬ 
tential  for  using  small-scale  specimens  to  characterize  the  life¬ 
time  properties  of  full-size  specimens.  The  scaling  factor  is 
a  definable  quantity  and  it  should  be  evaluated  using  a  realistic 
random  load  history  and  compatible  environment. 

6.6  JOINT  FAILURE  MODES 


Joint  failures  were  analyzed  by 

1.  Visual  inspection  of  the  broken  specimen 

2.  Microscopic  inspection  of  the  failure  surfaces 

3.  Visual  inspection  (by  naked  eye  and  by  magnifying  glass) 
of  photographs  of  broken  specimens  (Appendix  IV  and  V). 

The  test  types  studied  were 

1.  Static 

2.  Constant  Amplitude 

3.  Random  Fatigue 

a.  Residual  Strength  (1,  10,  50%  lifetime) 

b.  Fatigue-to-Failure. 

The  three  failure  modes  observed  were  (1)  laminate  failure, 
(2)  bond  failure,  and  (3)  complex  failure  (laminate  and  bond). 
Specimen  K904950  (Figure  86)  characterizes  the  laminate  failure 
mode,  specimen  K905023  (Figure  101)  characterizes  the  bond  fail¬ 
ure  mode,  and  specimen  K012476  characterizes  the  complex  mode. 
Different  failure  modes  were  observed  for  the  same  test  cases 
(both  1/5  and  1/2  scale).  For  this  reason, it  is  difficult  to 
confidently  relate  the  failure  modes  to  a  particular  type  of 
loading.  Generally,  however,  it  appears  that  the  laminate 
failure  mode  is  predominate  for  the  static  cases  and  that  the 
progressive  bond  failure  mode  prevails  for  the  fatigue  specimen. 
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6.7  EFFECT  OF  SAMPLE  SIZE 


The  question  treated  in  this  section  is:  what  size  sample 
should  be  taken  to  adequately  characterize  various  distribution 
parameters  associated  with  a  particular  element  or  component? 

Extremely  large  sample  sizes  are  usually  limited  in  number 
due  to  fabrication  and  testing  costs.  Without  doing  a  complete 
trade-off  study  between  the  cost  of  various  sample  sizes  and  the 
value  of  the  statistical  information  gained  from  such  samples, 
the  problem  is  looked  at  in  a  qualitative  manner.  This  allows 
for  selection  of  a  sample  size  that  is  relatively  efficient  in 
maximizing  statistical  information  gain  for  each  test  dollar  spent. 

In  general,  the  value  of  statistical  information  gained  does 
not  increase  linearly  with  the  sample  size.  Therefore,  dollar 
wise,  a  point  of  diminishing  returns  is  reached  by  increasing  the 
sample  size.  A  qualitative  measure  of  the  amount  of  statistical 
information  contained  in  a  sample  for  estimating  a  distribution 
parameter  is  the  number  that  the  point  estimate  (estimate  base 
only  on  observed  values)  must  be  multiplied  by  in  order  to  obtain 
a  new  estimate  with  an  associated  high  degree  of  confidence. 


6.7.1  Sample  Size  Effects  on  Establishing  Shape 
and  Scale  Parameters  Based  Upon  the  Par¬ 
ticular  Results  of  This  Program 


The  number  of  specimens  required  to  achieve  an  equivalent 
confidence  in  the  estimate  of  the  Weibull  shape  and  scale  param¬ 
eters  is  a  key  issue  in  the  design  of  a  cascaded  test  program. 

As  shown  in  Figure  37,  the  penalty  paid  to  obtain  a  90-percent 
confidence  lower  bound  on  the  shape  parameter  is  only  a  function 
of  sample  size.  The  improvement  in  the  estimate  is  steep  through 
17  specimens  where  the  slope  of  the  curve  dfqo/dn  is  greater  than 
1.0.  The  return  diminishes  rapidly  above  17  specimens. 

The  number  of  specimens  required  to  achieve  equivalent  con¬ 
fidence  in  the  Weibull  scale  parameter  for  the  case  of  a  90- 
percent  confidence  lower  bound  for  an  estimated  shape  parameter 
is  shown  in  Figure  38  and  known  shape  parameter  is  shown  in  Figure 
39.  The  number  of  specimens  required  to  estimate  the  scale 
parameter  is  a  function  of  the  shape  parameter. 

Examples  of  the  data  requirements  can  be  developed  from  the 
specific  data  obtained  under  this  program. 


102 


}' 

\ 

l 

l 

i 

l’ 

'• 

l 

l 

i 

t 

c 

E 

5 


f- 

\ 

V 


^f*gff^gar, y «^«*??**  r^nr  4 


\v»w>.j»»n.  ■.lllH,*(jgj 

h  WwWJMWW»W»Ht1«»*MTOI»aB«WOffiMHi 


1.0 


.9 


.8 


f90  -7 


.6 


.5 


.4 

0  5  10  20  30  40  50  60  70  80  90  100 

Number  of  Test  Specimens 

Figure  38  Factors  fgo  for  a  90%-Confldence  Lower  Bound  for  the  Welbull  Scale  Parameter 
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The  apparent  shape  parameter  for  the  static  data  shown  in 
Column  1  of  Table  XXXI  was  a0  =  10.6.  From  Figure  37  the  approp¬ 
riate  penalty  in  £90  to  achieve  a  90-percent  confidence  estimate 
of  «0  is  f90  =  0.66.  For  this  sample  size  the  corresponding  con¬ 
fidence  level  is  obtained  for  the  scale  parameter  with  a  penalty 
of  fgo  =  0.92.  Note  that  the  penalty  required  in  scale  parameter 
would  be  only  0.87  if  only  five  specimens  were  used.  For  high 
shape  parameters,  the  expected  convergence  to  a  mean  is  rapid  and 
can  be  heuristically  observed  in  the  data  presented  in  Column  1 
of  Table  XXXI.  The  data  set  converges  to  the  mean  within  the 
first  5  specimens  and  then  simply  oscillates  randomly  about  the 
mean  value.  The  50-percent  lifetime  data  was  the  slowest  to  con¬ 
verge,  and  as  noted  in  Table  XXVII,  it  had  a  shape  parameter  of 
4.94.  At  a  shape  parameter  of  4.94,  20  specimens  would  be  re¬ 
quired  to  achieve  the  same  design  penalty  (fqQ  =  .94),  and  two 
specimens  would  be  required  at  a  shape  parameter  of  10,  Figure  40. 

In  subsection  6.4,  it  was  noted  that  the  homogeneity  of  the 
fatigue-to-failure  population  could  not  be  established.  The  shape 
parameter  of  the  data  subsets  were  2.07  and  1.75  indicating  that 
a  shape  parameter  for  the  type  of  joint  probably  is  not  much  larger 
than  2.0.  Over  100  specimens  would  be  required  to  achieve  a 
penalty  of  fgg  =  0.94  at  a  shape  parameter  of  2.0,  Figure  39. 

This  problem  is  reflected  in  the  slow  convergence  (or  lack  of 
convergence)  in  the  fatigue-to-failure  data. 

A  summary  ot  the  observations  made  from  the  test  data  follows. 

1.  A  larger  sample  size  is  required  to  establish  the  mean 
lifetime  than  it  is  to  establish  the  mean  static  and 
residual  strength. 

2.  At  least  5  specimens  should  be  used  to  define  the  mean 
static  strength. 

3.  At  least  100  specimens  should  be  used  to  achieve  an  esti¬ 
mate  of  mean  fatigue  life  that  would  be  consistent  with 
the  static  strength  observations. 


Fifteen  specimens  appears  to  be  a  reasonable  number  to 
define  the  mean  residual  strength  for  the  50  percent 
lifetime  case  (1/5  scale). 


The  number  of  specimens  used  to  define  the  mean,  a  , 
coefficient  of  variation,  and  factor  of  safety  must  be 
consistent  with  the  confidence  limit  factor  desired. 

The  confidence  level  factor  increases  as  the  number  of 
specimen  increases  but  so  does  the  cost.  The  number  of 
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Table  XXXI  SUMMARY  OF  MEAN  VALUES  FOR  VARIOUS  SAMPLE  SIZES 
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In  same  order  tests  were  conducted.  Specimen  failed  before  reaching  50%  lifetime. 
Mean  values  based  on  only  those  specimen  reading  a  50%.  lifetime. 

(ti=k)  denotes  the  total  number  of  k  specimens  in  the  sample. 


specimens  required  to  do  a  credible  job  involves  a  trade¬ 
off  between  the  confidence  level  and  the  cost. 

6.  The  results  of  Table  XXXII  indicate  that  the  confidence 
level  factor  is  less  for  the  fatigue-to-faiJure  case 
than  for  the  static  or  other  residual  strength  cases. 

The  factor  of  safety  required  is  governed  by  the  fatigue- 
to-failure  requirement;  yet,  the  results  show  that  the 
confidence  level  was  less  for  this  case  (for  a  given 
sample  size)  than  the  other  cases.  Therefore,  careful 
attention  must  be  given  to  the  sample  sizes  and  corres¬ 
ponding  degree  of  confidence  in  the  statistical  param¬ 
eters  . 

A  summary  of  fqg  confidence  lower  bound  factors  for  the 
Weibull  scale  parameter  (assuming  a's  are  known)  with  respect  to 
sample  size  are  given  in  Table  XXXIII  for  1/5-  and  1/2-scale 
specimens.  The  fgo  factors  in  this  table  were  estimated  from 
Figure  39  assuming  the  computed  a  values  in  the  table  were  con¬ 
fidently  known. 

In  the  1/5-scale  fatigue- to-failure  case,  a  maximum  likeli¬ 
hood  value  of  1.44  was  computed  for  the  Weibull  shape  parameter 
(«  )  based  on  20  tests  (reference  Table  XXXIII).  The  90-percent 
confidence  lower  bound  for  a  ,  based  on  20  tests,  is  approximately 
0.74.  This  means  that  there  is  a  90-percent  confidence  that  a 
will  be  no  smaller  than  0.74  x  1.44  or  1.06.  The  coefficient 
of  variation  corresponding  to  a  of  1.44  and  1.06  is  approximately 
0.68  and  0.93,  respectively  (see  Figure  40  and  Appendix  VI). 

From  Figure  40,  it  is  seen  that  the  coefficient  of  variation  in¬ 
creases  as  a  decreases.  The  larger  the  coefficient  of  variation, 
the  larger  the  factor  of  safety  must  be  to  meet  the  reliability 
goal. 


6.8  EFFECT  OF  RESIDUAL  STRENGTH  AND  FATIGUE  LIFE  VARIA¬ 
BILITY  ON  REQUIRED  DESIGN  SAFETY  FACTORS 


The  effect  of  residual  strength  and  fatigue  life  variability 
on  the  safety  factors  required  for  design  are  dealt  with  in  the 
following  paragraphs. 

Safety  factors  are  commonly  used  to  account  for  the  varia¬ 
tion  in  residual  strengths  and  fatigue  lives.  The  results  of 
this  program  indicate  that  the  distribution  of  residual  strength 
broadens  with  lifetime.  This  implies  that  static  design  safety 
factors  for  bonded  joints  should  be  projected  to  the  one  lifetime 
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Table  XXXII  SUMMARY  OF  MEAN  VALUES  FOR  VARIOUS 

SAMPLE  SIZES  (1/2 -Scale  Specimen  Tests) 


Number  of 
Specimens  a 

1 

2 

3 

4 

5 

Sample  Mean 
Reference  Table 


>0  Lifetime 
(n-5) 

(KIPS) 

1  Fatigue-to- 
Failure  (n=5) 
Lifetime 

236.5 

0.57 

238.3 

0.68 

215.0 

0.59 

212.7 

0.51 

214.0 

l 

0.57 

214.0 

0.57 

XVIII 

XIX 

Notes : 

a  In  same  order  tests  were  conducted 

(n=k)  denotes  the  total  number  of  k  specimen 
in  the  sample 
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residual  strength  distribution  since  it  will  require  the  greatest 
safety  factor. 

The  flaw  growth  rate  exponent  and  the  static  shape  parameter 
affect  the  fatigue  lifetime  distribution  for  bonded  joints.  When 
the  primary  load  transfer  mechanism  is  through  the  adhesive  bond, 
the  flaw  growth  rate  exponent  is  essentially  constant.  The  static 
shape  parameter  can  be  controlled.  If  the  shape  parameter  is 
doubled,  the  resulting  joint  scatter  is  equivalent  to  notched 
aluminum  (Reference  3). 

The  mean  life  of  a  specimen  was  strongly  influenced  by  the 
specimen  scale.  The  scale  effect  has  not  been  analytically 
treated;  however,  it  appears  to  invalidate  the  concept  of  a  stress 
analysis  based  design  that  projects  large-scale  designs  on  the 
basis  of  experimental  data  gained  for  small  coupons.  Full-scale 
testing  of  major  bonded  joints  will  be  mandatory. 

Although  the  fatigue  scatter  factor  requirement  will  drive 
the  design  in  these  instances,  the  designer  must  have  this  param¬ 
eter  transformed  into  a  static  strength  requirement  in  order  to 
actually  design  the  part.  In  transforming  the  fatigue  scatter 
factor  into  a  static  safety  factor  that  accounts  for  fatigue, 
it  can  be  seen  that  this  safety  factor  will  depend  on  several 
statistical  considerations  including  the  materials  fatigue  coef¬ 
ficient  of  variation,  anticipated  random  environment,  the  relia¬ 
bility  goal  established  for  a  specific  lifetime,  and  the  statis¬ 
tical  relationship  of  time  and  residual  static  strength. 
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SECTION  VII 

RELIABILITY  ASPECTS  OF 
BONDED  JOINTS 

7 . 1  INTRODUCTION 


1 


The  reliability  aspects  of  boron-to- titanium  bonded  scarf 
joints  are  discussed  in  this  section.  The  concept  of  a  reliability- 
based  joint  design  and  the  data  requirements  are  also  considered. 


7 . 2  STRENGTH-TIME  RELATIONSHIPS 


The  primary  objective  of  this  program  was  to  develop  an  un¬ 
derstanding  of  the  material-related  parameters  affecting  the 
reliability  characteristics  of  a  large-scale  bonded  joint.  For 
this  study,  the  expected  service  history  was  assumed  to  be  constant 
within  the  fleet  and  described  by  a  set  of  mission-related  ex¬ 
ceedance  curves.  The  reliability  model  evaluation  reduces  to  the 
determination  of  G(x|t),  the  cumulative  distribution  of  loads,  and 
FR(t)  the  static  residual  strength  function. 

Tne  service  life  reliability  of  a  component  is  derived  here 
based  on  an  arbitrary  operational  load  distribution. 

Let  G(x|t)  denote  the  cumulative  distribution  of  environmental 
loads  at  tine  t,  i.e.,  G(xjt)  is  che  probability  that  a  single  load 
occurrence  will  have  magnitude  less  than  or  equal  to  x  at  time  t. 
Let  FR(t)  denote  the  cumulative  residual  static  strength  determina¬ 
tion  of  the  component  after  a  service  life  of  length  t.  Let  ti, 
t£,  .  .  .  ,  tfc,  ...  be  the  specific  times  for  which  probability 
of  survival  (reliability)  is  to  be  computed.  Also,  let  (f(F^(tjc)) 
denote  the  residual  strength  probability  density  function  charac¬ 
terized  by  a  set  of  time-dependent  parameters  and  assume  (for  com¬ 
putational  convenience)  that  this  function  is  fixed  over  each  small 
time  interval  (tfc-i,  tk) . 


g(x  =  FR(tk))  =  |  G(x  =  FR(tk)f  N 


(29) 
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which  is  the  probability  that  all  loads  in  the  block  of  time 
(t^-i,  t^)  are  less  than  the  value  Fr,  where  N  is  the  expected 
number  of  load  cycles  in  the  time  interval  (t^-i,  t^) .  Then  L(tfc) 
defined  by 

00 

L(tk)  =  f  g(x  =  FR(tk))f(FR(tk))  d  Fr  (30) 

o 

is  the  probability  that  during  the  time  interval  (t^.^,  tO  a  load 
will  not  exceed  the  residual  strength  of  the  component  ana  cause 
failure.  Finally,  let  R( t^)  be  the  probability  that  the  structural 
component  survives  each  time  interval  up  to  and  including  the 
interval  (t^-i,  tk) ;  then,  R(tk)  can  be  written  in  terms  of  L(tk) 
as  follows : 


R(t1)  =  L(tx) 

R(t2)  =  L(t1)  L(t2) 


(31) 


k 

R(t  )  =  n  L(t  ). 

k  i=l  L 

Equation  31  represents  the  service-life  reliability  as  a  function 
of  time. 


7.3  RESIDUAL  STRENGTH  FUNCTION 


A  positive  correlation  with  a  residual  strength  model  has 
been  obtained.  The  model  development  is  presented  in  subsection 
6.5.  The  model  defines  the  residual  strength  function  Fft(t).  The 
experimentally  observed  residual  strength  function  has  a  monotoni- 
cally  decreasing  residual  strength  and  increasing  coefficient  of 
variation  on  residual  strength  in  time.  The  parameter  estimation 
techniques  used  to  fit  the  experimental  data  at  strength  or  time 
traces  follows. 


112 


7.3.1  Maximum-Likelihood  Estimators 


The  density  function  of  the  two-parameter  Weibull  distribution 
with  shape  parameter  a  and  scale  parameter  b  is  given  by 


f(x  |  a,b)  =  |  (^) 


exp  (-  r  )  . 


(32) 


The  likelihood  function  for  a  sample  of  n  observation  is 


n  n  n  X,*  a-1  .  n  x.  a  v 

L  -  n  f(x.  )  a,b)  =  (§)  n  0  exp  (  -  £  0  )  .  (33) 
i=l  i=l  \  i=l  / 

Taking  the  logarithm  of  (33)  and  then  differentiating  with  respect 
to  b  and  a  in  turn  and  equating  to  zero,  the  following  equations 
are  obtained: 


(34) 


dln(L) 

da 


~  -  n  In  b  +  Y  In  x,.  -  YL  (  £)  =  0. 

3  1=1  1  i=l  b 


(35) 


The  maximum- likelihood  estimators,  a  and  b  ,  for  the  shape  and 
scale  parameters  a  and  b  are  obtained  by  solving  Equations  (34)  and 
(35)  for  a  and  b  yielding 

■  "  £  1/d 

/  \ 


"E  xi 


nE  .Xj  In  xL  -  E  x?  E  ln(xt: 


where  x^,  i-1,  2,...,n,  represent  a  random  sample  of  size  n.  Equa¬ 
tion  (37)  is  solved  by  an  iterative  procedure  using  a  "guess1'  as 
the  initial  value  of  cL  Equation  (36)  is  then  solved  using  the  a 
obtained  from  the  solution  of  Equation  (37).  An  IBM  System  370 
Computer  has  been  programmed  to  solve  Equations  (36)  and  (37).  The 
computer  program  is  designated  as  TW9. 

The  results  of  Reference  44  can  be  used  to  estimate  the  shape 
and  scale  parameters  a  and  b  at  yarious  confidence  levels  (see 
Reference  3).  The  random  variables  —  and  a  ln(^/b)  are  distributed 
independently  of  the  parameters  a  and  b  ;  thus,  the  distribution 
of  these  random  variables  is  the  same  regardless  of  the  "true" 
underlying  values  of  a  and  b  . 

Thoman  et  al.,  Reference  44,  have  empirically  evaluated  by 

Monte  Carlo  techniques  the  distributions  of  ^/a  and  S  ln(^/b) 
for  various  sample  sizes.  Using  these  distributions,  confidence 
intervals  for  the  parameters  a  and  b  are  easily  obtainable. 
Confidence  factors,  fjQ,.  f°r  a  90-percent  confidence  lower  bound 
on  the  shape  parameter  a  as  a  function  of  sample  size  are  given 
in  Figure  37.  Confidence  factors,  fqQ,  for  a  90-percent  confidence 
lower  bound  on  the  scale  parameter  b  as  a  function  of  sample  size 
are  given  in  Figure  38. 

Using  Figure  37,  a  90-percent  confidence  lower  limit  for  the 
shape  parameter  a  is  obtained  by  multiplying  a  by  the  abscissa 
value  of  Figure  37  corresponding  to  the  proper  number  of  specimens. 
A  90-percent  confidence  l^iwer  limit  for  the  scale  parameter  b  is 
obtained  by  multiplying  b  by  the  abscissa  value,  corresponding 
to  the  proper  number  of  specimens  in  Figure  38,  to  the  1/^  power. 


7.3.2  Estimation  of  Scale  Parameter  When 
the  Shape  Parameter  is  Known 


In  this  section,  a  two-parameter  Weibull  distribution  with 
known  shape  parameter  a  and  unknown  scale  parameter  b  is 
assumed  to  be  the  underlying  distribution.  In  this  case,  the 

x  a 

cumulative  distribution  is  given  by  F(x|a,b)  =  1  -  exp  -  (-g)  . 

The  maximum-likelihood  estimator  ^  ,  of  the  scale  parameter  b, 

obtained  from  Equation  (30)  is  given  by 


1 

a 


A 

D 
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(38) 


where  x.,  i=l,2,...,n,  represent  a  random  sample  from  the  above 
We ibu.ll  attribution. 

A 

b  3 

It  is  shown  in  Reference  3  that  the  quantity  2n(^)  has  a 

Chi-squared  distribution  with  2n  degrees  of  freedom.  If  an  inter¬ 
val  estimate,  by  ,  of  b  is  required  which  satisfies  the  equation 

Pr  (fcy  <  b)  =  r  ,  (39) 

* 

the  by  is  defined  using  the  Chi-squared  distribution  by 


2 

where  Xy  (2n)  is  the  y-fractile  of  the  Chi-squared  variate  with 
2n  degrees  of  freedom. 

2 

Confidence  factors,  f^^  =  (2n/Xy  (2n)),  for  a  90-percent  con¬ 
fidence  lower  bound  on  the  scale  parameter  b  as  a  function  of 
sample  size  are  given  in  Figure  39.  The  90-percent  confidence^ lower 
limit  for  the  scale  parameter  b  is  obtained  by  multiplying  b 
by  the  abscissa  value  (corresponding  to  the  proper  number  of  speci¬ 
mens)  in  Figure  39  to  the  1/a  power. 


7.4  SAFE-LIFE  ESTIMATES 


In  this  section,  a  procedure  is  stated  for  obtaining  safe- 
life  estimates  as  a  function  of  desired  reliability  and  fleet  size. 
The  strategy  proposed  in  Reference  3  for  safe-life  predictions  is 
the  basis  for  this  approach,  whic-  will  allow  for  the  estimation 
of  the  wing  tension  splice  joints  safe  operational  life. 

7.4.1  Individual  Components 


The  reliability  for  an  individual  fleet  member  is  developed 
for  various  levels  of  safe-life.  This  reliability  is  defined  as 
the  probability  that  service  failure  will  occur  after  some  speci¬ 
fied  safe-life  for  that  individual  member.  These  estimates  are 
based  on  the  estimated  service-life  distribution.  Once  a  service- 
life  distribution  for  a  component  is  decided  upon,  then  for  fixed 
levels  of  reliability  the  safe-life  of  the  component  can  be 


determined.  In  Figure  41,  the  safe-life  is  the  abscissa  value  of 
life  that  corresponds  to  one  minus  the  desired  reliability  on  the 
ordinate  probability-of-failure  coordinate. 


7.4.2  Safe-Life  Estimation  for  a  Fleet 


In  this  step,  estimates  are  obtained  for  the  reliability  of  a 
fleet  of  structural  components.  The  results  allow  for  the  specifi¬ 
cation  of  the  first  and/or  second  and/or  third,  etc.,  failure  safe- 
life  probabilities. 

If  the  reliability  of  an  individual  component  is  R  ,  then  the 
probability  of  less  than  n  failures,  Rr  ,  in  a  fleet  consist¬ 
ing  of  N  components  is  given  by  ’ 


£ 


(?)  (l-R)1  RN_i! 


(41) 


where  (^)  = 

vi 


N1 


1  (N-l)l 


and  Rn  ^  is  the  reliability  of  the  n-th 


weakest  in  a  fleet  of  N.  Using  Equation  (41),  safe-life  estimates 
for  a  fleet  are  obtained  in  the  same  manner  as  for  individual  com¬ 


ponents.  The  relationship  between  R  and  R 


n,N 


for  various  values 


of  N  and  n  is  illustrated  in  Figure  42.  Safe-life  as  a  func¬ 
tion  of  desired  fleet  reliability  is  illustrated  in  Figure  43. 
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SECTION  VIII 

CONCLUSIONS  AND  RECOMMENDATIONS 

8.1  INTRODUCTION 


The  objectives  stated  in  Section  I  were  accomplished  during 
this  investigation.  While  some  additional  questions  were  answered, 
others  were  also  raised.  The  conclusions  and  the  recommendations 
for  further  work  are  summarized  in  this  section. 


8.2  CONCLUSIONS 


The  conclusions  of  this  investigation  are  summarized  below: 

1.  Miner's  rule  is  inadequate  for  predicting  the  fatigue 
life  of  bonded  composite  joints  under  random  loading. 

The  predicted  fatigue  life,  based  on  Miner's  rule, 
was  about  three  times  longer  than  that  obtained  by 
test. 

2.  The  fatigue  life  characterization  of  bonded  composite 
joints  should  be  based  on  a  random  load  history  and  a 
compatible  environment.  The  load  history  effects  are 
not  fully  understood  for  composites,  but  there  is  no 
justification  for  constant  amplitude  fatigue  testing 
if  the  random  fatigue-environment  test  philosophy  is 
adopted. 

3.  A  scaling  effect  does  exist  between  different  specimen 
sizes  subjected  to  the  same  random  stress  history.  The 
scaling  factor  increases  with  lifetime.  The  value  needed 
for  design  is  governed  by  the  fatigue  life  requirements 
rather  than  the  static  strength  requirements.  The  scal¬ 
ing  factor  is  significant  for  composite  bonded  joints 
and  must  be  carefully  evaluated  to  obtain  realistic 
design  allowables.  The  results  of  this  investigation 
indicate  that  the  scaling  factor  is  a  definable  quantity 
and  confirms  the  feasibility  of  using  small-scale  speci¬ 
mens  to  characterize  the  static  and  fatigue  life  of  full- 
size  specimens.  The  specimen  scale  sizes  must  be  care¬ 
fully  uefined  to  achieve  a  scale  compatibility  between 
different  specimens.  More  statistical  data  is  needed  to 


determine  the  number  of  specimens  required  to  define 
the  scaling  factors  and  to  perform  a  credible  relia¬ 
bility  analysis. 

4.  Three  modes  of  failure  were  observed,  i.e.  bond,  lami¬ 
nate,  and  complex  (bond  and  laminate).  The  laminate 
mode  of  failure  was  observed  for  several  static  and 
short  life-time  residual  strength  tests.  A  bond  failure 
mode  was  observed  for  some  of  the  fatigue  tests.  Due 

to  the  mixture  of  failure  modes,  no  definite  conclusions 
were  reached  on  the  failure  process  or  the  type  of 
failure  mode  expected  for  a  particular  fatigue  test. 

5.  The  boron- to- titanium  adhesive  bonded  double  scarf 
joint  type  used  for  this  study  exhibited  a  relatively 
short  fatigue  life  and  considerable  scatter  in  the  test 
results.  In  an  actual  design  situation,  a  much  larger 
factor  of  safety  (or  reduced  stress  level)  would  be 
required  than  that  used  for  the  design  of  the  1/5-  and 
1/2-scale  specimens  to  achieve  a  satisfactory  service 
life. 

6.  Statistical  inferences  for  a  single  boron-to-titanium 
adhesive  bonded  double  scarf  joint  can  be  made  using 
the  test  results  of  the  double  scarf  joint  specimen 
for  this  program.  The  relationship  between  the  statis¬ 
tical  distribution  of  a  single  joint  and  the  distribu¬ 
tion  of  the  least-of-two  joints  (such  as  the  double 
scarf  joint  specimen  used  for  this  study)  are  discussed 
in  Appendix  VII.  If  the  underlying  statistical  dis¬ 
tribution  of  failure  is  Weibull,  then  the  distribution 
of  both  the  single  joint  and  least-of-two  joints  have 
the  same  shape  parameter  a  ,  but  the  scale  parameter 
for  the  single  joint  is  given  by  2*-'a  fi  where  (3  is 
the  scale  parameter  of  the  least-of-two  joints. 

7.  The  residual  strength-versus-lifetime  wearout  curve  for 
the  1/5-scale  specimens  did  not  exhibit  the  same  shape 
as  the  constant  amplitude  S-N  type  curve.  The  wearout 
curve  (residual  strength  versus  lifetime)  was  concave 
downward  with  slight  initial  residual  strength  loss 
followed  by  rapidly  degrading  residual  strength.  The 
S-N  curve  was  convex  downward,  yielding  increasing 
life  with  decreasing  stress  amplitude.  Positive  cor¬ 
relation  was  achieved  with  the  initial  flaw  field  - 
kinetic  crack  growth  based  wearout  mode. 
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|  8.  The  number  of  specimens  required  to  define  the  factor  of 

safety  must  be  consistent  with  the  expected  confidence 
:  levei  in  the  statistical  parameters.  Since  the  coefficient 

\  of  variation  increases  for  the  Weibull  distribution  as 

?  a  decreases,  careful  consideration  must  be  given  to  the 

sample  size  and  its  impact  on  cost  and  the  required  factor 
;  of  safety. 

9.  The  random  load  history  simulation  used  in  this  program 
preserves  the  composite  cumulative  exceedance  statistics 
and  the  assumed  waveform.  Different  random  load  simu¬ 
lations  are  possible  for  the  same  load  history;  all  of 
these  satisfy  the  cumulative  exceedance  statistics.  The 
question  arises  as  to  which  random  load  simulation  pro¬ 
cedure  should  be  used  for  conducting  random  fatigue 
tests.  Even  though  all  random  load  simulations  satisfy 
the  composite  exceedance  statistics  for  the  same  load 
history  will  the  same  fatigue  lines  be  obtained?  The 
effect  of  different  random  load  simulations  on  the 
fatigue  life  should  be  studied,  and  if  there  are  sig¬ 
nificant  differences,  a  standard  simulation  procedure 
should  be  developed  for  the  aerospace  community. 


8.3  RECOMMENDATIONS 


Recommendations  for  further  research  are  summarized  below. 

The  RMS  endurance  curves  are  described  in  Reference  45. 

1.  Determine  the  character  of  an  RMS-stress  endurance  curve 
for  at  least  two  joint  design  concepts  to  define  an  al¬ 
ternate  joint  design  process.  The  RMS-stress  endurance 
function  should  be  correlated  with  a  residual  strength 
model . 

2.  Investigate  design  concepts  that  either  negate  the  scale 
effect  or  reduce  the  static  variability  of  bonded  joints. 
Determine  if  the  concepts  result  in  improved  fatigue 
reliability. 


3,  Develop  an  integrated  random  load  and  environment  test 
capability  and  assess  the  effects  of  realistic  service 
environment  on  joint  performance. 


4.  Alter  the  joint  design  criteria  and  test  strategy  and 
establish  design  stress  versus  scale  factors  that  will 
yield  satisfactory  fatigue  performance. 


5.  Begin  a  rational  investigation  of  the  wearout  and  life¬ 
time  characteristics  of  bolted  joints. 
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Four  of  the  pertinent  elements  in  the  preparation  of  the  ran¬ 
dom  load  history  simulation  procedure  are  discussed  in  detail  in 
this  section.  Included  are  Press’  method  for  decomposing  exceed¬ 
ance  curves  into  a  series  with  a  Gaussian  distribution  basis, 
power  spectral  analysis,  random  number  mapping,  and  load  sequence 
generation. 


PRESS'  METHOD  FOR  ANALYZING  EXCEEDANCE  CURVES 


The  exceedance  curves  of  paragraph  4.2.2  can  be  characterized 
by  RMS  levels,  by  time  spent  at  each  RMS  level  (P^),  and  by  clip¬ 
ping  ratios.  This  information  can  be  used  to  develop  a  simulated 
random  load  history  and  a  load-by-load  basis.  Press  and  his  co¬ 
workers  (Reference  29)  developed  a  method  for  decomposing  exceed¬ 
ance  curves  into  a  series  with  a  Gaussian  basis.  This  method  is 
described  below  using  Figures  44  and  45. 

1.  Construct  a  straight  line  that  is  tangent  to  the  tail  of 
the  exceedance  curve  depicted  in  Figure  44.  Thi  line 
characterizes  the  first  RMS  level  to  be  used  for  approxi¬ 
mating  the  observed  distribution.  Additional  straight 
lines,  or  RMS  levels,  may  be  required  to  realistically 
approximate  the  observed  distribution.  The  value  can  be 
determined  from  the  slope  of  the  straight  line  constructed. 
The  number  of  straight  lines  required  depends  upon  the 
shape  of  the  exceedance  curve  and  the  accuracy  desired. 

2.  Construct  a  second  straight  line  characterizing  the  second 
RMS  level,  which  intersects  the  first  straight  line  at  a 
point  where  F(X^)  is  one-half  the  actual  F(X£A)  for  the 
same  X^  value  and  is  tangent,  at  the  upper  end  of  the 
line,  to  the  exceedance  curve.  The  intersection  point 
between  the  two  lines  can  be  determined  by  comparing  the 
ordinate  to  the  first  straight  line,  for  a  given  X£,  to 
the  corresponding  ordinate  to  the  exceedance  curve. 

Note:  F/X^  =  1/2  F(XiA) 
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Additional  straight  lines,  if  required,  can  be  constructed 
in  a  similar  manner  using  the  last  line  c  nstructed  and 
the  exceedance  curve  segment  above  the  point  where  the 
last  line  was  tangent  to  the  curve.  The  third  line  inter¬ 
sects  the  second  line  at  point  i  +  1  as  shown  in  Figure 
44. 

Note:  F(X.+1)  =  1/2  F(X£+1  ) 

Only  two  or  three  RMS  levels  are  usually  required  to  ap¬ 
proximate  the  actual  exceedance  curve.  The  RMS  level, 
characterized  by  each  line,  is  maximum  for  the  first 
line  developed  and  the  RMS  value  decreases  for  successive 
lines.  The  sum  of  the  ordinates  to  each  RMS  level  (char¬ 
acterizing  line)  for  a  given  value  is  equal  to  the 
ordinate  of  the  exceedance  curve  at  the  same  X^  (abscissa) 
value  (reference  Figure  45),  i.e., 

F(Xt)  =  F(XX)  +  F(X2)  +  F(X3). 

Compute  tr  values  from  the  slopes  of  the  straight 
lines  representing  each  RMS  level.  Using  the  notation 
of  Figure  44,  the  slope,  m,  of  line  i  can  be  determined 
using  two  arbitrarily  selected  coordinate  points  on  the 
line,  i.e., 


m  =  Slope 
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The  equation  for  aRMS^  follows: 


u 

RMS. 
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Compute  the  time  spent  at  each  RMS  level  (Pj).  N0Pi  is 
the  zero  intercept  value  for  each  RMS  level  and  can  be 
determined  using  similar  triangles  as  depicted  in  Figure 
44.  The  following  expressions  are  obtained. 


In  {F(xle)/F(Xlr)| 
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If  X^  is  equal  to  zero,  then 
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In  general, 
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The  term  N0P-[  describes  the  total  number  of  load  exceed¬ 
ances  from  the  mean  for  each  RMS  level.  The  portion  of 
time  spent  in  a  given  RMS  level  is  computed  by  dividing 
the  total  number  of  exceedances  from  the  mean  (for  a  given 
mission  segment  type  exceedance  curve)  into  the  total  num¬ 
ber  of  exceedances  for  a  given  RMS  level. 


POWER  SPECTRAL  ANALYSIS 


The  power  spectral  analysis  is  an  exact  method  for  determin¬ 
ing  the  random  load  history  waveform  for  a  stationary  random  pro¬ 
cess  in  which  the  power  spectral  density  (PSD)  is  constant  over 
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the  frequency  band  of  interest.  For  other  random  force/time  his¬ 
tories,  as  shown  in  Figure  46,  it  is  possible  to  approximate  the 
autocorrelation  function  &  (r)  and  the  power  spectral  density 
function  4>  (w) ,  i.e., 

T 

*Kt)  =  t^oo  2 r  f  y<t)y(t+7)dt  (A2) 


and 

*(“)  - 1 


00 
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where  y(t)  =  applicable  load  amplitude  or  load  factor.  Rice 
(Reference  46)  developed  the  following  expressions  for  the  number 
of  positive  slope  crossings  of  the  mean  load  level  per  unit  time 
(Nq)  and  the  number  of  peaks  per  unit  time  (N  ) 
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The  terms  N  ,  N  and  N  /N  define  the  waveform  content  of  the  ran- 
o  p  o  p 

dom  load  signal. 

To  develop  the  waveform  information  desired,  the  PSD  content 
must  be  matched  to  an  actual  service  history.  Available  flight 
recorder  data  for  fighter-type  aircraft  were  surveyed,  but  no 
suitable  waveform  data  was  found. 
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'’’he  load  spectrum  used  for  the  test  tape  is  based  on  an  F-lll 
A/E/D  mission  analysis  (Reference  47).  The  data  reduction  for  this 
analysis  was  based  on  the  following  assumptions: 

1.  The  peak  loads  y(t)  must  be  at  least  ig  above  or  below 
the  lg  flight  level  value. 

2.  The  load/ time  trace  must  have  a  rise  and  fall  of  at  least 
lg  and  50  percent  of  the  incremental  peak  value  (the  in¬ 
cremental  peak  value  is  equal  to  the  peak  value  minus  lg) . 

All  waveform  data  is  filtered  out  using  these  assumptions. 

Due  to  the  lack  of  waveform  data,  a  value  was  assumed  for  the 
irregularity  factor  (N0/Np);  consequently,  a  power  spectral  analy¬ 
sis  was  not  required.  Nevertheless,  a  power  spectral  analysis  is 
theoretically  feasible  for  computing  N0,  and  Nc/Np  values  when 
the  PSD  content  is  known  for  the  service  history.  For  this  reason 
the  framework  for  a  power  spectral  analysis  is  retained  for  the 
computer  program  of  Reference  26, 


RANDOM  NUMBER  MAPPING 


The  Rice-Bendat-Kowalewski  functions,  the  corresponding  cumu¬ 
lative  distribution  function,  and  the  random  number  mapping  pro¬ 
cedure  are  discussed  in  this  section.  Random  number  mapping  is 
required  to  generate  random  delta  loads  that  satisfy  a  prescribed 
statistical  distribution  for  the  population. 

Rice-Bendat-Kowalewski  Function 


The  random  loads  generated  must  fit  a  frequency  distribution 
compatible  with  the  possible  load  occurrences  for  the  population. 

An  aircraft  experiences  both  positive  and  negative  load  exceedances 
during  operation.  There  are  generally  many  more  positive  load  ex¬ 
ceedances  than  negative  ones.  A  frequency  distribution  function 
is  needed  which  describes  a  population  of  positive  and  negative 
loads  in  a  prescribed  proportion.  The  Rice-Bendat-Kowalewski 
function  (References  46,  48,  and  49  respectively)  satisfies  this 
requirement.  This  function  combines  the  Gaussian  and  the  Rayleigh 
distributions.  When  N0/Np  is  equal  to  0.0,  W(Z)  reduces  to  the 
Gaussian  distribution  and  when  N0/Np  is  equal  to  1.0,  W(Z)  reduces 
to  the  Rayleigh  distribution.  For  Intermediate  N0/Np  values,  W(Z) 
is  a  mixed  mode  distribution.  Since  the  Rice-Bendat-Kowalewski 
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function  preserves  the  irregularity  factor,  N0/Np,  for  incremental 
load  changes,  this  function  was  used  in  the  random  load  history 
simulation. 

Rice-Bendat-Kowalewski  Function  (Page  171  of  Reference  50)  is 
given  below 


2  2 

-Z  /2K,  - 

1  -Z  /2 

+  /M  Ze  [1-P(Z/K2>] 

\Np/ 
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Z/K2 


(Error  function,  y  being 
a  dummy  variable) 


W(Z)  =  Probability  density  function 

Z  =  Normalized  variable  for  the  density  function 

N  =  Number  of  positive  slope  crossings  of  the  mean  per 
unit  time 

Np  =  Number  of  positive  peaks  per  second. 

N^Np  =  Irregularity  factor 


Cumulative  Distribution 


A  cumulative  distribution  curve  |P(Z)  versus  Z\  is  generated 
by  summing  incremental  areas  under  the  normalized  W(Z)  curve. 
Normalization  involves  adjusting  the  ordinates  tc  the  W(z)  curve 
so  that  the  area  under  the  W(Z)  curve,  between  the  positive  and 
negative  clip  limits,  (Z+  and  Z~  respectively)  is  equal  to  1.0. 
P(Z)  is  described  in  Figure  47  in  terms  of  W(Z). 
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The  cumulative  distribution  curve,  P(Z),  is  needed  to  deter¬ 
mine  random  values  of  Z  compatible  with  the  distribution  for  the 
population.  Acceptable  values  for  Z,  which  fit  the  P(Z)  distribu¬ 
tion,  are  determined  using  the  mapping  procedure  described  in  the 
following  paragraph. 


Mapping  Procedure 


Each  random  load  in  the  service-history  simulation  is  based 
on  a  random  number.  A  mapping  procedure  is  used  to  "fit"  the 
random  numbers  drawn  to  a  given  statistical  distribution;  i.e., 
the  Rice-Bendat-Kowalewski  function.  Random  numbers  are  used  to 
characterize  the  RMS  values  for  each  RMS  level  of  each  mission 
segment  type  and  each  mission. 

The  random  number  mapping  procedure  used  in  the  BY4  computer 
program  (Reference  26)  is  described  below.  This  procedure  is  used 
to  generate  loads  on  a  load-by-load  basis. 

1.  Compute  400  values  of  the  Rice-Bendat-Kowalewski  function, 
(Z),  for  a  given  N0/Np  value.  Assume  a  clipping  ratio  of 
+6.0.  The  resulting  curve  is  depicted  in  Figure  48. 

2.  Generate  a  normalized  cumulative  distribution  curve  using 
the  results  above.  Since  the  area  under  W(Z)  between  Z 
equal  +6.0  is  less  than  1.0  the  incremental  areas  are 
normalized  to  make  the  total  area  under  W(Z)  equal  to  1.0. 
Such  a  curve  is  shown  in  Figure  49. 

3.  Define  the  clipping  ratios  for  each  RMS  level  for  both 
the  positive  and  negative  spectra  and  for  all  segments. 
Assume  symmetrical  clipping  about  ZZ  equal  to  zero. 
Determine  for  each  clipping  ratio  the  corresponding  seg¬ 
ment  of  the  cumulative  distribution  curve  for  which  the 
ratio  applies.  The  cumulative  distribution  curve  between 
points  i  and  j  depicts  an  applicable  curve  segment. 

Points  i  and  j  on  the  curve  are  defined  by  (ZZ  . ,  YY  . ) 
and  (ZZ ^ ^ ,  YY^),  respectively,  1  1 

4.  Each  mission  segment  has  a  positive  and  negative  (spectra) 
clipping  value  for  each  RMS  level.  Clipping  is  assumed 
symmetric  for  both  the  positive  and  negative  spectra. 

Two  vertical  lines  are  drawn  in  Figure  50  to  represent 
the  largest  clipping  ratio  for  a  given  RMS  level.  A  box 
is  formed  between  the  lines  k  -  j,  j  -  ZZgj ,  ZZgj  -  ZZg^, 
and  ZZr  .  -  k,  and  since  this  box  characterizes  the  largest 

•'j 
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Figure  48  Clipped  Rice- Benda t-Kowa lews ki  Function 
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Figure  49  Normalized  Cumulative  Distribution  Curve 
for  Random  Mapping 
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clipping  ratio,  it  is  called  the  "Big  Box."  In  a  similar 
manner,  a  smaller  box  can  be  formed  within  the  "Big  Box" 
to  characterize  the  smaller  clipping  ratio.  This  box  is 
called  the  "Little  Box."  Random  numbers  that  fit  both 
boxes  must  be  selected. 


A  sequence  of  pseudo-random  numbers  between  0.0  and  10,000 
are  generated  using  the  programed  procedure  of  Reference 
51  and  each  number  is  divided  by  10,000  to  scale  it  to 
the  YY  axis.  Two  modified  random  numbers  (ZZgn  and  ZZlii) 
are  generated  for  each  pseudo-random  number  in  the  sequence. 
The  first  scaled  random  number  in  the  sequence  is  resealed 
to  fit  the  "Big  Box"  and  the  first  acceptable  random  number 
for  this  box,  ZZgn,  is  obtained  from  the  library  data  gen¬ 
erated  by  the  program.  ZZjjn  corresponds  to  the  point  YYgn 
on  the  cumulative  distribution  curve  between  segment  i  -  j. 
YYgn  is  defined  by  the  following  equation. 


YYd  =  YY  .  +  (YY_. 
Bn  Bi  v  Bj 


-  YY  . )  *  RN 
Bi  n 


where 

YY  =  Scaled  random  number  applicable  to  curve 
segment  i  -  j  ("Big  Box") 

YYg^ ,  YYg .  =  Ordinate  of  the  cumulative  distribution 

J  curve  at  ZZgi  and  ZZgj  respectively  (largest 
clipping  ratio  for  given  RMS  level) 

RN  =  Scaled  random  number,  n,  from  the  sequence, 
(n  =  1  for  first  random  number). 


An  acceptable  random  number  for  the  "Little  Box"  (ZZ^n) 
is  determined  using  YYgn  and  the  following  equation, 


YYLn  =  YYT .  +  (YYt  .  -  YY_  . )  *  YYb 
Lj  v  Lj  Li  Bn 


where 


Scaled  random  number  applicable  to  "Little 
Box. " 


Ordinate  of  the  cumulative  distribution  curve 
at  ZZn  and  ZZ^j  respectively  (smallest  clip¬ 
ping  ratio  for  given  RMS  level). 


The  technique  described  is  symbolized  in  Figure  50. 

134 


Load-by-Load  Generation 


The  load-by-load  generation  procedure  of  Reference  26  is 
described  in  this  section.  This  procedure  is  based  on  the  idea 
implemented  by  Dr.  J.  C.  Kalpin  of  the  Air  Force  Materials  Labora¬ 
tory  (Reference  1). 


The  essential  steps  in  the  procedure  are 

1.  Draw  a  random  number  and  fit  to  cumulative  distribution 

curve  (fitted  R.N.). 

J 

2.  Compute  the  delta  load  measured  from  the  mean  using, 

APj  »  (-l^*1  (fitted  R.N.) j  *  (crRMS)i 

Note:  j  =  1  for  first  load,  j  =  2  for  second  load,  etc. 
i  =  RMS  level 


3.  Compute  the  load  peak  value  using 


P.  =  P  + 
j  mean  4 


4. 


Note: 


The  sign  of 


can  be  positive  or  negative 


Increment  the  time  by  1/2  N0  (sec)  and  generate  the  next 
load.  A  haversine  waveform  is  used  to  connect  successive 
loads.  This  process  continues  until  the  required  number 
of  loads  have  been  generated  for  the  service  history 
simulation.  This  procedure  is  characterized  in  Figure  51. 


The  loading  signal  starts  with  a  positive  excursion  from  the 
mean.  A  negative  excursion  follows  each  positive  one  until  the 
time  spent  in  the  positive  or  negative  RMS  level  has  been  consumed 
(Figure  52).  Usually  more  positive  load  excursions  occur  than 
negative  ones;  therefore,  more  positive  loads  must  be  generated 
than  negative  ones  to  preserve  the  exceedance  statistics. 


The  time  to  be  spent  in  each  RMS  level  (positive  and  negative) 
for  each  mission  segment  type  is  computed  from  the  corresponding 
exceedance  curve.  After  each  load  is  generated,  the  time  is  incre¬ 
mented  and  the  next  load  is  generated.  The  BY4  computer  program 
(Reference  26)  keeps  track  of  the  time  spent  in  the  positive  and 
negative  RMS  levels  and  preserves  the  distribution  of  positive 
and  negative  load  excursions  in  a  RMS  burst.  A  RMS  burst  is  the 
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load  time  history  simulation  for  a  given  positive  and  negative 
RMS  level  (reference  Figure  53). 

The  BY4  program  is  written  so  that  the  load  signal  will  return 
to  the  mean  after  each  positive  load  excursion  when  all  negative 
load  excursions  for  a  given  RMS  level  have  been  generated  or  vice 
versa.  The  time  is  incremented  for  a  forced  return  to  the  mean; 
however,  this  time  is  not  counted  toward  the  total  time  required 
to  generate  the  positive  or  negative  loads  for  a  given  RMS  level 
(reference  Figure  54),  A  forced  return  to  the  mean  is  not  included 
in  the  total  load  count  for  the  simulation. 

The  procedure  developed  to  support  this  contract  activity  is 
a  pilot  procedure  only.  The  difficulties  encountered  will  require 
further  exploration. 

The  method  of  creating  the  delta  load  increment  leads  to  an 
error  in  input  to  output  N0/Np.  A  polynomial  correction  function 
was  used  to  remove  this  error. 

A  capability  for  generation  of  an  asymmetric  history  for  only 
the  case  of  No/Np  =  1.0  was  developed  and  checked. 
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APPENDIX  II 


H  RP  AGE  PROBLEM 

STRAIN  SURVEYS 


Strain  surveys  were  conducted  using  L/5-scale  specimens.  These 
surveys  were  performed  to  (1)  determine  the  feasibility  and/or 
need  for  removing  or  reducing  the  warpage  effect  in  the  initial  25 
1/5-scale  specimen  and  (2)  evaluate  the  significance  of  the  warpage 
effect  on  the  test  results.  The  strain  surveys  conducted  are  sum¬ 
marized  in  Table  XXXIV. 

The  tension  and  compression  strain  surveys  were  performed  to 
determine  the  feasibility  and/or  need  for  reducing  the  warpage  ef¬ 
fects.  The  static  strain  surveys  were  conducted  to  evaluate  the 
effect  and  significance  of  induced  bending  stresses  due  to  speci¬ 
men  warpage.  A  strain  survey  was  also  performed  on  a  specimen 
that  had  been  subjected  to  a  prior  constant  amplitude  load  history. 
The  purpose  of  this  survey  was  to  evaluate  the  significance  of 
specimen  warpage  after  cyclic  loading. 

Five  typical  1/5-scale  specimens  were  selected  from  the  initial 
25  built.  The  dimensional  variations  are  shown  in  Table  IX.  Speci¬ 
men  K994568  was  used  for  the  tension  and  compression  strain  surveys 
and  specimen  K994566  was  used  for  the  constant  amplitude  strain 
surveys.  Specimens  K012439,  K012441,  and  K994588  were  used  for  the 
static  strain  surveys.  Four  strain  gauges  were  mounted  back-to- 
back  on  each  specimen  as  shown  in  Figure  55.  The  strain  survey 
descriptions  and  results  are  given  in  the  following  sections. 

Tension  Strain  Survey 


The  tension  strain  surveys  were  performed  using  a  1/5-scale 
specimen  mounted  in  fixture  6  (Figure  56).  Strain  surveys  were 
conducted  ?ith  and  without  a  center  support  or  jam  bar.  Figure 
57  shows  the  test  set-up  for  the  center  support  case  without  jam 
bar.  The  center  support  is  made  up  of  two  long  channels  (welded 
to  the  fixture),  two  short  channels,  two  teflon  rub  strips,  two 
support  angles,  nuts,  and  bolts.  A  Teflon  rub  strip  is  attached 
to  each  short  channel  with  countersunk  bolts.  The  specimen  is 


mounted  between  the  two  short  channels  with  the  Teflon  rub  strips 
next  to  tne  specimen.  Four  bolts  tie  the  two  short  channels  to¬ 
gether  as  shown  in  Figure  56.  Washers  were  used  with  these  bolts 


Table  XXXIV  SUMMARY  OF  STRAIN  SURVEYS  CONDUCTED 


Survey  I.D. 

Descrl ption 

Tension  Checkout' 

(a)  0  to  60K  to  0  w5 thout 
center  support 

(b)  0  tc>  20"  to  0  with  center 
support 

(c)  0  t'->  /.OK;  jam  bar  instal¬ 
led  for  equa Jibing 
strains;  20K  to  60K  to  0. 

Compression  Checkout 

(a)  0  to  20K  without  C-clamps; 
clamps  installed  &  load 
reduced  to  -10K;  clamps 
readjusted  to  equalize 
strains;  -10K  to  60K  to 
-10K  to  0. 

(b)  0  to  -10K  without  jam  bar; 
added  jam  bar  to  equalize 
strains,  ~10K  to  0. 

Static 

Static  tested  three  specimen; 
read  strains  for  10K  incre¬ 
ments;  loaded  to  failure. 

Constant  Amplitude 

(a)  Cycled  specimen  for  260 
cycles  using  5K  to  50K 
and  then  conducted 
strain  survey. 

(b)  Same  as  (a)  except  3300 

cycles  were  put  on 

before  strain  survey. 

' 
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Figure  55  Strain  Gage  Locations  for  Figure  56  Set  Up  for  Tension  and  Compression 

1/5-Scale  Specimens  Strain  Surveys  with  Center  Support 


Strain  x  10  (in. /in. 


to  separate  the  channels.  With  this  arrangement  there  is  no 
mechanical  attachment  between  the  specimen  and  the  center  support. 
The  center  support  is  removed  by  taking  out  the  two  short  channels. 
The  jam  bar  is  composed  of  a  1"  diameter  steel  rod  with  threaded 
ends  and  mating  nuts.  When  the  jam  bar  is  used  some  of  the  washers 
are  removed  between  the  two  short  channels  so  that  the  teflon  can 
be  forced  against  the  specimen  by  the  jam  bar.  The  jam  bar  is 
placed  between  the  left  hand  vertical  support  of  fixture  6  and  the 
face  of  the  corresponding  short  channel.  Pressure  is  put  on  the 
specimen  by  torquing  the  jam  bar  nuts.  The  strain  gauges  were  read 
for  a  given  loading  and  the  jam  bar  was  adjusted  to  equalize  the 
back-to-back  strain  gauge  readings. 


Figure  57  Tension  Strain  Survey;  without  Center  Support 


Using  the  various  support  options  the  specimen  was  loaded  in 
increments  and  the  strains  were  read.  To  preclude  lamina  yielding 
a  maximum  load  of  60K  was  used.  The  strain  gauge  readings  for  the 
cases  studied  are  summarized  in  Table  XXXV  and  XXXVI.  The  strains 
are  plotted  in  Figure  57  for  the  case  without  the  center  support. 
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Table  XXXV  TENSION  STRAIN  SURVEY  USING  1/5-SCALE 
SPECIMEN  WITHOUT  CENTER  SUPPORT 


Koto:  Specimen  K994563  used  for  survey 
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Table  XXXVI  TENSION  STRAIN  SURVEY:  WITH  CENTER  SUPPORT; 
WITH  CENTER  SUPPORT  AND  JAM  BAR 
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Compression  Strain  Survey 


A  compression  strain  survey  was  performed  using  specimen 
K994568  and  test  fixture  6  (Figure  57).  The  same  basic  apparatus 
used  for  the  tension  strain  survey  was  used  for  this  strain  survey. 
In  addition,  two  C-c lamps  were  used  to  clamp  the  two  short  channels 
against  the  specimen.  The  specimen  was  mounted  in  fixture  6  using 
a  given  support  option.  C-clamps  and  a  jam  bar  were  used  to  adjust 
the  strain  gauge  readings.  Loads  were  applied  to  the  specimen  in 
increments  and  the  strain  gauges  were  read.  The  maximum  tension 
and  compression  load  used  was  60K  and  -10K, respectively .  The  strain 
gauge  data  generated  is  shown  in  Table  XXXVII  and  XXXVIII. 


Static  Test  Strain  Survey 


Static  strain  surveys  were  conducted  using  specimen  K012439, 
K012441,  and  K994588  and  a  Baldwin-Tate-Emery  Universal  test 
machine.  The  specimens  were  mounted  as  shown  in  Figure  26.  Load 
was  applied  to  the  specimen  in  10K  increments  until  failure 
occurred.  The  strain  gauge  readings  at  each  load  level  and  the 
failing  loads  are  summarized  in  Table  XXXIX.  The  load  versus 
strain  is  plotted  for  applicable  specimens  in  Figures  58,  59,  and 


60. 


Constant  Amplituaa  Strain  Survey 


This  strain  survey  was  performed  to  evaluate  the  effect  of 
specimen  warpage  on  a  specimen  with  a  prior  constant  amplitude 
load  history.  Specimen  K994566  was  mounted  in  fixture  5  as  shown 
in  Figure  22, and  the  strains  were  read  under  zero  load.  Two  strain 
surveys  were  performed.  Two-hundred  and  sixty  cycles  of  5K  to  50K 
constant  amplitude  loading  were  first  applied  to  the  specimen. 

Static  loads  were  applied  in  increments  and  the  strains  were  read. 
After  reaching  a  maximum  load  of  50K  the  load  was  reduced  in  incre¬ 
ments  to  zero  and  the  strains  were  read  at  each  load  level.  An 
additional  3040  cycles  of  the  5K  to  50K  constant  amplitude  loading 
was  applied  to  the  specimen  and  a  new  strain  survey  was  performed 
as  before.  The  strain  gauge  data  generated  is  summarized  in  Table 
XL. 
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Table  XXXVII  COMPRESSION  STRAIN  SURVEY:  WITH  CENTER  SUPPORT; 
WITH  CENTER  SUPPORT  AND  C-CLAMPS 
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Table  XXXIX  STRAIN  SURVEY  USING  STATIC  TEST  SPECIMENS  (1/5-SCALE) 
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DATA  EVALUATION 


The  initial  set  of  1/5-scale  specimens  built  were  warped  in 
varying  degrees.  The  warpage  problem  was  eliminated  in  future 
specimens  by  modifying  the  manufacturing  procedure.  An  evaluation 
was  made  in  this  section  to  determine  the  significance  of  speci¬ 
men  warpage  on  the  test  results. 


Use  of  Warped  One-Fifth  Scale  Specimen 


The  effects  of  specimen  warpage  and  their  impact  on  the  ten¬ 
sion  and  compression  test  results  are  evaluated  in  this  section. 

A  discussion  is  given  for  justifying  the  use  of  selected  warped 
specimens . 

Evaluation  of  Strain  Surveys 

Feasibility  of  Negating  Warpage  in  Specimen  Constraint  - 
Tension  and  compression  strain  surveys  were  performed  using  various 
constraints  for  negating  the  effects  of  specimen  warpage.  It  was 
found  that  constraints  could  be  imposed  on  the  specimen  for  equal¬ 
izing  the  strains  rather  well  for  a  fixed  loading  condition  but 
the  effect  could  not  be  retained  for  different  load  levels  (ref. 
Tables  XXXV  through  XXXVIII).  It  was  concluded  that  the  effect  of 
specimen  warpage  could  not  be  satisfactorily  negated  by  imposing 
specimen  constraints. 

Significance  of  Induced  Bending  -  The  significance  of  induced 
bending,  due  to  specimen  warpage,  is  evaluated  in  this  section  using 
the  tension  and  static  strain  survey  results.  The  effect  of  the 
induced  bending  stress  is  evaluated  by  determining  the  contribution 
of  the  bending  stress  to  the  total  stress  in  the  specimen.  This 
contribution  is  evaluated  in  terms  of  a  stress  increase  factor,  F  . 
defined  below. 

2  € 
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i  j 
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where 


F  . .  =  stress  increase  factor  for  bending  moment 
mij 


e 

i 
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€ 
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=  back-to-back  strain  gauge 


readings  for  gauge  i  and  j. 


€ 

max 


e.  when  e.  >  e. 
i  1  J 


e .  when  e .  s  e . 

J  J  ^  i 

Fmij  is  based  on  the  assumption  that  the  extreme  fiber  stresses  on 
thespecimen  cross  section  (°i,  crj)  are  composed  of  a  uniform 
stress,  (Ta,  and  a  bending  stress,  anc^  the  tw0  can  be  determined 
by  superposition  as  shown  in  Figure  61. 


The  stress  increase  i.'aclor  for  induced  bending  is  computed  by 
dividing  the  average  stress,  <ra,  into  the  maximum  extreme  fiber 
stress  (tfj_  or  °j).  Since  the  increase  factor,  Fm£j ,  is  non- 
dimensional  strains  are  used  in  place  of  stresses. 
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Table  XLI  STRESS  INCREASE  FACTORS  FOR  INDUCED  BENDING 
MOMENT  BASED  ON  TENSION  STRAIN  SURVEY 
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;ote:  Loads  and  strains  taken  from  Table  XXXVI 
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HOTE:  Load c  and  strains  taken  from  Table 


Stress  increase  factors,  Fm^j ,  are  summarized  in  Tables  XLI 
and  XLII  for  the  tension  and  static  strain  surveys,  respectively. 
The  results  in  both  cases  show  that  Fm^-  generally  decreases  as 
the  load  increases.  This  indicates  that  specimens  tend  to  self 
align  with  increasing  load.  Fmjj  ranges  are  shown  in  Table  XLIII 
for  the  tension  and  static  strain  survey  specimens. 


i 

} 


Table  XLIII  Fmij  ^anSes  for  Tension  and  Static 
Strain  Survey  Specimen 


Specimen 

Strain  Survey 

Fmij  Range  (KIPS) 

K994568 

Tension 

1.074(5)  to  1.043(60) 

K012439 

Static 

1.0026(10)  to  1.0211(100) 

K0 12441 

Static 

1.0222(10)  to  1.0106(100) 

K994588 

. 

Static 

1.1257(10)  to  1.0506(80) 

Effect  of  Cyclic  Loading  -  A  strain  survey  was  performed  using 
a  specimen  subjected  to  a  prior  loading  history.  The  results  shown 
in  Table  XL  indicate  that  the  strain  readings  taken  after  260 
cycles  and  3300  cycles,  respective  ly,  were  essentially  unchanged. 

At  the  50K  load  level,  the  F^j  factor  was  approximately  1 .04  after 
260  and  3300  cycles,  respectively. 

Rejection  Criteria 

The  dimensional  variations  for  the  twenty-four  specimens  in 
Table  IX  were  surveyed  and  all  specimens  with  a  surface  measurement 
equal  to  or  greater  than  0.100  inch  at  any  of  the  eight  locations 
shown  in  Figure  11  were  rejected.  Using  this  criteria  seven  speci¬ 
mens  were  rejected  ^nd  seventeen  were  accepted  ,  . .  Table  IX). 

The  dial  gauge  sun  ace  measurements  at  points  1  c:. rough  8  were 
averaged  for  the  seventeen  specimens  accepted.  The  surface  dimen¬ 
sions  for  the  five  specimens  used  for  the  strain  surveys  a/e  com¬ 
pared  against  the  mean  surface  dimensions  in  Table  XLIV. 
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The  surface  dimensions  for  the  specimens  used  for  the  strain  sur¬ 
veys  represents  reasonably  well  the  average  specimen  depicted  in 
Table  IX. 

Justification  for  Specimen  Usage 

Specimen  tend  to  straighten  out  as  the  load  is  increased; 
therefore,  the  effects  of  warpage  at  failure  are  negligible.  The 
static  test  results  of  Table  XV  indicate  that  the  warpage  is  a 
second  order  effect  because  the  rejected  specimens  exhibited 
static  failing  loads  of  the  same  order  of  magnitude  as  the  ac¬ 
cepted  specimen.  Also, there  was  no  distinction  in  the  failure 
modes . 

Twelve  warped  specimens  were  used  to  conduct  the  constant 
ampiit  de  fatigue-to-failure  tests.  Seven  specimens  were  used  for 
the  tension-tension  tests,  and  five  were  used  for  the  tension- 
compression  tests.  The  tension-tension  constant  amplitude  test 
results  for  the  warped  and  unwarped  specimens  (built  after  the 
warpage  problem  was  corrected)  all  fell  in  a  reasonable  scatter 
band  (Table  XIV). 

The  tension-compression  constant  amplitude  test  results  for 
the  warped  specimen  did  not  differ  very  much  with  the  unwarped 
specimen  results.  For  example,  specimen  K900442  failed  at  2959 
cycles  of  a  -10K  to  50K  loading  and  the  unwarped  specimens  failed 
between  1790  and  3193  cycles.  The  scatter  in  the  data  does  not 
appear  to  be  unreasonable.  The  tension-compression  tests  were 
performed  to  (j.)  characterize  the  fatigue  life  of  the  specimens 
under  this  type  of  loading  and  (2)  permit  a  comparison  of  the 
relative  damaging  effects  of  a  tension-tension  constant  amplitude 
loading  versus  a  tension-compression  loading.  The  test  results 
based  on  the  warped  and  unwarped  specimens  satisfy  these  objectives 

Two  fatigue-to-failure  tests  were  performed  using  warped  speci 
mens  (K012443  and  K012464).  The  lifetimes  for  these  specimens  fall 
well  in  the  range  of  values  for  all  specimen  tested  (Table  XVII). 


157,158 


APPENDIX  III 


GLOSSARY  OF  TECHNICAL  TERMS 


The  principal  technical  terms  used  in  this  report  are  defined 
in  this  section. 

1.  A/A:  Air-to-air  cycle 

2.  A/G:  Air-to-ground  cycle 

3.  Angular  frequency  (omega):  The  angle  9  through  which  the 
body  rotates  divided  by  the  time  required  to  perform  this 
rotation,  or 


w  =  e/t 

4.  Asymmetric  spectrum:  An  asymmetric  spectrum,  as  used  in  this 
report,  is  one  which  has  different  positive  and  negative  ex¬ 
ceedances  for  given  mission  segment  types. 


5.  Autocorrelation  function:  A  function  describing  the  depen¬ 
dence  of  data  at  one  time  t  to  that  of  another  time  t  +  t 


*(0 


lim  ^ 
T— *  2T 


y(t)y't  +  r  )dt 


where : 


'H  ?  ) 

T 


T  = 


autocorrelation  function 
variable  time  delay  (sec) 
period  (sec) 
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6. 


Average:  The  sum  of  the  results  divided  by  the  number  of 
trials ,  i . e . , 


Average  = 


n 


£ 


n 


where  =  result  for  trial  i 
n  =  number  of  trials 


7.  Brittle  material:  A  material  for  which  strength  is  uniquely 
determined  by  the  distribution  of  flaws  and  toughness.  Rupture 
is  caused  by  raising  the  stress  intensity  at  the  greatest  flaw 
to  the  critical  value. 

8.  Broad  band  process:  A  broad  band  process  is  a  stationary 
process  whose  mean  square  spectral  density  has  significant 
value  over  a  band  or  range  of  frequencies  which  is  roughly 
the  same  order  of  magnitude  as  the  center  frequency  of 

the  band.  A  wide  range  ;  frequencies  appears  in  representa¬ 
tive  sample  functions  such  a  process.  (Reference  52 
pp.  17-18.) 

In  a  typical  aircraft  broad  band  random  load  versus  time 
signal  NQ/Np  is  less  than  one  (see  below). 


9.  Clipping  ratio:  The  maximum  value  of  the  load  signal  ampli¬ 
tude  (actually  achieved  or  a  truncated  value)  divided  by  the 
RMS  value.  Clipping  values  are  determined  for  each  RMS  level 
required  to  characterize  the  exceedance  curves  for  each  mis¬ 
sion  segment  type.  The  tail  of  the  exceedance  curve  is  trun¬ 
cated  as  shown  below.  The  truncation  level  can  usually  be 
determined  intuitively  because  the  tail  of  the  exceedance 
curve  tends  to  wander  at  the  higher  load  level  because  the 
frequency  of  occurrence  is  much  less  and  there  is  a  general 
lack  of  statistics. 
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10.  Closed  loop  test  system:  A  testing  system  which  uses  servo/ 
hydraulic  test  equipment  to  trace  a  set  of  desired  load  time 
excursions  with  fidelity. 


11.  Coefficient  of  variation  (CV) :  Defined  as  the  RMS  (<r) 
value  divided  b>  the  mean  or 


where  <r  -  RMS  value 

x  =  mean  value  of  x  for  sample 

The  coefficient  of  variation  can  be  used  to  estimate  the 
tail  of  the  distribution  curve  from  a  limited  sample. 

12.  Confidence  level:  A  measure  of  the  confidence  in  the  results 
with  respect  to  sample  size. 

13.  Cumulative  distribution  function:  The  cumulative  probability 
distribution  function,  F(x) ,  defines  the  probability  that 

a  random  variable  will  be  equal  to  or  less  than  x. 
Symbolically, 

F(x)  =  probability  |xn  =  X  j  . 
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F(x)  defines  the  area  under  the  probability  density  curve 
between  selected  limits.  The  range  of  F(x)  is  from  0  to  1.0. 
Refer  to  the  sketches  below  for  further  description. 


14.  Ensemble:  The  entire  load  history  in  time  and  amplitude  of 
a  random  load  signal  which  has  an  infinite  number  of  load 
samples  and  given  time  range. 

15.  Ergodic  process:  An  ergodic  process  is  one  in  which  the 
timewise  behavior  of  y(t)  over  a  long  period  of  time  will 
exhibit  the  same  average  properties  as  corresponding 
ensemble  averages  at  various  fixed  times. 

16.  Exceedance  curve:  A  curve  which  characterizes  statistically 
the  load  history  for  a  given  mission  segment  type.  An  ex¬ 
ceedance  curve  defines  the  number  of  times  per  mission  a 
given  load  exceeds  the  mean  value  for  a  giv'en  mission  segment 
type.  The  exceedance  curve  for  this  report  is  plotted  as 
the  (  AB.M.)2  (abscissa,  linear)  versus  the  cumulative 
number  of  exceedances  prr  mission  (ordinate,  log  plot). 
Exceedance  curves  are  generated  for  positive  and  negative 
spectra  using  measured  flight  data.  The  complete  fluctuat¬ 
ing  load  history  for  each  mission  segment  type  of  the 
mission  profile  may  be  characterized  by  an  appropriate 
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exceedance  curve.  A  typical  exceedance  cure  is  shown  below. 


Ln 

Cumulative 
Exceedances 
Per  Mission 


(A  B.M.)2 

A  composite  cumulative  exceedance  curve  characterizes  statis¬ 
tically  the  load  history  for  all  mission  segment  types  com¬ 
promising  a  mission.  The  composite  cumulative  exceedance 
curve  is  obtained  by  summing  the  exceedances  for  the  positive 
and  negative  spectra  (separately).  A  composite  cumulative 
exceedance  curve  can  also  be  generated  for  all  missions  in 
a  service  life  in  a  similar  manner. 

17.  Fatigue  failure:  A  structural  failure  (either  rupture  or  a 
dangerous  length  subinternal  flaw)  which  necessitates  repair 
or  replacement  of  the  vehicle  or  part.  The  failure  is  caused 
by  the  combined  mechanical/chemical  normal  service  environments. 

18.  Flaw  distribution:  The  statistical  distributions  of  flaw  sizes 
in  a  structure. 

19.  Frequences  of  occurrence:  The  number  of  times  a  given  event 
occurs  in  the  population. 

20.  GAG:  Ground-air-ground  cycle 

21.  Gaussian  Random  Process:  "A  random  process  1 3^  (t)}  is  said  to 
be  a  Gaussian  random  process  if,  for  every  set  of  fixed  times 
Jtn},  the  random  variables  ^R^n)  follow  a  multi-dimensional 
normal  distribution."  (Reference  52,  Page  91.) 

22.  Gust  loads:  Fluctuating  loads  imposed  on  the  aircraft  by 
air  turbulence. 

23.  Irregularity  factor  (N0/Np) :  Defines  the  number  of  positive 
slope  crossings  of  the  mean  per  positive  peaks  per  unit  time. 

N0  and  Np  can  theoretically  be  calculated  from  the  power  spec¬ 
tral  density.  Due  to  lack  of  aircraft  wave  form  data  an 
assumed  value  of  0.85  was  used  for  this  report.  NQ/Np  must 
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be  defined  to  evaluate  the  Rice- Bendat-Kowalewski  function. 
NOTE:  N0/Np  =1.0  (Rayleigh) 

No/Np  =  0.0  (Gaussian) 

24.  Level  crossings:  Refers  to  the  number  of  load  peaks  above 
(positive  spectra)  or  below  (negative  spectra)  a  given  load 
level  in  a  given  time  span.  The  number  of  load  peaks  above 
or  below  selected  load  levels  are  counted  for  the  simulated 
random  load  history  and  the  results  are  used  to  check  agree¬ 
ment  with  the  actual  load  history  exceedances. 


y(t) 


25.  Lifetime  distribution:  A  probability  frequency  distribution 
which  characterizes  the  time  to  specimen  failure.  Fatigue 
lifetimes  are  determined  experimentally  by  continuously  test¬ 
ing  specimen  under  a  simulated  random  load  history  until 
failure  occurs.  A  Weibull  distribution  fitted  to  the  test 
results  provides  a  physical  description  of  the  fatigue  pro¬ 
cess.  A  typical  lifetime  distribution  curve  is  depicted 
below. 


26.  Moan  load:  The  reference  load  level  from  which  the  load 

amplitude  is  measured.  Tn  most  aircraft  applications  a  Ig 
load  is  considered  to  be  the  reference  load  level.  The  mean 
load  level  at  Ig  may  vary  for  each  condition.  For  a  symmetric 
history  the  reference  is  the  algebraic  mean. 
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27.  Mission  profile:  Characterizes  the  activities  or  maneuvers 
the  airplane  is  expected  to  encounter  during  a  mission  and 
the  time  spent  in  each  activity.  The  activities  are  charac¬ 
terized  by  mission  segment  types.  A  typical  mission  profile 
is  depicted  below: 


®KD-+ — @ — H— iKH-dH 


Mission  Segment 


Time 


28.  Mission  Segment:  A  mission  segment  defines  the  type  of  activ¬ 
ity  the  aircraft  encounters  during  a  mission.  A  mission  is 
composed  of  different  mission  segment  types.  The  order  of 
the  segments  and  the  time  spent  in  each  segment  type  dis¬ 
tinguishes  the  missions  flown.  Typical  mission  segment  types 
are  climb,  cruise,  TFR,  A/G,  A/A,  GAG,  and  T&G.  The  load 
history  for  each  mission  segment  type  is  characterized  by 
exceedance  curves. 

29.  Narrow  band  process:  "A  narrow  band  process  is  a  stationary 
random  process  whose  mean  square  spectral  density  has  signifi¬ 
cant  values  only  in  a  band  or  range  of  frequencies  whose  width 
is  small  compared  with  the  magnitude  of  the  center  frequency 
of  the  band.  Only  a  narrow  range  of  frequencies  appear  in 
representative  samples  of  such  a  process."  (Reference  13, 
page  153).  A  typical  narrow  band  signal  is  shown  on  page 

17  of  Reference  52. 

1  65 
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30.  Negative  spectra:  Characterizes  the  load  exceedances  below 
the  mean. 

31.  Negative  load  excursion:  A  load  fluctuation  below  the  mean 


32. 


unit  time. 


33.  Np :  The  number  of  positive  peaks  (mathematical  maximum)  per 
unit  time. 


y(0 


34. 

35. 


P^:  Define  the  portion  of  time  spent  at  each  RMS  level. 

Population:  A  collection  of  all  possible  observations  for 
an  event.  For  example,  all  airplanes  in  the  fleet  compromise 
a  population.  A  portion  of  the  airplanes  in  the  fleet  is  a 
sample . 


36. 


Positive  spectra; 
the  mean. 


Characterizes  the  load  exceedances  above 


37.  Positive  load  excursion:  A  load  fluctuation  above  the  mean. 
(Reference  sketch  under  negative  load  excursion.) 


I  (»(> 


38.  Power  spectral  analysis:  This  is  an  analysis  in  which  N0 

and  Np  are  computed  from  the  power  spectral  density  function 
(PSD).  The  expressions  below  were  developed  by  Rice 
(Reference  42) . 


J  (w)  do) 

M  _1_  o _ 

N°  =  2  -p 

Z7r  J  0(0))  dO) 


(Reference  29,  page  26) 


w  =  angular  frequency  (radian/sec) 
f  =  cyclic  frequency  (Hz) 

39.  Power  spectral  density  function  (PSD):  A  frequency  analysis 
of  the  mean  square  of  the  random  variable.  The  PSD  charac¬ 
terizes  the  amount  of  energy  in  a  frequency  band.  The  total 
energy  is  obtained  by  adding  the  spectral  density  for  each 
frequency  increment  along  the  spectrum.  Using  the  PSD  the 
distribution  of  level  crossings,  peaks,  and  ranges  may  be 
computed  for  simulating  a  loading  history  (Reference  13, 
page  157  and  Reference  52,  page  22). 

40.  Probability:  The  ratio  of  the  expected  number  of  occurrences 
of  an  event  divided  by  the  total  number  of  possible  occur¬ 
rences  . 


41.  Probability  density  function:  A  function  which  characterizes 
the  frequency  of  occurrence  of  possible  events  in  a  popula¬ 
tion.  The  probability  density  functions  used  in  this  report 
are  summarized  below: 
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(References  46,  49  and  53). 


42.  Random  load:  A  load  from  the  total  load  population  which 
fits  a  prescribed  statistical  distribution. 


43.  Random  load  history:  An  ensemble  of  random  loads  which 
characterizes  the  load  history  for  a  single  mission,  for 
several  missions,  or  for  a  complete  service  life. 


44.  Random  number  mapping:  This  is  a  procedure  for  generating 
random  loads  by  fitting  random  numbers  drawn  from  a  rec¬ 
tangular  distribution  to  a  given  cumulative  distribution. 

For  this  report  the  Rice-Bendat-Kowalewski  Cumulative  dis¬ 
tribution  function  is  used.  A  detailed  description  of  the 
procedure  is  given  in  Appendix  I. 

45.  Ranking  data:  The  probability  of  survival,  P(x)  for  a  given 
set  of  test  results  can  be  expressed  using  the  equation 
below: 


P(x) 


r 


a 


n+l 


where:  ra=  rank  of  the  data 

n  =  number  of  tests  in  the  data  set 

The  data  is  ranked  in  descending  order  of  magnitude.  The 
largest  value  has  a  rank  of  1.0,  the  next  largest  has  a  rank 
of  2.0,  and  so  forth  with  the  smallest  value  having  a  rank 
of  n.  For  this  report  the  probability  of  survival,  based 
on  the  ranking  procedure,  is  compared  with  the  theoretical 
values  obtained  from  a  Weibull  distribution  fitted  to  the 
test  results.  This  comparison  shows  how  well  the  Weibull 
distribution  fits  the  observed  probabilities  of  survival 
based  on  the  ranking  procedure.  (Reference  42  ^  page  38  and 
Reference  37) 
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46. 


Reliability  function:  The  reliability  function,  R(x),  defines 
the  probability  that  a  value  of  the  random  variable  is  greater 
than  x.  Symbolically, 

R(x)  -  probability  |xn  >  xj^  =  1  -  F(x) 

v;here:  F(x)  -  cumulative  distribution  function 


R(x) 


47.  Reliability  goal:  A  reliability  standard  to  which  an  element, 
a  component,  or  assembly  is  to  meet. 

48.  Residual  strength:  The  static  strength  remaining  in  a  member 
after  the  member  has  been  subjected  to  a  given  load  history, 
environment,  or  other  damage. 

49.  RMS  burst:  A  segment  of  the  random  load  signal  which  charac¬ 
terizes  all  the  loads  for  a  single  RMS  level  (positive  and 
negative  spectra)  for  a  given  mission  segment  and  mission 
segment  type.  Two  or  three  RMS  levels  are  usually  required 
to  characterize  the  exceedance  curve  for  a  given  mission  seg¬ 
ment  type. 


Random  Load  Segment 


50.  RMS  level:  The  exceedance  curves  for  various  mission  segment 
types  are  characterized  by  different  RMS  values  according  to 
the  method  of  Press  (Reference  29).  Each  RMS  value  required 


\ 
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to  characterize  the  exceedance  curve  is  called  an  RMS  value. 
Two  or  more  RMS  levels  are  usually  required  to  characterize 
the  exceedance  curve  for  a  given  mission  segment  type. 


1. 


Root-mean-square  ( RMS ) :  The  RMS  value  is  a  statistical  term 
which  characterizes  the  load  deviations  from  the  mean.  For 
this  report  the  mean  is  at  the  lg  load  level.  The  mean  square 
value  is  given  by: 


2 

(T 

RMS 


lim 
T—  * 


T 

|  f  y2(t)dt 
o 


where:  T  =•  observation  time 


y(t)  =  ordinate  to  the  load  signal 


The  RMS  value  is  equal  to  the  positive  square  root  of  the 
mean  square  value,  i.e,, 


RMS  values  are  determined  from  exceedance  curves  for  given 
RMS  levels  by  the  Method  of  Press  (29).  The  RMS  value  can 
also  be  estimated  in  terms  of  discrete  loads  measured  from 
the  mean  (x)  using  the  equation  below. 


(r 

RMS 


-V 


r 


n 

E 

i=l 


-  2 

(xi  -  x) 


n-1 


where:  n  =  total  number  of  loads 
x^  =  value  of  load  i 
x  =  mean  value  for  all  measurements 
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Safety  factor:  In  structural  reliability  analysis  the  safety 
factor  is  the  ratio  of  the  mean  strength  divided  by  the  mean 
load.  Due  to  the  variation  in  load  and  strength,  the  safety 
factor  must  be  selected  with  respect  to  an  acceptable  proba¬ 
bility  of  failure.  The  load  and  strength  distributions  are 
depicted  below.  Failures  can  be  represented  by  the  shaded 
area  at  the  tails  of  the  distributions. 


Relative 

Frequency 


53.  Sample:  A  sample  is  a  particular  set  of  observations  out  of 
the  set  of  all  possible  observations  in  the  population.  One 
airplane  out  of  all  airplanes  in  the  fleet  constitutes  a 
sample  from  the  population. 

54.  Scaling  factor:  A  factor  which  allows  one  to  characterize 
the  fatigue  life  of  full-size  components  under  a  random  load 
history  using  small-scale  specimens.  Small  specimens  generally 
have  longer  ratigue  lives  than  their  full-size  counterparts 
subjected  to  the  same  stress  history.  The  scale  factor  is 
used  co  scale  small  specimen  fatigue  results  to  full-size 
structures . 

55.  Scale  parameter  (  /3  ) :  A  characteristic  value  of  the  Weibull 
distribution  below  which  63 . 2%  of  the  values  for  the  random 
variable  are  expected  to  fall.  It  can  be  determined  by  fitting 
the  Weibull  distribution  to  the  test  data.  For  residual 
strength  distributions  /?  is  expressed  in  pounds  and  for 
lifetime  distributions,  (3  is  expressed  in  lifetime.  An 
equation  for  the  probability  of  survival  (Weibull)  is  given 
below: 


-(—  )a 

P(cr  )  =  e  j3  ' 


whore  !’(<»  )  =  probability  of  surviving  residual  strength, 
<r  =  residual  strength 
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/?  =  scale  parameter  (characteristic  value) 
a  -  shape  parameter 

56.  Scatter:  A  group  of  specimen  built  to  the  same  specifications 
and  tested  with  the  same  fatigue  loading  will  have  different 
fatigue  lives.  If  all  things  were  truly  equal  all  the  speci¬ 
men  should  have  the  same  fatigue  life.  Scatter  is  a  term 
used  to  describe  the  variation  in  fatigue  lives.  The  amount 
of  scatter  is  often  described  by  the  lifetime  RMS  value;  e„g., 
design  life  equals  (scatter  factor  x  test  life). 

57.  Service  life:  The  service  life  of  an  airplane  is  a  specified 
time  which  the  airplane  is  expected  to  perform  its  intended 
purpose.  Service  life  projections  are  continually  made  on 
airplanes  in  the  fleet  to  evaluate  the  safety  for  continuous 
usage.  Fighter  aircraft  are  designed  for  4000  hours  of  usage. 

A  typical  F-lll  aircraft  is  expected  to  fly  about  1334  missions 
during  its  service  life.  For  test  purposes  in  this  program  the 
random  load  history  is  simulated  for  one  service  life  and  is 
stored  on  digital  magnetic  tape. 

This  load  history  represents  one  lifetime  of  usage.  Four  or 
more  lifetimes  are  usually  required  to  satisfy  the  fatigue 
life  requirements;  i.e.  a  scatter  factor  of  4. 

58.  Servo  amplifier:  A  device  for  adjusting  the  servo  valve  out¬ 
put  to  make  it  consistent  with  the  input  load  requirement 
from  the  magnetic  tape. 

59.  Servo  valve:  A  valve  which  controls  the  flow  of  oil  for  driv¬ 
ing  the  hydraulic  loading  ram.  It  receives  a  feed  back  signal 
from  the  servo  amplifier. 

60.  Shape  parameter  (a):  A  nondimensional  parameter  of  the 
Weibull  distribution.  The  shape  of  the  resulting  Weibull 
distribution  curve  is  governed  by  this  parameter.  Also  refer 
to  scale  parameter  ( /3  ) . 

61.  Stationary  random  process:  A  process  in  which  the  statistical 
parameters  are  constant  with  respect  to  time. 

62.  Statistics:  The  scientific  collection  and  analysis  of  data, 
and  the  projection  of  estimates  therefrom.  Statistics  is 
not  an  end  in  itself;  rather  it  is  an  aid  to  judgement  in 
arriving  at  valid  conclusions,  testing  theories,  measuring 
phenomena,  discovering  relationships,  or  projecting  estimates 
under  different  conditions.  Statistics  is  only  a  tool  and  it 
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must  be  remembered  that  it  is  not  a  substitute  for  profes¬ 
sional  knowledge  and  experienced  judgment.  Sound  judgment, 
however,  cannot  be  made  with  confidence  without  the  use  of 
this  tool. 

63.  Stress  intensity  factor  (K) :  The  stress  intensity  factor  is 
a  parameter  used  in  fracture  mechanics  for  characterizing 
the  elastic  stresses  in  the  area  of  a  crack  tip.  K  depends 
upon  the  loading  mode  and  the  geometry  of  the  crack. 

64.  Structural  reliability  philosophy:  In  the  structural  relia¬ 
bility  approach  the  probability  of  structural  failure  is 
predicted  considering  the  variations  in  strength  and  applied 
loading.  Conventional  structural  design  methods,  on  the 
other  hand,  evaluate  the  integrity  of  the  structure  on  a 
margin  of  safety  basis.  Since  the  structural  reliability 
approach  statically  characterizes  the  variation  in  strength 
and  loading  it  permits  a  more  realistic  assessment  of  the 
structural  integrity.  (Reference  safety  factor  discussion). 

65.  Symmetric  clipping:  A  clipping  ratio  is  computed  for  each 
RMS  level  for  both  the  positive  and  negative  spectra  exceed¬ 
ance  curves  for  each  mission  segment  type.  This  clipping 
value  defines  the  limits  for  the  tail  to  the  probability 
density  curve  on  the  positive  side  of  the  mean  or  character¬ 
istic  value.  When  the  density  curve  on  the  negative  side  of 
the  mean  is  clipped  the  same  way  as  the  opposite  side  sym¬ 
metries  clipping  is  said  to  be  used  (refer  to  sketch  below). 


66.  Symmetric  spectrum:  A  symmetric  spectrum,  as  used  in  this 
report,  is  one  in  which  the  positive  and  negative  exceedances 
are  the  same  for  given  mission  segment  types. 

67.  TFR:  Terrain  following  radar 

68.  T&G:  Touch  and  go  landings 
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69.  TIP:  Transonic  Improvement  Program;  predessor  to  USAF  "TACT" 
program  being  directed  by  the  AFFDL. 

70.  Truncation  level:  The  number  of  large  load  occurrences  in 
the  load  history  population  are  infrequent  compared  to  the 
number  of  smaller  load  occurrences.  Unless  a  very  large 
sample  of  the  load  history  is  taken  there  is  generally  a  lack 
of  statistics  to  adequately  characterize  the  large  load 
occurrences.  In  such  cases  the  loads  are  cut  off  or  truncated 
at  some  reasonable  maximum  load  level.  In  this  report  the 

(  A  B.M.)2  loads  from  the  exceedance  curves  and  the  maximum 
positive  and  negative  load  excursion  from  the  mean  random 
load  signal  are  truncated.  Truncation  of  the  exceedance 
curves  is  discussed  under  clipping  ratios.  Truncation  of 
the  random  load  signal  is  depicted  below.  The  only  centered 
truncation  used  resulted  from  input  data  or  test  machine 
accuracy. 

y(t) 


71.  Variance  (cr^):  a  statistical  term  defined  as  follows: 


n 

2  £  <xi  “  *) 

<r  =  i=l 


n- 1 


where:  n  =  number  of  measurements 
Xj  =  value  of  measurement  i 
x  -  mean  value  for  all  measurements 


72.  Waveform:  T>*  a  random  process  the  waveform  content  charac 
terizes  the  load  peaks,  the  frequency  of  load  occurrence. 
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and  load  reversals.  The  wave  form  is  characterized  by  the 
irregularity  factor  (N0/Np). 

73.  Wearout  curve:  A  curve  in  which  the  residual  static  strength 
is  plotted  against  lifetime.  This  curve  graphically  charac¬ 
terizes  the  residual  strength  of  the  specimen  with  respect 

to  service  life. 

74.  White  noise:  A  term  used  to  indicate  that  the  power  spectral 
density  (PSD)  is  constant  over  the  frequency  range  of  interest. 

The  term  "White  noise"  was  coined  due  to  the  fact  that  all 
frequencies  contribute  to  the  random  process  in  a  manner 
similar  to  the  contribution  of  all  colors  to  white  light. 
(Reference  13,  page  157.) 
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APPENDIX  IV 


PHOTOGRAPHS  OF  1/5-SCALE 


SPECIMEN  FAILURES 


Specimen  photographs  in  this  section  are  grouped  according 
to  test  types  with  the  specimens  for  each  group  given  in  numeri¬ 
cal  order.  Joint  failure  modes  are  discussed  in  subsection  6,7. 


Figure  63  One-Fifth  Scale  Static  Specimen  K012441 
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Figure  64  One-Fifth  Scale  Constant  Amplitude  Fatigue  Specimen 
K994566 


Figure  65  Full  Size:  One-Fifth  Scale  Static  (KO'12441)  and 
Constant  Amplitude  Fatigue  (K994566)  Specimens 
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Figure  66  Ten  Times  Size:  One-Fifth  Scale  Static  (K012441)  and 
Constant  Amplitude  Fatigue  (K994566)  Specimens 


Figure  67  Twenty  Five  Times  Size:  One-Fifth  Scale  Static  (K012441) 
and  Constant  Amplitude  Fatigue  (K994566)  Specimens 
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Figure  69  Constant  Amplitude  (Tension-Tension)  Fatigue  Test 
One-Fifth  Scale  Specimen  K012477 


Figure  70  Constant  Amplitude  (Tension-Tension)  Fatigue  Test 
One-Fifth  Scale  Specimen  K012478 


Figure  71  Constant  Amplitude  (Tension-Tension)  Fatigue  Test 
One-Fifth  Scale  Specimen  KOI 2483 
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Figure  73  Constant  Amplitude  (Tension-Tension)  Fatigue  Test 
One-Fifth  Scale  Specimen  K994587 
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Figure  76  Constant  Amplitude  (Tension-Compression)  Fatigue  Test 
One-Fifth  Scale  Specimen  K012462 
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Figure  77  Constant  Amplitude  (Tension-Compression)  Fatigue  Test 
One-Fifth  Scale  Specimen  K012475 


Figure  78  C 
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Figure  79  C 
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Random  Loading  (17,  Lifetime)  Test:  One-Fifth  Scale 
Specimen  K900441 


Figure  80 


Figure  81  Random  Loading  (170  Lifetime)  Test:  One-Fifth  Scale 
Specimen  K900448 
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Figure  84  Random  Loading  (1%  Lifetime)  Test:  One-Fifth  Scale 
Specimen  K905016 


Figure  85  Random  Loading  (10%  Lifetime)  Test:  One-Fifth  Scale 
Specimen  K900441 
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Figure  92  Random  Loading  (50%  Lifetime)  Test:  One-Fifth  Scale 
Specimen  K904951 
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Figure  93  Random  Loading  (507  Lifetime)  Test:  One-Fifth  Scale 
Specimen  K904952 


Figure  96  Random  Loading  (50%  Lifetime)  Test: 

One-Fifth  Scale  Specimen  K904956 


Figure  97  Random  Loading  (50%  Lifetime)  Test: 

One-Fifth  Scale  Specimen  K905015 
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Figure  98  Random  Loading  (50%  Lifetime)  Test 
One-Fifth  Scale  Specimen  K905018 


Figure  99  Random  Loading  (507  Lifetime)  Test 
One-Fifth  Scale  Specimen  K905022 
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Figure  102  Fatigue- to-Failure  Tests:  One-Fifth  Scale 
Specimen  K012464 


Figure  103  Fatigue-to-Failure  Tests: 
Specimen  K012479 


One-Fifth  Scale 


Figure  105  Fatigue-to-Fai lure  Tests:  One-Fifth  Scale 
Specimen  KOI 2481 
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Figure  106  Fatigue- to- Failure  Tests:  One-Fifth  Scale 
Specimen  K0 12484 


Figure  107  Fatigue- to-Failure  Tests:  One-FiiV:>  Scale 
Specimen  K012485 
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Figure  108  Fatigue- to-Fai'lure  Tests:  One-Fifth  Scale 
Specimen  K90Q443 
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APPENDIX  V 

PHOTOGRAPHS  OF  1/2-SCALE 
SPECIMEN  FAILURES 


Photographs  of  1/2-scale  specimens  are  given  in  this  section 
for  the  10%  lifetime  tests  and  the  fatigue-to-failure  tests. 
Failure  modes  are  discussed  in  subsection  6.7. 
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Figure  110  Random  Loading  (10%  Lifetime)  Test: 

One-Half  Scale  Specimens  (Back  Side) 


Figure  111  Fatigue-to-Failure  Test:  One-Half  Scale 
Specimen  F504619  (Front  Side) 


Figure  112  Fatigue-to-Failure  Test:  One-Half  Scale 
Specimen  F504619 (Back  Side) 
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Figure  115  Fatigue-to-Failure  Test:  One-Half  Scale 
Specimen  F504622  (Front  Side) 
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Figure  116  Fatigue-to-Fai lure  Test:  One-Half  ScaL 
Specimen  F504622 (Back  Side) 
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APPENDIX  VI 


STATISTICAL  PROPERTIES  OF 
THE  WEIBULL  DISTRIBUTION 


The  Weibull  cumulative  probability  distribution  function  is 
represented  mathematically  as  follows. 


F(t|a,b)  =  1  -  exp( — (-—)  ), 


where  a  and  b  are  parameters  greater  than  zero  to  be  ’determined 
from  test  data,  and  where  t  is  the  random  variable  involved,  e.g., 
time  to  failure,  cycles  to  failure.  The  Weibull  probability  den¬ 
sity  function  is  obtained  by  differentiating  Equation  (44) 


f(t|a,b)  exp( — (—•)  ). 


This  function  is  plotted  in  Figures  119  and  120  for  various  values 
of  a  and  b.  The  failure  rate  function  defined  by 


p / f.  I  g  M  _  f  tt  |a,b;_  __ 
R(t|a,b)  1  _  F(t|a>  b) 


becomes 


R(t|a,b)  = 


Jl  (JL) 
b  vb  ' 


This  function  increases  with  time  for  most  fatigue  situations  since 
in  these  situations  the  parameter  a  is  greater  than  one.  Equation 
(47)  is  plotted  in  Figure  121  for  various  values  of  the  shape 
parameter  a. 

The  k-th  moment  of  the  distribution  is  given  by 


E(tk)  =  bk  l'(^-  +  1). 
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Weibull  Density  Function  -  Scale  Parameter  Fixed  to  One 


where  r  is  the  standard  gamma  function.  Thus,  the  mean  and 
variance  of  the  distribution  are 


e(c)  =  b  r<-r  +  i). 


Var(t)  -  b2  r£-  +  1) 


r2(T  +  1>] 


In  Equation  44  if  t  =  b  then  F(t=b|a,b)=l  -  exp(-l)  =  .632,  which 
is  constant  regardless  of  the  value  of  the  parameter  a.  Thus  t  =  b 
is  called  the  characteristic  value  of  the  random  variable.  The 
parameter  b  determines  the  location  of  the  distribution  and  the 
parameter  a  being  inversely  proportional  to  the  variance  of  the 
random  variable  t  (see  Equation  (50))  is  called  the  shape  parameter 
of  the  distribution.  The  Weibull  distribution  is  positively  skewed 
for  values  of  the  shape  parameter  a  less  than  3.57,  and  is  nega¬ 
tively  skewed  for  values  of  the  shape  parameter  greater  than  3.57. 

A  very  useful  parameter  of  the  We...  .11  distribution  is  the 
coefficient  of  variation  which  is  defined  by 


CV(t)= 

f  ^  E(t) 

Thus  using  equations  (49)  and  (50)  equation  (51)  becomes 


CV(a)  = 


nxj-  +  i)  -  r2(~7  +  i )11/2 


The  coefficient  of  variation  is  only  a  function  of  the  shape  param¬ 
eter  a  and  reflects  relative  variation  independent  of  the  underlying 
scale  parameter  of  the  distribution.  A  plot  of  Equation  (52)  is 
I  given  in  Figure  40. 


APPENDIX  VII 

LEAST-OF-TWO  STATISTICAL  ANALYSIS 


The  test  specimens  used  for  this  program  had  two  nominally 
identical  sides  (Figures  5  and  6).  Due  to  symmetry,  one  side  of 
the  specimen  fails  before  the  other  side  when  fatigue  tested. 

The  service  life  of  symmetric  specimens  may  be  thought  of  as  the 
time  to  first  failure  of  two  nominally  identical  specimens  or  the 
time  to  failure  is  the  least-of-two  (Reference  54). 

The  service  life  distribution  parameters  may  be  determined 
by  estimating  the  service  life  variability  parameter  from  small 
element  tests,  and  estimating  the  central  location  parameter  of 
the  service  life  distribution  from  full-scale  fatigue  tests. 

Thus  the  major  part  of  least-of-two  statistical  analysis  is  being 
able  to  relate  the  distribution  of  the  least-of-two  variate  and 
its  parameters  to  the  distribution  of  a  single  variate  and  its 
parameters . 

A  least-of-two  statistical  analysis  is  developed  in  the 
following  paragraphs  for  the  log-normal  and  Weibull  distributions. 


LOG-NORMAL  LEAST-OF-TWO  STATISTICAL  ANALYSIS 


Based  on  the  fact  that  the  logarithm  of  the  number  of  cycles 
to  failure  for  a  single  component  is  normally  distributed,  the 
probability  distribution  for  the  log  life  of  a  full  scale  symmetric 
component  made  up  of  two  nominally  identical  halves  is  the  least- 
of-two  normal  variates. 

Let  F(xj/x,cr^)  and  f(xbu.,<r  )  denote  the  cumulative  distri¬ 
bution  and  density  function  respectively  for  a  normally  distributed 
variable  with  mean  ^  and  variance  cr2.  Then  the  logarithm  of  the 
least-of-two  log-normal  variables  has  the  same  distribution  as  the 
least-of-two  normal  variables  which  each  have  distribution 
F(x|^,<r2) . 

If  y  denotes  the  least-of-two  normal  variables,  then  the 
probability  density,  g(y|M,or2),  of  y  is  given  in  terms  of 

F(x|/x,o-2)  ancj  f(xjAt,o-2)  by 

=  2f(y|/',<r2)  (1  -  F(yj/i,o-2)  )  (53) 


I 

* 

s 
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and  the  cumulative  distribution,  G(y| /^,o-2) }  of  y  is  given  by 
G(yU,cr2)  =  ^  g(yU,cr2)dy 


— OC 


=  2F(y|^,cr2)  -  F2(y|/t,cr2) 


(54, 


Equations  (53)  and  (54)  are  easily  derivable  from  the  Theory  of 
Order  Statistics  (see  Reference  55). 

The  mean  and  variance  of  y  are  derived  to  be 


E(y)  =  P- 


(?r) 


1/2 


(55) 


Var(y)  =  or2 (■-■■■  *-) 


(56) 


Equations  (55)  and  (56)  allow  estimates  of  the  mean  and  variance 
of  y  to  be  transformed  into  estimates  of  the  mean  and  variance 
of  the  single  normal  variable  and  vice  versa. 

With  the  properties  derived  here,  standard  statistical  tech¬ 
niques  can  be  applied  to  the  least-of-two  distribution  for  making 
inferences  regarding  service  life  characteristics  of  the  least- 
of-two  log-normally  distributed  variables. 


WE I BULL  LEAST-OF-TWO  STATISTICAL  ANALYSIS 


In  this  section  the  Weibull  distribution  with  the  shape  para¬ 
meter  a  and  scale  parameter  b  is  assumed  to  be  the  parent  dis¬ 
tribution.  The  cumulative  distribution  function  in  this  case  is 
given  by 


F(x|a,b)  =  1  -  exp(-(£)a)  (57) 

The  cumulative  distribution,  G(y|a,b),  of  the  least-of-two  Weibull 
variates  is  given  in  terms  of  F(x|a,b)  by 

G(y|avb)  =  2F(y|a,b)  -  F2(y|a,b)  (58) 
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Equation  (58)  is  easily  derivable  using  Theory  of  Order  Statistics 
(see  Reference  51).  Equation  (58)  can  be  rewritten  as 


F(y|  a,b)  =  F(y|a,b)  (2  -  F(y|a,b)  )  .  (59) 

Then  substituting  Equation  (57)  into  (59)  and  simplifying,  the 
distribution  of  y  becomes 

G(y|a,b)  =  1  -  exp(-  2(^)a)  ,  (60) 

which  is  a  Weibuil  distribution  with  shape  a  and  scale  parameter 
b(2~S).  That  is 

1 

G(y|a,b)  ^  F(y ja,b(2  a)  )  .  (61) 

Thus  the  only  difference  between  the  distribution  of  the 
least-of-two  Weibuil  variables  and  a  single  Weibuil  variable  is 
the  scale  parameter. 

Statistical  inference  about  the  least-of-two  Weibuil  variates 
can  be  based  on  existing  statistical  techniques  for  making  infer¬ 
ences  about  the  Weibuil  distribution. 
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